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CLOUD 


ABSTRACT 


The critical supercooling, required produce visible condensation 
pressure-limited adiabatic expansions presumably ion-free saturated vapor 
has been measured small glass cloud chamber. The transient gas temperatures 
were measured during expansions using carbon resistance thermometer 
Wheatstone bridge. oscilloscope was used place galvanometer, and its 
trace was photographed during expansions. Low-amplitude temperature oscil- 
lations the gas, which might not normally have been detected, were frequently 
observed the early stages this work. These oscillations either appeared spon- 
These oscillations were subsequently eliminated. Condensation thresholds were 
then measured using starting temperatures from 4.2 1.7° AT, dropped from 
about mdeg above the point less than mdeg below the point. 
The critical supersaturations required produce condensation were thus always 
less than 105%, The supersaturations were calculated without making the usual, 
but highly implausible, assumption that the expanding gas ideal. Assuming 
that the condensation nuclei are embryonic droplets arising accidentally from 
density fluctuations the supersaturated vapor, the critical droplet radius 
found about above 2.5° although the theoretical treatment here 


not rigorous. 


INTRODUCTION 


The physical processes involved the initiation condensation from 
supersaturated vapor are means clear. Most cloud chamber experiments 
have involved the use some condensable vapor mixed with another ‘carrier’ 
gas which does not participate the condensation. Both gases are always 
assumed ideal attempts theoretical analyses the processes. 

The authors undertook examine condensation single-component 
systein, using pure low temperatures. Furthermore, the measurements 
supercooling and consequent calculations supersaturations were made 
without using the totally implausible assumption that vapor the verge 
condensation behaves ideal gas. 


EXPERIMENTS 


2.1 Apparatus and Method 
The condensation cloud chamber used (Fig. consisted glass chamber 


about 110 cm’ volume, which had offset lower tip The chamber was 
received October 13, 1959. 
Contribution from the Department Physics, Royal Military College, Kingston, Ontario. 
Supported part the Defence Research Board Canada under Grant No. 9510-10. 
Project D44-95-10-10. 
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Fic. The apparatus. 


immersed liquid helium bath and connected room temperature 
5/8-in. outside diameter Inconel tube about meter long, and then through 
fast-acting Veeco solenoid vacuum valve (SV) large tank approxi- 
mately liters volume. Tube was soft-soldered into the brass top the 
chamber. The seal between and was effected means Teflon bushing 
which was machined press fit room temperature. light film 
vacuum grease was wiped the parts before each assembly. Upon cooling, the 
Teflon contracts much more than either brass glass that seal tight 
liquid helium was regularly, but not invariably, formed. bakelite spacer 
was used support the 56-ohm 1/2-w Allen Bradley carbon composition 
resistance thermometer Holes were drilled through permit easy gas 
flow past it. The leads were No. copper wire. 

The cotton batten plug was loosely stuffed into the lower end 
eliminate temperature oscillations the gas (see 2.2 later). 

The chamber was illuminated 6-volt microscope illuminator used with 
three heat-absorbing filters. The chamber was viewed visually through the 
slits the silvering the liquid air and helium Dewars, darkened room. 
The helium clouds formed could best seen when looking almost directly 
towards the light source. 

small amount liquid was condensed from external filling system, 
into the offset tip the cloud chamber, using He‘ gas which had been purified 


oa 
Liquid 
Helium — I 
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liquid-air-cooled charcoal trap. The initial saturated vapor pressure, 
the chamber was then read and the closed. The large tank was 
adjusted some pressure less than Then the was opened that 
expansion essentially constant final pressure (the tank pressure) occurred. 
After preset time, usually between 0.4 and 1.5 seconds, the was closed 
again. The cloud chamber was watched determine whether not fog 
cloud appeared the supersaturated vapor. The temperature the vapor was 
measured during expansions using the suspended carbon resistance thermo- 
meter one arm Wheatstone bridge. medium persistence DuMont 
Type 403 oscilloscope with d-c. sensitivity 100 per was used this 
bridge instead galvanometer. The drop temperature the vapor 
during expansion was recorded photographing the oscilloscope trace. 
the temperature fell the resistance rose, and the bridge unbalanced, pro- 
ducing deflection the trace the oscilloscope. This method following 
the cooling during expansions was developed Stachérska (1954, 1956). The 
power dissipated was kept below and our sensitivity was then 
between and per deflection the oscilloscope trace. 

Figure shows typical photograph the oscilloscope trace, taken during 
expansion from initial temperature 2.725° The sweep speed 
was 250 milliseconds per division. The camera shutter was opened for total 
about seconds showing parts three sweeps. Starting near the upper 
right, the bridge was balance. After about second, the was opened 
(at about 1.8 divisions from left) and the temperature fell about mdeg 
about 100 milliseconds, and then remained nearly constant. After about 1.7 


Fic. Temperature record typical expansion. 
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seconds the was closed and the temperature began rise towards the bath 
temperature. After further 2.7 seconds, the shutter was closed again. 


2.2 Temperature Oscillations the Gas 


(a) Spontaneous Oscillations 

When the apparatus was first set up, without the cotton batten plug (see 
Fig. place, low-amplitude temperature oscillations the gas frequently 
occurred. Temperature oscillations gas column, one end which 
room temperature, and the other low temperature, have been known 
exist for many years (Keesom 1942) and are usually referred 
The exact nature the phenomenon not clear (Taconis 
1949; Kramers 1949; ter Haar 1955). Audible resonances often are able 
transport impressive amount heat (up several watts) into cryostat. 
Such large amplitude oscillations are very easy detect when they occur 
because the sound and/or the increased boil-off rate the bath liquid. 

The low-amplitude temperature oscillations which observed might very 
easily have been overlooked without cathode-ray oscilloscope the bridge 
circuit the resistance thermometer the gas. Figure shows such 
low-amplitude resonance which occurred when the chamber was pressure 
21.5 3.24° (S.V.P. 255.5 mm). The temperature amplitude 
+22 mdeg frequency cycles/sec. the same temperature with 
the amplitude was reduced +12 mdeg, and above 

some runs temperature oscillations occurred even the saturated vapor 
pressure and below the point. 


Fic. Spontaneous temperature oscillation 
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Such oscillations, undetected apparatus similar geometry experi- 
ments measure the film thickness flow properties, might have very 
serious effects upon those measurements, because the extreme sensitivity 
such experiments small temperature differences. would easy install 
carbon thermometer circuit such have described, any film experi- 
ment. This could then act low-amplitude Taconis resonance detector. 


(b) Shock-excited Oscillations 

Even when had chosen initial conditions for expansions where tem- 
perature oscillations occurred, regularly found that the expansion itself 
could excite damped oscillation, shown Fig. taken 3.2° with the 
kept open for 0.3 second. The frequency this oscillation was about 


Fic. Temperature oscillation shock-excited expansion. 


cycles/sec. After the cotton batten plug was stuffed into the lower end 
the tube (see Fig. both the spontaneous and shock-excited temperature 
oscillations the gas were damped out and did not reappear any subsequent 
run.* All the later expansions were then similar that shown Fig. 


2.3 Results 

each seven temperatures, from 1.7 4.2° from expansions 
were performed. any such run the initial temperature was kept constant. 
The cooling was subsequently measured from enlargements the photo- 
graphs the C.R.O. trace. any temperature, large expansions produce 
persistent fogs helium from the supersaturated vapor, whereas small 
expansions produce condensation. Each series expansions was performed 


*We are grateful Dr. Jackson for this suggestion. 
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such way try determine the critical supercooling required 
produce visible condensation. this threshold supercooling more nearly 
realized, the condensation produced sometimes seemed only few large 
drops like rain, else brief flash fog. The measurements finally give two 
numbers for any initial temperature viz. maximum temperature drop 
giving visible condensation, and minimum temperature drop 
ATyes giving visible condensation. Presumably the true threshold lies between 
these two numbers. Figure graph these measurements showing A7,, 
circles and plus signs. 1.7° the point believed rather 
high because the vapor density low (1/50 its value 4.2° that the 
first condensation may not have been seen. The accuracy the other points 
probably better than mdeg. 


2.0 3.0 4.0 


Fic. Maximum temperature drop ATno which gives visible condensation, and 
minimum temperature drop which gives condensation, function the temperature 
before expansion saturated He‘ vapor. 


ANALYSIS 
Supersaturation Imperfect Gas 

follow Wilson (1897) defining supersaturation, the ratio the 
actual density the vapor immediately after expansion the saturation 
density the lower final temperature 7;. Only when the vapor ideal 
gas does this definition become equivalent where the actual 
pressure and the saturation vapor pressure over plane surface. 

wish find method calculating from the measured cooling 
without having assume that the expanding gas ideal. The initial conditions 
before expansion are: the saturated vapor pressure, the 
saturated molar volume, and After adiabatic expansion through the 
critical temperature interval have P;, and 


O + 
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wish evaluate the critical supersaturation 


where the saturated molar volume the final temperature. Now 


where the average rate change molar volume along the 
saturated vapor curve. And where the actual change 
molar volume during the adiabatic expansion. 

now calculate AV, for imperfect gas. First obtain relationship 
between and adiabatic change. The equation state any gas 


Then 


Cy oT Vv 


one Maxwell’s relations. Now for any imperfect gas the equation 
state may written volume virial expansion 


Then 
for imperfect gas (and for ideal gas also). 
And hence 


Thus the supersaturation 


This the expression desired. 
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3.2 Temperature Dependence the Critical Supersaturation 

order apply eq. (3.1), need detailed information the equation 
state the saturated vapor. (Cy/P) values were obtained graphical 
averaging the measurements van Itterbeek and Laet (1958). The values 
used were calculated from the refractive index data 
Edwards (1957) for vapor. The Lorenz—Lorentz law was used, 
taking the polarizability vapor, the liquid value (Edwards 1958) 
Thus, 


where the refractive index, was used calculate the saturated molar 
volumes from Edwards’ (1957) data. 

Table shows the resulting critical supersaturations taking the critical 
supercooling the average the and points Fig. The 


TABLE 


Condensation data for showing final temperatures, critical supercooling, experimental 
and theoretical critical supersaturations, and critical radii, and number molecules per 
droplet calculated two different methods 


4.181 1.018 1.077 104 61.7 


3.147 1.033 1.297 137 105 66.3 104 
2.671 53.6 1.047 1.549 135 2.21X105 
2.167 32.2 1.045 2.129 202 118. 
2.139 13.6 1.020 465 280. 2.03 
1.980 3.8 1.007 2.586 1570 108 950. 
1.705 8.5 1.022 3.593 579 392. 


results show that down about 2.2° lower temperature requires greater 
supersaturation produce condensation. Below about there seems 
drop the critical supersaturation required. Because the rapidly 
decreasing vapor density and therefore the added difficulty observing con- 
densation, the critical supercooling seems more likely have been over- 
estimated than underestimated these lowest temperatures. Such error, 
present, would then tend make the critical supersaturations the lowest 
temperatures appear too high, rather than too low. 

Volmer and Flood (1934) and later Farley (1952) have discussed the kinetics 
the condensation process and have predicted that the temperature variation 
the critical supersaturation should given 


3/2 


where constant, empirically equal 0.557 according Volmer and 
Flood. the surface tension the liquid, the absolute temperature, 
and the molar volume the liquid. 

only 0.7% less than 4.18° and still closer 
lower temperatures; may verified using the second and third virial 


EDWARDS AND WOODBURY: CONDENSATION Het VAPOR 343 


coefficients for gas. Table the values calculated from 
eq. (3.3) are shown (theory). clear that the predicted temperature 
variation far greater than observed experimentally. 

these calculations used the data Allen and Misener (1938) for 
and obtained from Edwards (1958) data liquid density. 


3.3 The Mechanism Condensation 

has long been recognized that condensation requires nucleation (see, for 
example, Wilson 1951). Dust particles water-saturated air are known 
very effective condensation nuclei. are also known effective— 
this fact leads the usefulness cloud chambers tracking ionizing particles. 
Even the absence foreign particles ions, however, condensation can 
take place, self-nucleation, what sometimes referred the ‘cloud 
limit’. 

The possibility that the condensation our experiments occurred foreign 
nuclei ions cannot entirely ruled out. However, care was taken 
avoid vacuum leaks which might produce solid air the walls the Inconel 
tube (see Fig. which might become dislodged during expansion and fall 
into the chamber. The cotton batten plug and the bakelite spacer would 
course tend stop any such particles before they could enter the chamber. 

experiments were undertaken determine the effect ionization the 
condensation thresholds. However, the high energy (46 ev) required form 
ion pair helium gas low temperatures (Hanauer 1959) and the 
reasonably high density the clouds produced our chamber are taken 
evidence that were studying self-nucleated condensation ion-free vapor. 
this case the condensation nuclei are presumed embryonic droplets 
neutral molecules arising accidentally from density fluctuations the vapor. 
The range sizes these embryonic droplets assumed larger for higher 
supersaturations, until finally some the embryonic droplets exceed critical 
size and then persist. When this size reached, the vapor then unstable 
equilibrium, since smaller drops will tend evaporate and larger drops will 
tend grow. The critical radius usually calculated using Kelvin’s equation 
(Thomson 1870) 


20 V, 


where the critical radius, the gas constant per mole, the specific 
surface energy the droplet (Frenkel (1946, 368) has stated that this 
equivalent the ordinary surface tension corresponding the 
liquid molar volume, the vapor pressure over drop radius and 
the saturated vapor pressure over plane surface the absolute 
temperature Equation (3.4) has been derived also Frenkel (1946, 368) 
and Farley (1952). 

its derivation, the density the vapor neglected comparison the 
liquid density, and the vapor treated ideal gas. hardly likely that 
the assumption ideality justified here, since the condensation process 
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clearly requires interactions between molecules. Moreover, should noted 
that for substances which condense small supersaturations, eq. (3.4) 
suggests very large critical radii, tending infinity the supersaturation 
tends unity. order make some comparison this theory with our 
experiments, show the result Table the critical radii, calculated 
according eq. (3.4), but replacing our experimental values the 
supersaturation 

Frenkel (1946, 373) has also given expression for the critical droplet 
radius, which may written 


2o0ViT 
(3.5) = LAT 


where the latent heat per mole, and the undercooling required 
produce condensation. Equation (3.5) derived from eq. (3.4) with the aid 
the Clapeyron—Clausius equation (which does not hold exactly here because 
the Gibb’s functions are not exactly equal inside and outside the drop 
(Pippard 1957)). Table shows the critical radii, calculated using eq. (3.5) 
and our directly measured critical supercoolings, will noted that 
these values are between 1.5 and 2.1 times smaller than those given eq. (3.4). 
The latent heats, used were those Berman and Mate (1958) for the three 
highest temperatures, and the smoothed values collected van Dijk and 
Durieux (1957) for the lower temperatures. 

The number molecules, droplet critical radius also shown, 
Table for and for These numbers are very high compared with other 
substances, e.g., for water, where less than 100. (Das Gupta and Ghosh 
(1946, 238) calculated for water, but their calculation involved 
numerical slip.) 

Farley (1952, 538) has calculated the number drops per cubic centi- 
meter, supersaturated system which may expected reach the 
critical size gives 


(3.6) 


where the number vapor molecules per and Boltzmann’s 
constant. From the data Table equal drops per 
4.181° and drops per 3.147° For lower tempera- 
tures the values are even smaller. This implies that helium droplets 
large could appear statistically. 

therefore conclude that the present state the theory homogeneous 
nucleation inconsistent with our data. Whether this defect the under- 
lying assumptions, or, seems more likely, due the gas’’ and 
other approximations made the theoretical treatment, not present clear. 
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CONDENSER MEMORY UNIT FOR IMPROVING 
SIGNAL-TO-NOISE 


ABSTRACT 


The unit contains low-leakage condensers which can store signal for 
several hours necessary. the signal repeated over and over, the successive 
scans can added and the signal-to-noise ratio builds the square root 
the number repetitions. principle, the final signal-to-noise ratio only 
slightly better than would obtained from single scan stretched out 
the same total time, but practice the result may considerably better, 
especially the signal fluctuates slowly. has been used successfully several 
investigations twilight spectra with photoelectric 
spectrometers. The original version took minute per scan and was rather 
bulky; recent modification can scan channels seconds required. 


INTRODUCTION 


problem which arises many lines research the measurement 
repetitive signals the presence noise. have mind especially the 
photoelectric spectra aurora and twilight, similar faint and somewhat 
variable sources, and will discuss the matter with reference this particular 
case; thus, one repetition will called scan. the signal-to-noise (S/N) 
ratio inadequate can improved increasing the time observation, 
and this normally done slowing down the rate scan and increasing 
the time constant the filter. But the source aurora, the total intensity 
may vary much during the scan that the results are useless, even the 
length only minute; and this time often insufficient. Although twilight 
does not vary nearly quickly, great difficulty would still present 
interpreting 10-minute scan. The other possibility take large number 
short scans and average them, what essentially the same thing, add 
them together. This has sometimes been done manually from number 
recorder tracings, but automatic method obviously much pre- 
ferred. not difficult show that for constant source this scheme gives 
the same S/N ratio single scan taking the same time, but for fluctuating 
source the averaging method far superior. addition, the method 
described lends itself naturally time integration the signal each channel, 
and thus provides further small factor improvement over the conventional 
one. Any device which can add successive scans must able store informa- 
tion and recall when needed; convenient call and our 
device will therefore called memory”’, since the information 
stored the form voltages bank low-leakage condensers. The 
stored spectrum divided into number channels, one for each condenser 
the bank scanned synchronism with the scan the spectrometer. 

received November 11, 1959. 
Contribution from the Department Physics, University Saskatchewan, Saskatoon, 
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suitable device adds the new spectrum the old sum and stores the new 
sum back the memory. This adder the principal novelty the instru- 
ment described. The natural field application this memory 
scans several seconds longer, and fairly small number channels 

interest consider some previous methods which have been used 
for similar purpose; the survey does not pretend complete, however. 
photographic time exposure the repeating trace oscilloscope de- 
scribed Bloch and Garber (1949); they were observing nuclear induction 
signals, but the same technique fairly common with radar displays. While 
simple, the method inefficient because the individual signals are always 
added even when they should, part the time, subtracted. Suryan (1950) 
used magnetic recording drum rotating synchronism with the signal; 
the recording head was moved slowly parallel the axis that the signals 
were impressed along helical track. The reading head covered the whole 
drum and added all the signals together. The biggest difficulty with this 
scheme its mechanical complexity; but not suitable any case for 
signals which repeat slowly ours, since the drum would rotate slowly 
that pickup from would too small. Harrington and Rogers (1950) stored 
radar signals storage cathode-ray tube; sense, this memory the 
condenser type, but the condensers are small that the signals leak off 
few seconds the most; also, the fidelity not high. dielectric recording 
drum described Anderson (1957) much better both respects but was 
not available when our apparatus was being developed. 

direct ancestor our system the memory commutator Beard and 
Skomal (1953). Though rotated much faster, the important difference 
that they fed the condensers through relatively small resistor and allowed 
each condenser charge exponentially. Thus the input signal constant, 
the voltage the condenser will approach asymptotically. our system, 
the voltage continues rise until the circuit overloads, and can many 
times larger than the input; the stored voltage the true time integral 
the input. The device Beard and Skomal gives the same approximation 
this integral that circuit does. The Kozak 
(1958) rather similar, though used for entirely different purpose. 

The principle our method the following. each condenser switched 
into the circuit, the voltage duplicated second con- 
denser special amplifier unit gain. Then the transfer condenser 
switched series with the incoming signal and the sum voltage duplicated 
the memory condenser the same amplifier. The action accurate 
enough that stored signal can read off and back again large 
number times with little degradation. The condensers used 
styrene dielectric that the effect small and signal can 
stored for hours without leaking away. The unit has been found valuable 
making accurate measurements twilight spectra, adding from 
one-minute scans for total much hour. Application 
aurora has been prevented the scanning mechanism which could not 
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made around seconds; scans are too long useful with 
aurora. The new switching scheme described below removes this limitation. 


OPERATION AND CIRCUITS 


The most natural way add series voltages condenser feed 
them succession through constant-current device, such large resistor 
suitable type feedback amplifier. But for bank condensers this 
has the serious disadvantage that the voltage generated inversely propor- 
tional the capacitance, that all the condensers must closely matched. 
The method described operates voltages instead currents and 
therefore independent capacitance. Moreover, gives readout the 
stored signal which may observed during operation. based pre- 
cision voltage follower shown the triangle Fig. This device accepts 
signal very high-impedance level, drawing essentially current, and 
reproduces low impedance with negligible error, estimated less 
than 0.01%. 


readout 


Fic. (a) Effective circuit the memory while reading from one condenser the bank. 
The triangle represents voltage follower with gain very close 
Circuit for adding the input and storing the result back the bank. 
(c) Actual circuit shown with the relay the position. 


First let the follower connected Fig. represents one the 
condensers the memory bank and condenser. Whatever 
voltage will reproduced and also available for monitoring. 
the circuit then switched into the arrangement Fig. the voltage 
connected series with the input from the spectrometer, and the 
sum will appear Cy. The next time this action repeated, new 
the circuit. The switching done the three-pole relay shown Fig. 
which drawn the position. The input voltage does not come directly 
from the spectrometer, but passes first into feedback integrator (Korn and 


(a) 
(b) 
” 
(c) 
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Korn 1952) which reset zero just after the end part the memory 
cycle. Thus, the voltage each condenser the bank the integral the 
signal for the total time has been the circuit. 

The follower itself essentially high-gain direct-coupled amplifier with 
100% negative feedback. The circuit shown Fig. was developed after 
brief qualitative description Bousquet (1956). can regarded 
three-stage cathode follower; the output fed back the first tube the 


J 
Ground 


Fic. Circuit the follower. 


common cathode connection. For some purposes also useful regard the 
circuit the following way: transfer the ground from the end the source 
the common cathode connection, leaving the source floating between the 
input and the plate the third tube. The essentials this arrangement are 
shown Fig. The 220-kQ resistor now acts the plate load the third 
stage, even though the other side its power supply. The amplifier 


NORMAL GROUND 


The circuit Fig. with different ground point. 


thus ordinary three-stage one except that feedback applied 
through the floating source. The loop gain will about 10° for the circuit 
shown; thus the external gain should 0.99999. This large loop gain means 
that there great danger high-frequency oscillation; the condensers the 
plate circuits control the shape the high-frequency cutoff prevent this. 


160 
100 
ol 
400K 200K 
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The high loop gain realized only the plate and bias power supplies can 
float along with the output indicated Fig. have used single 350- 
volt supply feeding two miniature regulator tubes; the complete circuit thus 
includes five tubes and rectifier. Normal grounded supplies provide the plate 
and cathode-return voltages the last stage; they need not regulated. The 
rheostat the cathode the 6C4 allows the level the output set 
equal that the input; this must done rather accurately the voltages 
stored the memory will drift twice the offset each time the information 
cycled. Often small drift this sort not important, since all channels 
are affected the same, but large drift might take the circuit out its range 
correct operation. The first stage operates with low heater and plate voltages 
that the grid current less than ampere. 

the output the follower connected condenser considerably 
different potential, may not able supply the required current and the 
output will not able follow the input for moment. During this moment 
the input may draw large burst current. eliminate this trouble, the 
series resistors shown Fig. 1(c) are included; they must small enough 
let the condensers charge fully during the time available. More charging 
current could made available replacing the third tube with power 
tetrode; this would useful memory with very short time per channel. 

SWITCHING 

The very low leakage polystyrene condenser realized only there 
are leakage paths parallel with it. This dictated the design the com- 
mutator the first version, since ordinary switches adequate performance 
were available. The condensers were arranged circle and contact springs 
attached and supported the inside terminal each; the arrangement 
was similar the sketch Fig. arm carrying another contact was 
rotated steps make connection each turn. The intermittent motion 
was obtained from synchronous motor driving through six-point Geneva 
movement and then 1:5 gear ratio. This operated well rate one step 
seconds, but would not quite keep time three steps per second. Probably 
slightly more powerful motor would have been able run the fast speed, 
but this was not pursued. 

The shaft which drove the Geneva movement also carried cam which 
operated microswitch and thereby the relay shown Fig. 1(c). was left 
the position while the commutator was switching, and for short 
time after, that the transfer condenser could become fully charged; the 
part the cycle occupied the rest the time, little less than half. 
Just after the end the operation, the integrator the input was 
reset zero another relay making brief contact. first this was operated 
second cam, but later transient generated when the transfer relay 
moved the position. Not only this simpler; also guarantees 
synchronization between the two relays. 

much more compact and faster arrangement possible taking advan- 
tage the properties the dry-reed switch, suggested 
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Helmer and Hemmendinger (1957). This switch has the high insulation 
resistance needed, and the commutator may ordinary one occupying 
little space. The high speed the Glaswitches suggests their use also for 
resetting the integrator. The commutator available required the use 
channels instead 30. 

worth while pay some attention the method setting all the 
condensers zero initially. The easiest build switch ground the 
wiper the commutator, the equivalent point the Glaswitches are 
used. But often inconvenient have wait for complete cycle the 
memory, and separate shorting bar which can connect all the condensers 
once preferable. 

switch should provided lock the transfer relay the position. 
Then the contents the memory remain undisturbed and may desired 
observed one more times. This switch can also used start the memory 
operating any suitable time after has been set zero. 


OPERATION 


The contents the memory are easily recorded 
ometer; the output the follower low-impedance point and can drive 
calibrated voltage divider helical potentiometer give suitably 
low voltage. The recording memory should made high voltage 
level possible reduce the importance drift and other stray voltages; 
the output may then attenuated the magnitude required the recorder. 
often useful follow the accumulation the signal the recorder; 
any large interfering pulses are stored, they can allowed for. Records 
this sort show two steps per channel, corresponding the old and new con- 
tents. record made the end the run with the transfer relay off shows 
only one. The final recording may repeated several times desired until 
the best settings gain and time constants are found. 

example run the sodium lines the twilight shown 
Fig. (Lytle 1958). The way which the S/N ratio builds more and 
more spectra are included may observed. The lines appear above con- 
tinuum caused nonresonant scattering the atmosphere. The width 
the lines and their triangular shape arise entirely from the spectrometer; 
there detectable degradation introduced the memory except for the 
steps corresponding the individual channels. The area line recorded 
this way good measure its intensity; easily found adding 
all the steps within the line. The fact that the intensity ratio the two 
lines much less than the normal 2:1 results from self-absorption the 
sodium layer the upper atmosphere (Chamberlain, Hunten, and Mack 
1958). The measurements could not have been made with nearly much 
accuracy without the memory; single 12-minute scan would have been 
hopelessly distorted variations intensity and atmospheric transparency. 

The memory was also used Harrison study the band discovered 
him the twilight 1.58 microns (Vallance Jones and Harrison 1958). 
was able use much better resolution than would otherwise have been 
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Fic. Part run the lines twilight. 


possible. The intensity decays with half-life about minutes, that 
single long scans are again not practical. But scans minute each added 
the memory gave useful result and allowed the band identified. 

search for the potassium resonance line 7699 twilight has recently 
been carried out Lytle (1958). this case the radiation was not detected, 
but use the memory enabled the limit detection set much lower 
than would otherwise have been possible. appears that other studies 
twilight may also made possible this new technique, which offers 
considerable increase available sensitivity and permits photoelectric spectro- 
meters compete nearly equal terms with fast spectrographs. expected 
just valuable for studies the auroral spectrum now that can 
used with 10-second scan. Until now, only the brighter auroral forms gave 
large enough signal for quantitative studies; the fainter forms are not only 
much more plentiful, but may also give different results. 

Many studies would profit from the use memory with more channels 
than the used here; and faster repetition rate might also useful. 
doubtful the present adding unit would operate well beyond three 
four channels per second; and might economic make condenser 
memory with channels even 96. likely that the capacitance the 
condensers could reduced 0.1 microfarad without impairing the operation; 
but there point reducing their cost much below the cost the switching 
system. With this mind, the best size probably about 0.2 microfarad. 

great increase both speed and number channels possible 
adopting digital-computer techniques, has been done widely recent years 
for multichannel pulse-height analyzers (Hutchinson and Scarrott 1951; 
Schultz, Pieper, and Rosler 1956). The memory requirements for the two 
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purposes are almost identical, except that fast random access any channel 
important most cases nuclear physics, but importance what- 
ever here. The ultrasonic delay-line memory then just good the more 
expensive magnetic-core type. has also the advantages that the stored 
spectrum continuously displayed cathode-ray tube, and that informa- 
tion may read and out arbitrary rates, simultaneously desired. 

speed five channels per second sufficient, another possibility 
record the individual spectra digital form punched paper tape, similar 
teletype tape. This may later read into general-purpose digital com- 
puter and the averaging done suitable program. This requires much 
less equipment, long the computer available, but has the disadvan- 
tages that the averaged spectrum not immediately visible; anything has 
gone wrong, will not known until the next day. the observations are 
experiment that can repeated, this little importance; but 
may very important for observations uncontrolled natural phenomena. 


pleasure acknowledge the help several colleagues with the 
construction and testing the instrument: Harrison, Lyon, 
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IMPURITY CONCENTRATION PROFILES ICE 
ANTHRONE METHOD! 


ABSTRACT 


Solutions sodium carboxymethyl cellulose (CMC), 
ethylene glycol, and ethyl alcohol water were prepared and frozen under con- 
trolled conditions. The ice each case was divided into horizontal sections, 
melted, and the amount each impurity measured. The organic liquid and salt 
concentrations were found from refractive index values using dipping refracto- 
meter, and the CMC content, after treatment with anthrone solution, was 
measured with spectrophotometer. Developments the anthrone method 
are described. Qualitative tests were also made freezing solutions soluble 
dyes. 

the freezing aqueous solutions CMC and the dyes, rejection the 
impurity almost complete. When solutions salt alcohol glycol water, 
CMC and alcohol glycol water, are frozen surface concentration 
peak and second peak about below the surface, corresponding 
previously observed discontinuity the crystal structure impure ice, are 
found. The ability alcohol glycol carry CMC dye places the 
ice structure not available them the absence the organic liquid noted. 


INTRODUCTION 


has been shown Perey and Pounder (1958) that when water con- 
taining added foreign material frozen, transition layer indicating radical 
change crystal structure clearly visible distance the order 
below the surface. was assumed that such change associated with 
changes concentration the added foreign material. The present work 
was undertaken verify the assumption quantitative measurement 
the distribution the foreign material the ice. 

seems probable that impurity distribution major factor deter- 
mining the mechanical properties any form impure ice (see Pounder 
1958). This particularly true the extremely important case sea ice. 

According theory suggested Anderson (1958) and elaborated 
Assur (1958) the properties sea ice are very largely dependent the 
distribution salts and brine within the ice. the present work additive 
materials other than those found sea water have been studied but one 
example included where pure sodium chloride was the additive and 
hoped extend the work this direction soon. 


METHODS MEASUREMENT 
The impure melts studied Perey and Pounder (1958) contained various 
alcohols and other organic solvents, and variety carbohydrates, par- 
ticularly sodium carboxymethyl cellulose (grade D435 supplied Chemical 
received June 29, 1959. 
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Developments Canada, said contain minimum 66% sodium carboxy- 
cellulose, the other salts present being principally sodium chloride and 
sodium glycolate). 

was decided confine the present quantitative analyses largely 
D435, ethyl alcohol (Reliance Chemicals 650P 823% volume), and 
ethylene glycol (Purified Fisher Laboratory Chemical). Since the behavior 
freezing CMC and sodium chloride quite different far rejection 
concerned, unfortunate that grade D435 contains some sodium chloride. 
grade CMC free sodium chloride was available, however, and any 
event the concentrations D435 used freezing tests were low (about 
weight) that the small amount sodium chloride contained 
introduced negligible errors the tests for CMC and alcohol glycol. All 
the tests were made with D435 samples from single batch that the 
composition the D435 was assumed uniform. 

Ethylene glycol was chosen provide, with the alcohol, examples both 
relatively non-volatile and volatile solvents; they serve also examples 
both heavier than water and lighter than water solvents. 

The method which was found the most suitable for D435 estimation 
was spectrophotometric technique which will referred the anthrone 
method. This method consists essentially adding anthrone solution 
the test samples. When carbohydrate (CMC these tests) present, 
green color develops through chemical reaction between the anthrone and the 
carbohydrate. The optical density the solution, measured with the spec- 
trophotometer, can used calculate the concentration CMC and hence 
D435 the test sample. was found convenient use the concentration 
D435 calibrating the method well citing results. 

The anthrone method was suggested Dreywood (1946) the Eastman 
Kodak Laboratories qualitative general test for carbohydrates. has been 
carefully developed for quantitative use Black (1951) the DuPont 
Laboratories and frequently quoted research literature. 

Its use, however, presents several difficulties; the results vary for example 
with the temperature which the reaction takes place, with the age the 
reagent and the solution, with the acid concentration, and with the degree 
substitution the sodium cellulose. The presence various 
other carbohydrates their derivatives interferes with the reaction. Black 
states that calibration curves should made every weeks and that the 
reagent used should between and hours old. Under these conditions 
Black achieved relative accuracy 2%, the degree substitution was 
known. The anthrone dissolved sulphuric acid and the quantities used 
are large (50 60% acid for each test), that day testing might 
necessary use, and dispose of, about 1500 60% sulphuric acid. 

modification the method has been suggested Loewus (1952). Loewus 
used acetate solvent for the anthrone and states that this prevents the 
gradual darkening the reagent which occurs when the solvent sulphuric 
acid, darkening which varies with any slight changes the purity the 
sulphuric acid. The expected accuracy not stated his paper. 
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used tests this laboratory Loewus’ modification was very disappoint- 
ing accuracy, consistency, and with regard the aging the reagent. 
Further the results were affected the alcoholic content the samples 
tested, making its use impossible for this investigation. 

two recent papers, Scott and Melvin (1953) and Mokrasch (1954) re- 
turned dissolving anthrone sulphuric acid but stored cool dark 
place and emphasized the importance controlling the temperature 
samples and reagent before mixing. Scott and Melvin used 90°C for the 
development color; Mokrasch used 80° 

Scott and Melvin gave detailed account sources errors and felt that 
standardization necessary for each set samples. they could not com- 
pare one day’s absorbance with another indication precision, they 
measured the standard deviation sets duplicate blanks and 345 sets 
duplicate sample determinations 0.48% absorbance 0.600 absorbance. 
They found their most important source error was variation readings 
obtained blank tests. Mokrasch found his reagent stable for months 
when stored the dark. gave estimate accuracy. 

Roe (1955) has suggested adding thiourea the reagent means in- 
creasing its stability and found that this way could used for weeks 
when stored refrigerator. states that when the heat mixing used 
means developing color, the accuracy expected only 
Roe’s (1954) careful experiments with dextran the coefficient variation, 
using duplicates, was +0.36%. 

The following procedure has been found most satisfactory this laboratory 
and was used the experiments reported. 


Preparation the Reagent 

volume concentrated sulphuric acid (Fisher Reagent 
diluted with distilled water. Both volumes are measured into 
beakers room temperature, and the beakers are then cooled crushed ice 
before mixing and again after mixing. With the beaker diluted acid re- 
maining the ice, 0.200 anthrone (Fisher Reagent grade) and 2.000 
thiourea (Fisher Certified Reagent) are added and stirred until the solution 
clear and dissolving seems complete, which takes about half hour. The 
used within week preparation. 


Test Procedure 

Samples each the solutions tested for carbohydrate content 
are pipetted into screw-top pyrex test tubes. quantity the re- 
agent measured into each several additional test tubes, using burette. 
All test tubes are then cooled for minutes crushed ice, after which the 
test solution and reagent samples are mixed pairs (by pouring backwards 
and forwards twice) and returned the crushed ice. The test tubes, with 
their tops loosely, are arranged rack and plunged into boiling water for 
minutes, returned the ice bath again for minutes, and then allowed 
stand (covered exclude light) room temperature for minutes. After 


this time, the absorbance each sample measured 620 using Beck- 
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man Model spectrophotometer. This procedure, developed after many 
trials, insures precisely controlled and reproducible temperature regime 
for the color-producing reaction, and gives quite reproducible results for the 
carbohydrate content. 

each experiment one blank included using distilled water instead 
test solution. The optical density the blank, compared with the reading 
given distilled water (without any reagent), serves check the 
experimental conditions and the condition the reagent. great im- 
provement consistency obtained using the optical density the 
blank zero, rather than that pure distilled water. 

The use thiourea one the components the reagent not only in- 
creases the length time for which the reagent can kept without deteriora- 
tion, but also decreases the variation between successive tests the same 
solution. 

The measured optical density for D435 solutions was not affected the 
presence ethyl alcohol the concentrations tested, alcohol. 

The method very sensitive. concentration about 200 
possible distinguish differences optical density corresponding 
differences concentration one part million, but the 
reproducibility not high. There are the uncertainties introduced the 
measurement the volumes test samples and reagent and slight 
variation the heat treatment, etc. Several authors mention the erratic 
behavior reagent aging, and color stability, and unaccountable varia- 
tion blanks. For single tests the possible error estimated 

The calibration curve given Fig. All D435 solutions were filtered 
through No. filter paper, this was found necessary obtain 
reproducible results, probably because D435 does not form completely 
clear aqueous solution. However, the concentration figures plotted the 
graph are the original values before filtration. The curve was checked over 
period months, using number preparations both reagent and 
D435 solutions. Each value plotted thus the average 
measurements. The deviation individual readings from average values 

Tests were made with four different commercial grades sodium carboxy- 
methyl cellulose (D235, D435, L245, T253 made Chemical Developments 
Canada), various concentrations and using solutions both filtered and 
unfiltered. 

Using filtered solutions, and plotting optical density against the actual 
amount sodium carboxymethyl cellulose each preparation, according 
the specifications, was found that the values for all four 
preparations lay along one straight line. This shows that the anthrone test 
for sodium cellulose not affected the presence varying 
amounts other salts such are found these preparations. 

Figure illustrates the effect filtering, using D435. Points marked 
crosses were using unfiltered solutions; points marked circles 
represent filtered solutions. The concentration values plotted are the values 


before filtering, each case. 


358 


DENSITY 


OPTICAL 


44 
| 

vA 

| 50 100 150 200 25 

Fic. 
5+ UNFILTERED SOLUTIONS 
O-FILTERED SOLUTIONS 

> 

z 

a 

° 


+ 


CANADIAN JOURNAL PHYSICS. VOL. 38, 1960 


50 100 1s0 200 250 
CONCENTRATION OF 0435 IN /mi 


Fic. 


SMITH AND POUNDER: IMPURITY CONCENTRATIONS ICE 359 


the beginning this investigation alcohol concentrations were measured 
means distillation and specific gravity measurements, method both 
lengthy and uncertain. For most the work here reported Bausch and 
Lomb dipping refractometer was used (with prism range 
refractive index, and occasionally No. auxiliary prism). it, rays light 
(from incandescent lamp) pass through beaker containing the sample, 
and enter the prism. The position the grazing critical ray forming the 
boundary between the light and dark parts the field view may 
measured. 

calibration curve for alcohol concentration was prepared using solutions 
containing both D435 and ethyl alcohol the proportions used ice prepara- 
tion, since the D435, probably the small amount sodium chloride it, 
does raise the values trifle over those found when using the alcohol alone. 
This effect very slight, being undetectable for D435 concentrations less 
than per ml. 

Using water bath and measuring temperatures 1/10 Centigrade degree 
when equilibrium was established, measurements could made 1/100 
scale division. using average determinations these were usually 
reproducible +0.01 scale division. When the auxiliary prism was used 
(in cases when only small amounts liquid were available) the uncertainty 
became least +0.02 scale division. Similar calibration curves were 
prepared for aqueous solutions ethylene glycol and sodium chloride, 
was found that the refractive index these solutions changed appreciably 
with concentration. This method much more sensitive, more reliable, and 
far quicker than the specific gravity method which had previously been used 
for sodium chloride determinations. variation +0.02 refractometer scale 
divisions corresponds, 25° +0.016% ethyl alcohol, +0.008% ethylene 
glycol, sodium chloride. these three cases the refractive index 
varies linearly with concentration the ranges tested. All concentrations 
were converted standard temperature 25° 


PREPARATION ICE 


For each ice study 1.000 D435 was dissolved 975 distilled water. 
This was filtered through Whatman’s No. filter paper and the volume made 
975 ml. this was added alcohol and 5100 distilled water, 
ethylene glycol and 5100 water, 5125 distilled water. 

About 100 was retained for testing, and the other liters solution 
were placed rectangular plastic vessel about deep, surrounded 
sufficient insulating material ensure essentially one-dimensional cooling. 

The solution was frozen constant ambient temperature, usually —28° 
until solid. The ice was then taken cold room kept and the 
sides the ice block were trimmed off with band saw, eliminate edge 


effects. 

some cases the ice was then divided into horizontal sections with the 
band saw, the sections were placed screw-top bottles and then melted 
room temperature and analyzed. The minimum thickness sections obtain- 
able with the band saw (the usual section cut this way was mm) 
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and about 1.5 ice lost between sections the saw cuts. the ice 
soft the uncertainties thickness increase. 

Where the concentration impurities the ice changed rapidly with 
depth, thinner sections were prepared either planing the surface the ice 
melting. Planing very tedious, slow process, but very fine sections 
may obtained this way. For melting, the block ice was held vertically 
horizontal pyrex tray which was gently warmed from below, the melt 
being poured off small batches. The thickness the layers melted was 
determined from knowledge the dimensions the block and the weights 
the samples. Great care needed keep the block vertical melting 
proceeds. 

The advantages melting and planing are that they enable very thin layers 
separated, and that material lost between sections with saw 
cuts. might expected that the volatile ethyl alcohol would partially 
lost these methods, but the values obtained were consistent with samples 
prepared sawing. 

RESULTS 

Concentration profiles are given for four blocks ice containing water, 
D435, and ethyl alcohol; for two blocks ice containing distilled water and 
D435 only; and for one block containing pure sodium chloride and water. 
all these cases the blocks were frozen ambient temperature —28° 

addition profiles are shown for three blocks ice containing D435 and 
ethylene glycol, which were frozen ambient temperatures 
—19° and —5° respectively. 

Some samples alcohol and water and glycol and water were frozen 
insulated vessels (cylinders about diameter and 7.5 deep) and the 
additive (alcohol glycol) concentration the ice measured freezing 
proceeded. 

visual indication distribution, water-soluble dyes (methylene blue, 
violet, and eosin) were dissolved pure water and water containing 
trace ethyl alcohol, and the solutions frozen. 

The concentration profiles for ices II, III, and each containing D435 
and ethyl alcohol, are shown Fig. where crosses indicate that the samples 
were prepared planing, circles sawing, and squares melting. The 
results for ice and ice which were made from the same stock solution 
have been combined. This was done because the alcohol values were few and 
unreliable, some them having been measured distillation and specific 
gravity measurements before the dipping refractometer was available. 

will seen that there inconsistency between values for samples 
prepared the different methods. Values for ice and ice fit well into the 
same curve. 

The graphs show that the profiles are similar for both D435 and alcohol. 
approximately the position the visually observed transition layer there 
marked peak the concentration both D435 and alcohol. will also 
seen that there high concentration both added substances the 
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surface layer. There considerable variation the values obtained the 
surface and this appears result from variations the thickness the 
surface layer examined. seems probable that the high concentration con- 
fined very thin layer, perhaps even monomolecular layer, and that 
the concentration drops sharply below this. 

Below the transition layer the concentration decreases rapidly value 
which remains fairly constant throughout most the depth the block. 

ensure that the variation the concentration with depths conse- 
quence only the freezing process, solutions containing D435 and ethyl 
alcohol, and D435 and ethylene glycol, both the same dilution used 
the freezing tests, were left undisturbed room temperature for some time. 
The vessels were covered and insulated guard against evaporation and 
sudden changes temperature. After periods weeks, samples care- 
fully pipetted from various levels showed surface concentration nor any 
variation concentration with depth. 

Two tests and VI) were made with ice containing D435 alone. The 
results were almost identical and are illustrated for ice Fig. slight 
concentration was observed the surface layer but below this the D435 was 
almost totally rejected from the ice until near the bottom the vessel where 
the concentration rose rapidly. Visual inspection foretold this result before 
the block ice was sectioned. both tests, most the ice was clear and 
glassy and there was bottom layer which was white and opaque. 
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Figure shows the concentration profile for ice VII. This example 
the results obtained freezing pure sodium chloride brine concentration 
(parts per thousand). this case the freezing was done much 
deeper vessel that the brine was initially 21.6 deep and the very high 
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concentrations salt found towards the bottom the vessel after freezing 
not appear over the depths here plotted. this case also there was 
transition layer clearly visible above the level and the graph shows 
corresponding concentration the sodium chloride. There is, however, 
evidence surface layer concentration. 

Tests ice frozen from melts containing D435 and ethylene glycol were 
undertaken avoid the uncertainty losses evaporation, which are 
involved the study ice containing alcohol. Three concentration profile 
pairs are shown: for ice VIII (Fig. the ambient freezing temperature was 
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—28° for ice (Fig. was —19° and for ice (Fig. 
with the alcohol samples discussed previously, the D435 concentration follows 
closely that the organic solvent, although the shape these profiles differs 
considerably from those obtained the alcohol tests. With alcohol the con- 
centration throughout the bulk the ice was uniform, whereas with glycol 
rose more less uniformly with depth. This indicates some rejection 
glycol (and D435) the ice, rejection which more pronounced when 
freezing takes place more slowly. each the three tests some glycol was 
expelled completely from the ice, remaining liquid film the bottom 
the vessel. The profiles all show evidence surface concentration and also 
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change slope the curve. ice and ice the transition layer was 
observable visually vertical sections through the blocks and was seen 
quite jagged. the analysis thin horizontal sections covering fairly broad 
area this irregularity would have the effect flattening the transition peak. 

Visual examination ice and ice also showed patches clear crystal 
growth surrounded regions containing numerous air bubbles. appears 
probable that thin sections may give erratic values such non-uniform ice. 
This may account for points lying off smooth curves. Note that Fig. the 
same irregularities occur both glycol and D435 concentrations. 

study possible evaporation losses, concentration measurements were 
made for large number samples alcohol and water and glycol and 
water, frozen small cylindrical vessels diameter 7.5 deep). 
Freezing was done with both covered and uncovered vessels, and number 
the tests the surface layer ice was removed repeatedly for analysis while 
freezing was proceeding. 

These experiments showed that the first ice that forms the surface 
solutions the form fairly compact layer consisting 
horizontal thin ‘‘sheets’’ which are well bonded together when the alcohol 
content the solution low (up about 0.7%), but are weak and easily 
separated when the alcohol content the solution high (about more). 

The alcohol content this surface ice layer, and successive ice layers that 
form after each turn removed, lower than the average solution con- 
centration and surprisingly constant (about 0.4% with initial concen- 
tration 0.7%). 

For closed vessels, analysis the ice and the unfrozen melt confirmed 
that alcohol had been lost evaporation. For open vessels the alcohol 
concentration the surface ice does decrease slowly with time. (For example, 
two similar samples ice were stored, with and without covers, cold 
room for week. The concentration alcohol the top 1-cm layer 
was then 3.23% and 2.90% respectively. Below about 1.5 the concentra- 
tion remained alike the two vessels.) 

After the ice layer reaches thickness the order cm, the ice begins 
grow long pendent vertical plates. these experiments, the plates were 
1.5 width, very thin (probably about 0.1 thick), and crystal clear. 
They were observed vertical length and were spaced about 
0.5 apart, groups least roughly parallel. The alcohol concentration 
these plates, strained from the solution (so that inevitably little liquid 
adhered), did not differ from the value found for the surface ice layers. 

When the solution allowed freeze completely, alcohol extruded from 
the block and remains liquid film the inside the vessel. 

The concentrations alcohol found throughout the ice, when the sample 
was allowed freeze completely, were somewhat higher than the concen- 
trations found ice formed similar solutions, when only part the melt 
was allowed freeze. This probably owing the fact that horizontal 
layers cut from solid block ice contain both ice crystals and some con- 
centrated alcohol solution trapped between them. The sections analyzed were 
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relatively thick, that the fine detail Fig. was not obtained but the 
concentrations were consistent with the profiles Fig. 

Melts containing glycol and water were frozen and analyzed similar 
fashion. The freezing pattern found was very much like that for 
water mixtures; the horizontal layer and its pendent vertical plates contained 
equal concentrations glycol, and this concentration was surprisingly similar, 
quantitatively, that the alcohol case. When glycol melts are frozen 
solid, however, the concentration profile below the 2-cm level markedly 
different; comparing ice (Fig. with ice VIII (Fig. evident that 
below the transition layer the alcohol values remain constant throughout 
most the ice, whereas the glycol values increase steadily with depth. 

the ice crystals form, both alcohol and glycol are rejected. But the 
alcohol, free move between the pendent ice crystals, will tend rise, since 
lighter than water, whereas the glycol being heavier than water will 
tend sink. the case alcohol then, the rejection and gravity effects are 
opposition, and throughout the bulk the volume, measurements the 
alcohol content horizontal layer give almost constant value. 

the case glycol, where the effects rejection and gravity move the 
glycol the same direction, the trapped solution the lower layers will 
have higher concentration than those above. 

Solutions containing water-soluble dyes (methylene blue, methyl violet, 
and eosin) were next investigated, with most striking results. When the solu- 
tions contained only dye and water the dye was strongly rejected 
freezing proceeded that the ice was quite clear and colorless except the 
bottom where the dye was concentrated sharply defined region. But when 
these solutions contained both dye and ethyl alcohol, the ice appeared 
uniformly colored throughout. 

Examination this ice was attempted under the microscope the cold 
room. When sufficiently thin sections could set the right angle, 
was possible see that the ice consisted clear platelets colorless ice 
separated thin highly colored films. 


DISCUSSION 


The concentration profiles shown this paper give strong evidence for 
peak concentration values the position below the surface found Perey 
and Pounder (1958) the position transition layer the crystal 
structure. this transition layer rapid change occurs both crystal sizes 
and orientations. most the ice samples tested concentration peak was 
also found the surface layer. The absolute size this concentration un- 
certain, because the value found each case must function the thick- 
ness the layer ice which was measured. The thinner the layer the 
higher was the measured concentration, suggesting that the concentration 
maximum very thin, possibly even monomolecular, surface film. The 
same argument applies the value found for the transition layer peak, with 
added complication. Visual examination vertical section usually showed 
the transition layer very thin but not horizontal plane. Sometimes 
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was quite jagged, occasionally varying depth much half centi- 
meter. the transition layer itself could isolated, one would expect find 
very much higher concentrations than any here recorded. 

has been shown that when pure aqueous solutions sodium carboxy- 
methyl cellulose (CMC), various water-soluble dyes, sodium chloride, ethylene 
glycol, and ethyl alcohol are frozen, the amount solute trapped the 
ice varies widely. For the concentrations and freezing rates used, the CMC and 
dyes are virtually completely rejected the ice. Some salt trapped; 
Fig. the region from about depth might considered normal 
and shows that about 30% the salt was retained the ice. The concentra- 
tions alcohol and glycol vary too much with depth speak normal 
concentration but reasonable guess that 90% the alcohol and possibly 
60% the glycol are trapped the ice. These results suggest the 
possibility that the solubility water the additive factor— the CMC 
and the dyes are slightly soluble water compared salt and the two 
organic liquids are infinitely soluble water. 

The freezing solutions containing both alcohol (or glycol) and CMC (or 
dye) interesting, since the presence the CMC does not appear alter 
the distribution the alcohol glycol the ice but the distribution the 
CMC now conforms closely that the organic liquid. Apparently the 
CMC carried the alcohol glycol positions the ice not available 
pure CMC. some physical chemical process the small alcohol molecule 
(or possibly several them) takes the huge CMC molecule its 
and persuades the growing ice lattice treat them approximately equals. 
perhaps important note that the number alcohol molecules greater 
several orders magnitude than that CMC molecules, that even 
some the alcohol enters into chemical reaction there plenty left 
alter the refractive index the solution. possible chemical reaction 
esterification the CMC. Although esterification ordinarily reduces the 
solubility water, this not necessarily the case. Noller (1951, 303) points 
out that cellulose insoluble water because strong van der Waals forces 
and proton bonding between chains. Etherification CMC separates the 
chains because the size the carboxymethyl group and reduces these 
forces greatly. Possibly esterification the CMC carries this process step 
farther. possible physical mechanism for retaining the CMC the ice 
through lattice defects. There evidence that alcohol can fit permanently 
into ice lattice, but alcohol affects the properties ice markedly that 
tempting suppose that alcohol can attach temporarily growing ice 
crystal, giving rise imperfections the lattice. these imperfections are 
large the CMC molecules may trapped them, although not course 
forming part the lattice. The present experiments not suggest which 
these speculations the more likely, or, indeed, either them the correct 
explanation. 

The results illustrated Figs. and show that the amount solute 
trapped the ice increases with the decrease ambient temperature and 
hence increase freezing rate. This result well known field studies 


sea ice. 
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NEGATIVE RESISTANCE SEMICONDUCTOR 


RONALD BURGESS 


ABSTRACT 


Negative resistance can appear the static and high-frequency characteristics 
devices which the current determined both voltage and temperature. 
The properties the a-c. impedance arising from the interaction thermal and 
non-thermal effects are discussed and criteria for the appearance negative 
resistance over certain ranges frequency are derived. The application this 
analysis devices which and which not exhibit isothermal negative 
resistance considered. Appendix, relation depend- 
ing two activation energies shown provide quantitative model for 
interpreting the observed turnover behavior germanium diodes. 


The appearance ‘‘thermal i.e. voltage maximum followed 
negative resistance region the characteristics semiconductor devices, 
has been discussed previously (Burgess 1955a, c). the present paper the 
combination thermal and voltage-controlled processes considered generally 
and the results are applied discussion the attainment negative 
resistance. exclude the effects obtained transit time phenomena and 
also the spurious negative resistance which seemingly occurs junction diodes 
under avalanche breakdown conditions and fact arises from the voltage- 
controlled statistics the breakdown process (Burgess 1959). 

The voltage-controlled effects most commonly involved are: impact ioniza- 
tion leading avalanche multiplication, carrier injection and extraction, 
space-charge-limited flow, internal field emission (or Zener effect), and non- 
linear carrier mobility strong electric fields. Except the simple case 
ohmic device low electric fields, one more these effects will present. 

general the current may written and the differential 
resistance which the small-amplitude a-c. resistance zero frequency 


given 


thermal conductance between device and ambient 

The dimensionless quantities and are respectively the representative 
parameters the isothermal characteristic the operating temperature and 
the temperature the current the operating voltage. 

The four alternating sectors positive and negative differential resistance 
corresponds devices which are linear the isothermal condition (e.g. 
thermistors). Usually semiconductor devices will have positive values 
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although some exceptional cases such those involving high electric fields 

positive then the current through the device increases attainment 
negative resistance the current-controlled type obtained. When and 
are the same sign the differential resistance must lie between and 
and this indicates the limits magnitude the attainable negative 
resistance any operating point beyond turnover. 

may also inferred from equation (1) that diode whose isothermal 
reverse does not saturate the differential resistance smaller 
magnitude (whether positive negative) than otherwise identical 
diode exhibiting saturation. Thus avalanche multiplication 
breakdown process compared with that due thermal effects alone but con- 
versely reduces the attainable negative resistance beyond turnover. 

diode displaying isothermal negative resistance the 
forward characteristic narrow germanium junctions (Esaki 1958), the 
differential resistance modified thermal effects once longer small 
compared with both and thermal effects are sufficient for exceed 
and then they will dominate over the Zener effect and produce 
thermal turnover the voltage maximum type. Circumstances can arise 
such device which the characteristic proceeding from the origin could 
exhibit current maximum followed current minimum followed turn 
voltage maximum. 

The above remarks about the characteristics apply frequencies sufficiently 
low for the temperature always steady state with the applied power. 
higher frequencies thermal inertia causes lag the temperature variation 
and the differential resistance derived from the static characteristic must 
replaced the complex impedance. the thermal time constant this 
impedance derived replacing the thermal conductance the complex 
thermal admittance whence 


(2) Z(w) 


This function has semicircular locus the plane with the ends the 
immediately seen that the criterion for thermal effects absent 
not that shall large compared with unity but large compared with 
and y/x. 
The condition for negative resistance appear 


(3) 


this inequality may derived 

(i) the upper frequency limit which device without isothermal negative 
resistance may display negative resistance due thermal effects, 

(ii) the frequency below which isothermal negative resistance 


suppressed thermal effects. 
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linear device the upper limit for purely thermally derived 
negative resistance simply 


. 


(4) 


The presence thermal effects introduces reactive component into the 
impedance which inductive positive and capacitive y(1 +.) 
negative. This may useful embarrassing depending upon the manner 
which the device being used. example benefit occurs very low- 
frequency oscillators using thermistors where the inductive reactance produced 
the thermal inertia permits frequency control simply shunting the 
thermistor with pure capacitance (Burgess 1955c). Similar inductive behavior 
and oscillations point-contact junction germanium diodes 
but here greater than unity due avalanche multiplication reducing 
the upper frequency limit for negative resistance appear (equation (3)). 

any given circumstances the effective time constant differs from 
extent which depends and and also the external circuit that the 
dynamic behavior can only evaluated considering the complete system. 

One type dynamic behavior practical importance the switching time 
bistable negative resistance device between two stable points. The transi- 
tion involved essentially large amplitude embracing the whole the 
negative resistance region and thus outside the scope the present small- 
amplitude theory. However, there still validity the observation that the 
switching time not simply related even when determined thermal 
inertia, since and and the external circuit determine the fashion which 
the variations current, voltage, and temperature are interrelated. 

Thermally produced negative resistance commonly associated with 
relatively long relaxation times the order milliseconds. If, however, the 
current-controlling region small dimensions and material has 
thermal diffusivity (especially true low temperatures) the thermal time 
constant can easily less than microsecond. Conditions this type obtain 
for instance the contact region point-contact diode. 

Furthermore must stressed that the upper frequency limit for negative 
resistance occur may appreciably greater than the plausible value 
conditions where the temperature control the current sharp, 
found from equation (3) that this type situation seen 
reference equation (2) that the heat capacity and not the thermal 


conductance 

possible, though not yet conclusive, that the negative resistance 
and Rediker 1959) and oscillations (Koenig and Brown 1959) 
observed germanium filaments low temperatures have thermal origin. 
The latter authors have suggested that the negative temperature coefficient 
the breakdown field associated with impact ionization leads local 
thermal instability some small region the filament. The former pair 
authors point out that switching time second possible 
negative resistance region. was remarked above the thermal time 


| 
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will not small this switching time when the parameter large com- 
pared with unity and, indeed, this phenomenon. 

Even devices with long thermal time constants those used frequencies 
high enough for thermal effects unimportant the behavior the low- 
frequency differential resistance given equation (1) still importance 
determines the d-c. and low-frequency stability the device conjunction 
with the associated d-c. bias circuit and must heeded obtain monostable 
operation free from low-frequency oscillations. 

the discussion above, the only inertial process considered was thermal and 
the isothermal admittance was assumed real and frequency independent 
and equal the slope the isothermal characteristic the operating point. 

The treatment can generalized replacing which 
the small-amplitude admittance and fixed temperature The 
impedance the device including thermal effects then becomes 


which behaves for frequencies sufficiently high for and 
The condition for negative resistance occur becomes 


where 


high frequencies where the isothermal admittance may become complex 
considered above also necessary examine the parameter which 
determined Here denotes the total current through the device 
and real only the current variation constant voltage cophasal with 
the temperature variation producing it. This would not the case, for instance, 
when the capacitive component current temperature modulated varia- 
tion the permittivity the material the geometry space-charge 
region. Hence when becomes complex high frequency account must 
taken also becoming complex though with phase angle which will 
general differ magnitude and even sign from that 

complex and furthermore the temperature coefficient 
negative occurs ferroelectric crystals near the Curie point. For instance, 
barium titanate can made semiconducting adding cation impurities 
suitable valency (Saburi 1959), and the resulting material can have room 
temperature resistivity the range ohm cm. the Curie point 
increase temperature may produce 100-fold increase resistivity 
Thermistors constructed from such materials would have very large negative 
values for example, the applied power causes 10° rise temperature 
above the ambient —10. 

such ferroelectric material near the Curie point the temperature co- 
efficient the permittivity would large and thus would significantly 
complex. For have 
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with both and negative. The comparison 
with the thermal time constant determines the importance the ferro- 
electric behavior modifying the dispersive properties the impedance. 


APPENDIX 
THERMAL TURNOVER THE TWO-ACTIVATION-ENERGY MODEL 

germanium diode the consistently observed decrease turnover power 
with increasing ambient temperature proves exacting the construction 
satisfactory phenomenological model for the characteristics. The suggestion 
(Cutler 1957) that dependence based 
involving two activation energies here shown yield the correct behavior 
but with interpretation the experimentally derived parameters which 
differs from that commonly used. 


Consider current the form 


temperature constant voltage when plotted logarithmically against 
will display two asymptotic linear regions with slope low temperature 
and high temperature and intersecting temperature where the 
two components current are equal: 


Now under steady-state operation and for linear thermal conduction the 
device temperature and ambient temperature are related 
where the thermal conductance between the current-controlling region 


the device and the ambient. 
The turnover condition independent and given 


where suffix denotes the values turnover and the slope function given 


For turnover occur 47, necessary. 

The function has the value low temperatures, and high tempera- 
tures, and where close approximation reaches its 
maximum value. 

The dependence the turnover power the ambient tempera- 


ture found from 
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germanium diodes this quantity found experimentally negative 
over extensive range temperature. 
For this possible 


The most favorable temperature for this realized close 
for which the inequality becomes 


This condition does not depend the thermal conductance and deter- 
mined only the two ratios and a/b. 
Introducing the parameter 


which must greater than zero for the above inequality hold, have the 
following results the condition 


Thus the intercept the tangent the 7,) plot the axis 

this intercept temperature which the literature (Benzer 1949; 
Yamaguchi and Hamakawa 1959a, has been identified erroneously with the 
contact temperature with the temperature which the onset intrinsic 
conduction occurs. fact the intercept larger than the contact temperature 
the fraction 1/C (which may typically lie between 0.1 and 0.3) and has 
relation intrinsic conductivity. Similarly the slope usually much 
less than whereas often assumed provide direct measurement the 
thermal conductance. 

Typical values for germanium diode are: 


2000°, giving 
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will noted that the thermal conductance way affects the 
temperature relations only appears expressions involving the power. 
Furthermore the value not great influence the results and could 
zero (corresponding current component independent temperature) 
and yield the same qualitative features for non-zero value. 
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THE ELECTRICAL RESISTANCE DILUTE MAGNESIUM 
AND ALUMINUM ALLOYS LOW TEMPERATURES! 


ABSTRACT 


Experimental data the influence type and concentration impurity 
the low-temperature resistance the metals, magnesium and aluminum, 
presented. particular, comparison made the electrical behavior alloys 
formed using transition element impurities the divalent metal magnesium, 
and the trivalent metal aluminum with the behavior similar impurities 
some the noble metals. 


INTRODUCTION 


The electrical resistance minimum metals was discovered years ago 
gold Haas and van den Berg (1934), and since that discovery 
considerable number pure metals and alloys have been investigated for 
similar low-temperature anomalies. Most these experimental studies have 
involved the noble metals and their alloys. present some order beginning 
emerge from systematic investigations the influence type and con- 
centration impurity the temperature occurrence and relative size 
the resistive anomaly.* 

Electrical resistance minima magnesium were first reported Meissner 
and Voigt (1930) and later work MacDonald and Mendelssohn (1950), 
Garfunkel, Dunnington, (1950), Thomas and Mendoza (1952), and 
Spohr and Webber (1957) confirmed these observations. Reports elec- 
trical resistance minimum occurring aluminum (Keesom 1933) were not 
confirmed later measurements made Boorse and Niewodniczanski (1936), 
Garfunkel, Dunnington, and Serin (1950), and later Thomas and Mendoza 
(1952). Although low-temperature resistive anomalies have been reported 
magnesium and aluminum, detailed study the influence type and 
concentration impurity the minimum has been carried out. Interest 
this laboratory the magnetic properties magnesium and aluminum 
prompted investigation the electrical properties number dilute 
alloys magnesium and aluminum low temperatures and wish report 
briefly these measurements. 
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EXPERIMENTAL METHOD AND SPECIMEN PREPARATION 


2.1 Apparatus 
The low-temperature cryostat similar the helium desorption unit 


described Rose-Innes and Broom (1956) and will not described here. 
With freshly activated charcoal the desorption unit, temperature run 
from 4.2 30° could made very slowly, the run usually taking from 
hours. This meant that number resistance specimens could run 
one time. All temperatures above 4.2° were measured with Simon-type 
gas thermometer. Relative temperatures could measured with accuracy 
better than 0.3% whereas the absolute temperature the higher tem- 
perature regions should within 1%. Between and 4.2° the temperature 
was determined using the known vapor-pressure scale. Temperature agreement 
between the gas thermometer and vapor-pressure measurements was better 
than 0.2%. 

current-potential method was used measure the resistance the 
samples. Pressed contacts were used and were either screw- knife-edge 
design. The voltages were measured with galvanometer amplifier described 
MacDonald (1956). Scatter the experimental points (due electrical 
noise and stray thermal voltages) permitted the resistance ratios the 
aluminum alloys determined with accuracy least whereas 
the magnesium alloys which exhibited large anomalies were not measured 


better than 1%. 


2.2 Magnesium Specimen Preparation 

The alloys were made Dow Chemical Company metallurgical laboratories 
and supplied either chill-cast ingots 0.25-in. rods. The resistance speci- 
mens were made hot-rolling* the ingots 400° into strips approximately 
thick and cutting these strips into wires rectangular section. 
After fabrication, the samples were annealed 450° for hours 
atmosphere helium gas. Table lists the alloys investigated their 
resistance ratios All the alloys show resistance ratio 
higher than the pure stock that, not all, least some the impurity 
solution. That alloys had been made the case magnesium—manganese 
and systems can seen Table where the resistance 
ratios show approximate linear increase with nominal impurity content. 
Hot-rolling sample 732 produced extremely large impurity gradients 
the samples that values for this specimen given Table are thought 
averages. The room-temperature resistivity alloys 702, 730, 731 are 
4.62 ohm-cm, 4.67 ohm-cm, and 4.82 ohm-cm, respectively, 
and are good agreement with the published values Salkovitz, Schindler, 
and Kammer (1957) for these concentrations and type solute. Micrographic 
examination the alloys showed one-phase system with minimum 


oxide present. 


*We are indebted the Department Mines, Ottawa, for carrying out the hot-rolling. 
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Aluminum Alloy Preparation 

The aluminum alloys were prepared the Aluminum Company Canada 
and spectrographic analysis these alloys shown Table II. The alu- 
minum samples were cold-rolled the same size the magnesium alloys 
and usually given strain-relieving anneal for hours 300° atmos- 
phere helium prior making measurements. Table also shows the 
resistance ratios for the aluminum alloys function nominal impurity 
content. should pointed out that some difficulty was experienced 
keeping the iron solution. The maximum solubility iron aluminum 
the samples were quenched from 620° after having received homoge- 
nizing anneal this temperature for approximately hours. The purpose 
this anneal was try and remove the large impurity gradients that existed 
the original ingots. estimate the gradient that existed the original 
ingot may seen Table where measurements the resistance ratio 


TABLE 
Spectrographic analyses and electrical properties aluminum alloys 

Alloy Major 


16.4X 


are given for each end the ingot which contains nominal impurity 
content 


EXPERIMENTAL RESULTS 

Electrical Resistance Magnestum Alloys 

series temperature-resistance runs magnesium are shown Fig. 
Sample 728 the pure material and exhibited resistance minimum. Sample 
727-1 iron alloy which presumably has too small minimum 
detected with the present method measurement. Alloy 757-3 has both 
iron and present and the minimum presumably due the presence 
iron. None the alloys listed Table containing major impurities 
tin, antimony exhibited resistance minimum. The rest the allovs 
are magnesium manganese The depth the minimum 
magnesium previous authors (see above). Both the depth the minima 
and the temperature which the minima with increasing 


phase 


are indebted Mr. Marchand the Aluminum Laboratories for making 
diagram aluminum iron available 
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the impurity concentration, Fig. shows graph and the normalized 
depth the minima function impurity content. The behavior both 
these parameters similar that described for certain copper alloys 
MacDonald and Pearson (1955). 


3.2 Electrical Resistance Aluminum Alloys 

Figure shows the temperature dependence relative resistance some 
the and aluminum-iron alloys described Table II. 
None these alloys exhibit low-temperature resistive anomaly. 
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3.3 Room Temperature Resistivity and Residual Resistance the and 
Aluminum Alloys 
Mathiessen’s rule states that the total resistivity alloy considered 
consist thermal component which the same the pure solvent 
(pm) and temperature independent contribution due scattering from 
impurity ions then the total resistivity can written 


Although the statement limited value, suggested that since all 
the resistance temperature curves approach parallelism above the region 
where impurity scattering dominates the present alloy, systems may said 
obey the Mathiessen rule approximately. 

Nordheim (1931), assuming the Mathiessen rule, and using first-order 
perturbation calculation, has shown that the electron concentration and 
structural properties the base metal are not modified then for 
low concentrations solute (and hence should increase linearly 
with concentration. Using the data Tables and may seen that 
the Nordheim rule not obeyed well for these dilute allovs. has been 
discussed previously MacDonald and Pearson (1955), the nominal com- 
position dilute alloys often bears very little relation the amount 
impurity solution and for this reason the resistance ratio more reliable 
indication the actual amount impurity dissolved. 

Values the room-temperature resistivity for pure magnesium (alloy 728) 
and pure aluminum (alloy are ohm-cm, and 
ohm-cm respectively. Both the values room-temperature resistivity and 
resistance ratio for the pure metals aluminum and magnesium are good 
agreement with the most recent values reported Thomas and Mendoza 
(1952) and Salkovitz (1957). 


CONCLUSION 


Five general conclusions may drawn from the alloy systems mag- 
nesium and aluminum studied far: 

(i) Nosample nominally pure magnesium showed resistance minimum 
within the accuracy the present measurements. 

Only the transition elements, manganese and iron, appear induce 
magnesium, manganese being the more effective the two. 
Solutes such antimony, and tin not cause minimum. 

resistance maximum that manganese induces copper, gold, and 
silver (see article, 1956) temperatures lower than the minimum 

(iv) The concentration and temperature dependence the minimum 
the exhibits behavior similar that the 
copper alloys with group solutes gallium, indium, germanium, 
studied MacDonald and Pearson (1955). seems reasonable conclude 
therefore that the general behavior the resistance minimum divalent 
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metal such magnesium similar the behavior the minimum 
monovalent metal such copper containing group solute. 

(v) The present results made more complete set alloys than those 
previously studied are agreement with the conclusion Thomas and 
Mendoza (1952) that the transition elements iron and manganese not 
induce minimum aluminum. 
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THE ASSOCIATION VISIBLE AURORAL FORMS WITH 
RADAR ECHOES! 


Lyon 
ABSTRACT 


peak 48.2-Mc/sec echo occurrence observed Saskatoon corresponding 
time the period breakup quiet arcs into active rayed structures. This 
also the time most frequent occurrence characteristic forms 
the aurora. If, Gartlein suggests, the forms are formed instabilities 
sheet beam then the primary particles are positively charged. 


INTRODUCTION 


Many observers have noted close statistical correlation between the 
occurrence very high frequency radar echoes and visible aurora. Little, 
Rayton, and Roof (1956) have summarized much the previous work and 
concluded that the radio wave scatterers are closely associated space and 
time with visual aurora, but that because aspect sensitivity, only aurora 
occurring low elevation angles toward the north seen from the radar 
station may expected give echoes. general, auroral forms elevation 
angles below about 15° will appear consist quiet arcs diffuse forms, 
whereas the same aurora when viewed near the zenith may exhibit more 
structure and motion. Thus because the radio-observing site must some 
500 from the scattering region, zenith observations the visible form 
corresponding given radar echo are rarely obtained and few workers 
have attempted correlate radio echoes with specific visible forms. Currie, 
Forsyth, and Vawter (1953) Saskatoon, using simultaneous camera and 
radar observations from the same site, noted that echoes were usually ob- 
served from moderately bright aurora showing some structure and that 
structure was apparently more important criterion for echo occurrence than 
luminosity. 

During recent radar investigation auroral motions (Lyon and Kavadas 
1958) carried out Saskatoon (latitude 52.1° N., longitude 106.6° W.), 
excellent statistical correlation was obtained between echo occurrence 
and the occurrence visible aurora recorded all-sky camera 
situated Flin Flon (latitude longitude 101.8° W.), approxi- 
mately under the scattering region. This paper concerns attempt detect 
more detailed correlation between echo occurrence and visible auroral 
forms observed the Flin Flon camera. 


OPERATIONAL DETAILS 
The radar frequency was 48.2 Mc/sec and the radar was operation during 
the period January April 28, 1958. The equipment has been described 
previous paper (Lyon and Kavadas 1958). 
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detailed study was possible only when certain conditions were all fulfilled, 
namely that (1) the radar was operating, (2) the camera was not obscured 
cloud snow, (3) visible aurora occurred, and (4) the picture quality 
was such permit identification the auroral form. These conditions 
were all fulfilled nights during the period under investigation. For 
analysis each night was divided into 10-minute periods. all there were 910 
periods which visible aurora was present. This was considered 
sufficient number allow statistical treatment. 

The useful field view the camera covered the ranges 
between 350 and 650 from the radar Saskatoon. Few echoes were 
received from ranges less than 350 km. Some echoes were received from 
ranges greater than 650 and these cases the form the visible aurora 
beyond the camera was taken the same that within the camera field 
view. Most the echoes were received from ranges within the camera 
field. attempt was made correlate the visible aurora and the echoes 
exactly range usually the form present the exact echo range was also 
present lesser and greater ranges the same time. 

RESULTS 
Echo-Visible Form Correlation 

general, the pattern the visible aurora for any one night followed 
that described Heppner (1954). The form before local midnight was usually 
quiet arc and the form after midnight was usually active form, often 
exhibiting rayed structure. During the period around midnight which 
the breakup arcs into rayed structures occurred, the intermediate form 
which most often developed was that are with spiral dis- 
continuity. typical camera record arc with “curl” shown Fig. 

For the purpose analysis the visible aurora was divided into three main 
groups, namely: 

(a) quiet arcs, 

(6) ares with “curl” spiral discontinuity, 

(c) active rayed structures. 
this basis the form present any one 10-minute period was recorded. The 
radar records were then analyzed using the same 10-minute periods had 
been used for the camera records. echo occurring any one period was 
associated with the visible form the same period displayed Table 


TABLE 
Echo-form correlation 
No. 10-minute periods 
Percentage 


Visible 

aurora Echoes echoes giving echoes 
Quiet arcs 365 347 
Curls 174 120 
Rayed structures 371 319 
910 124 786 


Total 


PLATE 


Some Flin Flon all-sky camera records, showing ares with structures, for 
the night March 1958. 


Note that prints are taken from positive film that the auroral form appears dark. 


> 
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Thus seen that while 14% all visible aurora could associated with 
echo, only the quiet arcs could associated with echoes, the figures 
for and active rayed structures being 31% and 14% respectively. 
Further, may noted that 

(a) quiet arcs occurred 40% the time but accounted for 14% the echoes; 

occurred 19% the time but accounted for 42% the echoes; 

(c) rayed structures occurred 41% the time and accounted for 41% 


the echoes. 


Echo and Form Diurnal Variation 

Using the same 10-minute periods, the total number periods which 
visible aurora and/or echo occurred any one heur throughout the night 
was computed and plotted against local time shown Fig. This diurnal 


QUIET ARCS 


wo 
ACTIVE 
RAYED 
a 
uw 


LOCAL TIME 


Fic. Diurnal variation occurrence visible auroral forms and echoes. 


variation illustrates clearly the general pattern visible aurora: quiet arcs 
followed and active rayed structures. The similarity between 
the occurrence curve and the echo occurrence curve striking. 


The Direction the 
The appearance arc developing spiral structure bears 


marked resemblance the characteristic instability profile sheet beam 
charged particles under the influence longitudinal magnetic field. 


/ 
/ 
/ 
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This fact has also been noted Gartlein (1959), who has examined qualita- 
tively the expected appearance auroral arc assuming the causative agent 
beam protons moving downward along the lines force the 
earth’s magnetic field. states that the beam not quite uniform, the 
more intense parts will curve into shape figure and the intensity 
raised these curves will into spirals which, when viewed from the ground 
looking upward, spiral inward clockwise direction. These spirals finally 
break into narrow rays. Further, this pattern not solely intensity 
effect but will also occur weak beam after more extended travel. Gartlein 
notes that the southern hemisphere the spirals observed spiral inward 
anticlockwise direction would expected. concludes that aurora, from 
the time arc formation until well after the time ray formation, the 
expected effect sheet protons moving along magnetic field. 

All the spiral forms examined this investigation spiraled inward 
clockwise direction when viewed from the ground, thus indicating that the 
primary causative agents for such forms are positively charged particles. 
Examination the records some the visual observations aurora made 
Chesterfield, N.W.T., during the 1932-33 (made available Dr. 
Currie, Institute Upper Atmospheric Physics, University Saskat- 
chewan) again reveal that when such forms were observed they spiraled 
the same direction. Further, Hector (1959), who has examined many all-sky 
camera records taken from Canadian stations during the 1957-58 I.G.Y. 
reports over 100 cases such forms all without exception spiraling the 
same direction. 

SUMMARY AND DISCUSSION 

While the results reported here have been obtained for one station and one 
frequency (48.2 Mc/sec) only, would appear that 

(a) the probability echo from quiet arc small; 

the most probable form which will give echo 

(c) active rayed forms occur more frequently than there 
equal probability any one echo coming from either active 
rayed form. 

There remarkable similarity between the diurnal variation echoes 
and that (Fig. 2). The spiral form the one most often 
seen during the period breakup quiet arcs and the formation rayed 
structures and evidently intermediate form. Thus, there would seem 
very real case for associating echoes with the period breakup quiet 
arcs into active rayed structures. this experiment the receiving antenna 
was fixed azimuth direction 45° east geomagnetic north. Forsyth 
(1959) has suggested that spiral form presents greater volume 
ionization pulse electromagnetic radiation than does quiet arc. 
Thus, might expect receive more echoes from forms, but the 
same token had the radar antenna been directed along the geomagnetic 
meridian greater number echoes might have been obtained from quiet 
arcs than was obtained with the antenna directed 45° the meridian. 
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All the available evidence from the direction the ‘‘curls’’ and spirals 
points positively charged incoming particles. 

Gartlein (1959) has suggested that particles entering the northern por- 
tion the auroral zone (in the northern hemisphere) traverse longer paths 
than those entering the southern portion the zone, then auroras observed 
the northern part should more frequently composed rayed structures. 
Possibly more should seen also, though the longer path travelled 
the particles the structure may have developed before the excitation 
the visible aurora takes place. interesting note that Collins and 
Forsyth (1959), bistatic radio investigation auroral ionization over 
different paths, report that for type events (signal enhancements related 
auroral ionization) occurrence increases with increasing latitude, least 
the auroral zone. This just what might expected echoes are pre- 
dominantly associated with rayed structures and the occurrence 
such structures increases with increasing latitude also. 

Forsyth and Green (1959) report that the peak the diurnal variation 
occurrence type events occurs progressively earlier with increasing 
latitude. possible that the average time breakup quiet arcs into 
rayed structures also occurs progressively earlier with increasing latitude. 
hoped report further this point later paper. 
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INFRARED TEMPERATURE AND LINE STRENGTH 
MEASUREMENTS CARBON MONOXIDE EXCITED 
RADIO-FREQUENCY DISCHARGE! 


ABSTRACT 


Relative line intensities the sequence bands carbon monoxide 
observed emission from radio-frequency discharge have been measured. 
enhancement the branch lines and diminution the branch lines due 
interaction has been observed agreement with theoretical 
predictions. Vibrational and rotational temperatures have been calculated and 
when allowance made for the interaction these prove 
equal within experimental error indicating the existence thermal equilibrium 
the discharge. 


INTRODUCTION 


The problem the relative strengths individual lines 
rotation band has been the subject number recent papers. The theory 
the effect the interaction vibration and rotation the line strengths 
has been developed Herman and Wallis (1955) and extended other 
authors. The predicted effect has been observed means absorption 
experiments the fundamental bands (Benedict, Herman, 
Moore, and Silverman 1956, 1957), (Kuipers 1958), and DCN (Eaton 
1959). the case the effect the fundamental predicted too 
small for observation and this also agreement with experiment (Herman, 
Silverman, and Wallis 1954; Eaton and Thompson 1959). 

The results referred above have all been obtained from absorption spectra 
but there are some advantages gained carrying out measurements 
emission rather than absorption. particular the effect the finite slit 
width such that the relative line strengths can obtained more easily 
and directly from the actual measurements and the elevated temperature 
the gas allows the extension the measurements higher values and 
vibrationally excited states. most molecules the discrepancy between 
and line strengths becomes pronounced only high values. This type 
measurement also provides method probing the closely related problem 
the existence otherwise thermal equilibrium among the various rota- 
tional and vibrational energy levels gas flame discharge. 

The present measurements have been carried out the sequence 
bands carbon monoxide. This molecule has very small dipole moment 
that the parameter also small and the nature the vibra- 
interaction rather different from that molecules such 
and for which The theoretical treatment Herman and 
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Wallis (1955) and Herman, Rothery, and Rubin (1958) predicts that the 
interaction effect the band will become appreciable high values and 
that will the opposite sense that observed more polar diatomic 
molecules. This prediction has now been experimentally verified. The intensity 
distribution these bands also some astrophysical interest since they 
have been observed the sun and used determine solar temperatures 
(Goldberg and Miiller 1953). 


EXPERIMENTAL 


The sequence carbon monoxide has been previously observed 
emission from flames Plyler, Benedict, and Silverman (1952) and under 
lower resolution from R.F. discharge Wilkinson, Ford, and Price (1955). 
The discharge tube used for the present work was pyrex, meter long and 
around diameter. was equipped with system multireflection 
mirrors the type first suggested White (1942) for use absorption work 
and later adapted for emission experiments Welsh, Cumming, and Stans- 
burg (1951). was found that the optimum signal was obtained with path 
length about meters. Power was provided 5-kw oscillator operating 
Mc/sec. The output from this oscillator could controlled adjusting 
the position coupling coil. Power was transferred the discharge tube 
either placing the latter inside tuned coil means aluminum foil 
electrodes. The tube was air-cooled but still tended overheat and crack 
when the foil electrodes were use. For this reason the coil method coupling 
was used spite the rather greater efficiency the electrodes. Cylinder 
carbon monoxide was streamed through the tube constant pressure 
mercury. Under these conditions stable discharge giving re- 
producible line intensities could maintained for periods hour. 
The vacuum infrared spectrometer has been described previously (Douglas 
and Sharma 1953). The detector was lead sulphide cell cooled dry 
acetone temperature and the slit widths employed corresponded resolu- 
tion about 0.2 The amplifier time constant was 0.3 second and the 
scanning speed was sufficiently slow allow maximum throw the 
strongest lines. Readjustment the fore-prism was made every 

The variation the sensitivity the instrument over the frequency range 
3850 4400 was obtained calibration with tungsten filament 
bulb using the tungsten emissivity curves Vos (1954) and the black- 
body radiation tables Lowan and Blanch (1940). The temperature the 
filament was measured with optical pyrometer and was verified calcu- 
lation that the corrections applied are insensitive possible errors 
this temperature. 

RESULTS 

The observed spectrum consisted the sequence containing the 
the and bands were also observed but not sufficient number 
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Line 2-0 3-1 
Pi 19.3 
P2 21.0 
Ps 

Ps 

Ps 

Ps 43.1 

Ps 50.5 
Ps 25.7 

Pio 28.3 

Pu 32.6 66.6 
Pr 35.2 

Pis 36.8 78.3 
Pu 72.7 
Pw 85.5 
Pus 45.1 81.6 
Pu 80.3 
Pris 43.0 

Pis 43.0 

Px 

Pu 50.7 

Px 78.1 
Pes 62.7 78.9 
Px 61.4 74.3 
Ps 60.9 76.9 

Pu 

77.2 
Po 74.9 
Px 72.7 
Pa 71.9 
Px 62.0 73.1 

Ps 59.8 67.4 
Pu 

Px 

Px 55.5 69.8 
Px 65.0 
48.5 57.0 
Pas 57.1 
Po 51.9 
Pa 46.0 651.2 
Pe 45.5 45.5 
Pa 

Pu 

Pw 36.4 37.2 
Ps 34.6 37.4 
Pa 35.5 33.5 
Ps 30.6 25.6 
Po 27.6 

Pw 24.8 23.0 
Px 29.6 
Ps 22.3 

Ps 19.0 21.9 
Pu 15.0 

Ps 13.2 
Pu 

Ps 15.2 
Ps 14.1 
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2-0 3-1 
13.2 
24.7 
35.4 
15.4 54.2 
59.6 
63.1 
72.9 

26.5 
39.0 90.1 
38.9 97.3 
40.6 94.6 
41.4 101.7 
41.3 104.4 
44.6 108.8 
53.9 112.4 
54.6 110.6 
55.5 112.4 
56.3 112.3 
61.3 116.8 
64.6 118.5 
67.9 111.3 
66.9 118.3 
130.8 
75.8 122.7 
112.0 
76.6 117.3 
7.5 116.4 
74.8 114.7 
74.8 109.3 
71.4 103.0 
74.6 101.2 

73.8 
70.4 94.1 
72.8 99.4 

71.1 
65.3 89.5 
64.5 86.8 
63.6 83.2 


96. 


89 


to 


67-4 7-5 
18.3 
18.4 
17.8 
10.2 
25.9 17.2 
26.8 18.7 
28.5 18.8 
27.8 18.9 
31.0 19.0 
33.6 
33. 
18.6 
33.3 
18.8 
34.4 
36.2 
36.3 17.6 
40.5 16.9 
28.3 18.5 
17.8 
27.6 
15.6 
26.9 
27.8 17.2 
14.9 
26.2 
23.7 
24.6 
24.6 


11.6 
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30.7 
15.7 
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al 21.6 10.1 Ris 89.0 
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| Ros mi 
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with sufficient intensity for the present type measurement. portion 
the spectrum shown Fig. indicate the degree resolution. About 
60% the total number lines was sufficiently free from overlap 
allow reliable measurement. The measured relative intensities these lines, 
corrected for changes instrumental sensitivity, are given Table The 
absolute values are arbitrary. 

The amount light emerging from the exit slit the spectrometer when 
set the center emission line given expression the type 


0 


the integrated intensity the mth line and and are functions 
describing the line shape and spectrometer slit function respectively. They 
involve the line width parameter and the spectral slit width 
Equation (1) much simpler than the corresponding expression for absorp- 
tion that directly proportional Under the present experimental 
conditions the line width arising from Doppler broadening will around 
0.01 and that from pressure broadening <0.001 (Eaton and Thomp- 
son 1959). Natural line width negligible the infrared. Since the slit width 
~0.2 the slit function will therefore approximately constant 
over the range frequencies for which appreciably large and can there- 
fore removed from the integral equation (1). Variation slit width with 
frequency therefore absorbed into the instrumental sensitivity correction. 
Further, when has been removed the remaining integral must the 
definition integrated line intensity constant independent the line 
width The values given Table will therefore used direct measure 
the line intensities. 

There source error emission intensity measurements which should 
also considered, namely that arising from self-absorption. Errors from this 
source will greatest lines arising from transitions heavily populated 
lower states, i.e. the lower-lying vibrational energy levels and values. 
Fortunately the high values which are most significant the present 
case. Preliminary experiments showed that self-absorption from gas the 
cooler parts the tube was appreciable for the band and revealed 
itself distinct curvature the temperature plot. The measurements from 
this band have not therefore been used any the subsequent calculations. 
There is, however, indication significant effect the higher bands and 
the results from these latter are thought reliable. 


DISCUSSION 
expression the type 


The function F(m) (where for the branch and J+1 for the 


branch) introduced allow for the interaction. The 
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other symbols have their usual significance. F(m) unity the band origin 
but becomes increasingly different from unity higher values. F(m) 
were unity all values, i.e. interaction, the rotational 
temperature could found plotting log against 
and obtaining the slope. If, however, the values for the corresponding lines 
the and branches are not equal the points for the two branches will 
fall separate curves. The separation the curves will give the ratio the 
factors given value but the method temperature determination 
will fail. Figure shows this plot for the band. The results for the 
other bands are very similar. The frequencies and energy levels were calcu- 


"T RBranch" = 3570°K 


I 


log 


R Branch 


ie) 1000 2000 3000 4000 5000 6000 


lated from the molecular constants given Plyler, Benedict, and Silverman 
(1952). apparent that the factors are appreciably different the two 
branches, being largest the branch which accord with the theoretical 
treatment Herman, Rothery, and Rubin (1958). Straight lines drawn 
through the points for the two branches gave approximate apparent tem- 
peratures 3570° and 2780° and since the theory predicts that F(m) 
will greater than unity and less than unity roughly equal amounts 
the true temperature will between these values and may estimated 
around 3175° not possible, however, obtain either the temperature 
the factors independently from these plots. 

vibrational temperature may obtained from the relative intensities 
corresponding lines the different bands the sequence. many pairs 
possible from both and branches were chosen and from the average 
ratios the intensities the other bands relative the band were 
calculated. also necessary know the vibrational transition probabilities 
and the values given Herman (1952) were taken. These are stated 
uncertain the extent 10%. Figure shows the plot log 
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against (vib.). This leads vibrational temperature 3120° Bearing 
mind the approximations involved both cases this value probably 
sufficiently close the estimated rotational temperature 3175° in- 
dicate the existence thermal equilibrium between the vibrational and 
rotational energies. interesting note that from the flame spectrum 


2.5 


2.0 


Vibrational Temp. = 3120°K 


5000 10,000 15,000 


Vibrational temperature plot for the sequence bands. 


Plyler, Benedict, and Silverman obtained vibrational temperature 
2360° and rotational temperature 2800° This latter temperature, 
though, was obtained only from branch lines and the apparent discrepancy 
very similar the corresponding quantity obtained from the present dis- 
charge spectrum (3570° rotational, 3120° vibrational). correct allow- 
ance for the interaction could well produce agreement 
between the two temperatures. would seem from this that the conditions 
energy transfer carbon monoxide flames and discharges are quite similar 
and that both cases there sufficiently close approximation thermal 
equilibrium give physical meaning the temperature values obtained. 

order further substantiate this conclusion the vibrational value 
3120° was taken rotational temperature and used obtain the 
values from equation (2). The results are illustrated Fig. points from the 
and bands being included. The solid curves are obtained from 
the theory Herman, Rothery, and Rubin. The theoretically calculated 
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Pekeris 
1.3 


Band fe) gt 


Anharmonic 


Fic. experimental and theoretical values. 


values depend the form the potential field assumed and curves for both 
anharmonic oscillator and for Pekeris oscillator are included. The 
values also depend whether the signs the dipole moment and its 
derivative for are the same opposite and the curves labelled 
and represent these alternatives. can seen from Fig. the experi- 
mental factors agree well with the theoretical values within the experi- 
mental error but this error too large enable choice made between 
the different types oscillator between the alternative signs The 
factors for the higher bands show systematic deviation from those for the 
lower bands. This also would seem agreement with theoretical 
expectations. 

may concluded that the Herman and Wallis theory gives satisfactory 
account the effects interaction for the present case 
where small and the anharmonic effects are predominant well for 
cases previously studied where large. allowance for this effect must 
made meaningful temperatures are obtained from intensity measure- 
ments and the present case the results show that thermal equilibrium 
between the vibrational and rotational energies has been attained. The results 
obtained also indicate that some cases infrared intensities can with profit 
studied emission rather than absorption. This may prove par- 
ticularly true the case unstable radicals such OH. 
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TABLES THE MODIFIED BESSEL FUNCTIONS THE 
SECOND KIND FOR PARTICULAR TYPES 


ABSTRACT 


Tables the modified Bessel functions the second kind for arguments 

the form and where integer and for all integral and 

half-integral orders the range are presented. For all but small range 

the argument the tables can relied upon the figures given; the 

omitted range error unity the 10th place may occur. The tables should 

useful for the computation many lattice sums arising the theoretical 
calculation crystal properties. 


INTRODUCTION 

For the theoretical calculation the properties ionic crystals impor- 
tant able evaluate, reasonably accurately, various summations over 
the lattice sites, where the quantities summed are usually potential functions 
involving inverse powers the distances from the ions the point which 
the potential required. functions may represent the relatively 
long-range attractive interactions arising from the electrostatic and dispersive 
forces within the crystal, not surprising that often their sums not 
converge very rapidly. Occasionally for inverse powers greater than about 
three four, summations this type have been carried out direct addition 
terms for large number lattice points, the summation being completed 
integration over the remainder the lattice (Lennard-Jones and Dent 

1928; Shuttleworth 1949). Alternatively the summation may effected 
transforming the sum more rapidly converging form and summing the 
resulting series directly (see, for instance, Madelung 1918; and van der Hoff 
and Benson 1953). The latter approach has been used fairly recently 
papers dealing with the Coulomb and van der Waals potentials associated 
with the problem the adsorption polarizable atoms ions onto the 
surface infinite hemicrystal (Hove and Krumhansl 1953; Hartman 1956; 
Durham and Soderberg 1958) and for the calculation edge energies (Schreiber 
and Benson 1955) and surface energies (van Zeggeren and Benson 1957; Ben- 
son, Balk, and White 1959; Benson, Dempsey, and Balk, unpublished). All 
these applications have the common ground that the series resulting from 
the various transformations involve the modified Bessel functions the 
second kind, although frequently other transcendental functions also 
arise. For the definition used here see Watson (1948). 

the work originating from this laboratory (references above Benson 
and co-workers) has been found that the arguments the Bessel functions 
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which have occurred far have integral and half-integral orders the range 
0-8 and that the arguments may all expressed one other the general 


forms: 


(1) 


problems involving distortion the lattice surface, the potential energies 
required usually have evaluated points off the regular lattice array. 
Since the displacements from the regular sites are, general, small, the 
summations required may, these cases, expanded Taylor series about 
zero displacement. This allows the original summation expressed 
linear combination number summations involving functions the 
form (1), but with orders which, general, are higher than those the 
original summation. For problems concerned with the interaction ions 
situated arbitrary distance above infinite hemicrystal (as the 
references adsorption problems given above) the argument the Bessel 
functions will usually include many these cases should also 
possible remove the variable from the argument suitable series 
expansion. 

order facilitate the evaluation lattice summations our own 
work, tables the modified Bessel functions for arguments the forms (1) 
have been computed for fairly wide range orders and arguments. The 
purpose this paper present these the hope that they may use 


other workers. 


CALCULATION THE TABLES 

The computation the tables was carried out entirely Bendix 15D 
digital computer operating double-precision floating point mode. Although 
this gave figures the final results, the entries the tables were eventually 
rounded off figures. simplify the presentation, the tables have been 
left the floating point notation, which the two figures the left the 
point indicate, excess basis, the power which the remaining 
decimal part must multiplied. This notation should evident from the 
following examples: 


46.123 
50.123 0.123 


The basis for most the computation was the asymptotic expansion for 
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MODIFIE 


TABLES 


SON 


N 


AND BE 


DEMPSEY 


TEEMEZQOST "EN 
LOGZOLESOT* EN 
QLOLentegt 
STSSTHECOe 
nSStnngste eh 
TEOLESLETE 
QTZEG2ESHE “EN 
QnE22LoeLe en 
SSEMEQTEON 
TTE9LLESEN 
ZTLS62ETSS 
6n2n0SSL6S 
ngeengt ng’ en 
99SS902TOL*En 
99620£009L° 
E66960THEQ 
T296SQOO0L6* 
QEOONTESOT 
ZELSLLENTT 
62LOELEHeT 
LLTETQOSET 
LETONERIAT Th 
096980865 T* 
TEELSESGELT 
EQELIBLEQT 
LSQT9ZEQON * Hh 
T606eLENce * Hh 
Hh 
“Hh 
* 
LETELTQQTE 
LSOQ6TEQHE 
ZTLASHLOBE 
S620SON9TH 
* 
LEEGLOEGES * 
699081959 *trh 
SETLLS9OEL * 
2TELIOSTEL 
QTROOLTE9Q 


GLOLSHLONG * 
THHLELLTOT 
TLAQ9EN6OT 
6S9HTTLLIT 
QETEQOESET “En 
EENOOLORST 
QTLS9EZOLT "En 
QE2OQELLET* EH 
OLESQLE6eC * En 
OSTETHTQHe 
QLLLOSHLEL 
ETEONSENEE 
ETS9OTSTON 
SLHS6ES66n 
6Q25060NS 
ZEBSETESQS 
"En 
En 
neonTenns6 
LLOBLENECT* Hh 
QTLOORCEET 
QLTOSHGLST* rh 
6Q6TOSETLT 
QTEELEELET 
* Hh 
QnSLEQLece 
* 
"nh 
* 


nSTHSLnZeg* en 
$9L6092999 
QETSHEQOLL*2h 
SozeL_tegeg*en 
TLL2ZL206Q° 2h 
NSZLQETLSE 
Tetéckézort 
"En 
QLOBTSTETT’ En 
En 
QT69LNZO9T "En 
TSQ9SOLELT*En 
th 
negLoeggte* 
En 
Ot nagtese 
Lesgtonnée* 
ELLHeNETLE 
LEQOTEEZON* En 
6nTOELNSEN 
En 
SE999C90TS 
TE999SeESS 
HTTOZR96ES 
9596292059 "En 
SLEEQEESOL* En 
LEBTETSS9L* ER 
En 
OOOLOTLE06* Eh 
S9SSLEQOR6 
SO6TEL99OT “Hh 
6296696STT 
QneezQtget 
669S696e9T 
NZEQHLOLLT 
OTISLLLEGET * 
LEeTLEGOES * 
£960QE TeSe* Hh 


ELESLHOG2H 
OLTOZENEEN * 2h 
LSQSOLNGES 
1656922995 
LONE9Q6609 *en 
NTONQEENSD 
O6L062SE0L* 2h 
LOLEOLESSL* 2h 
eH 
eH 
TIGETTOTOT 
TOTTLE9BOT 
“En 
METIOLTONT 
STHEOBNLST*EN 
O6LSTLOLQT 
gostonoote 
OoTsgongne 
L6602966ge "En 
gntegségtt 
S6nSTAIIIE 
oggnTeggét 
“En 
LOOT9SHZES 
"EH 
HITCLETLSL* EH 
oggesngt6g "En 
200ESS9SOT 
J60QTECIET 
LHESHOSAAT 
SHOLESSOLT “Hh 
STT29LICET 


teezoszort *2 

BLNENTI9SE 
TUMELLOTBE 
en 
LTTE0SETHS 
HIEQTELTQS 
O9HOOHOTLI 
eh 
66SS99GNLL* 2h 
Eg 
en 
nn62ST9b96 "2h 
HSIOLGILOT 
eGnQtogttt 
OOLLTEN6ET 
LeogeLnggt 
ORISCELLO? 
MHECTS9ECE 
629SR968LE *En 
LIQQETEOTH’ En 
OTe 
QTLTSTOZeS 
TLT0960995 
Tzésgennse 
ONT * En 
SSETOSOTOT 
* 
ETSQLELETT 
LESLOTROET 


t 


on 
en 


nSOSeingon 
992S02669n eh 
9etoésSenos *en 
Genegsetns 
QTELEQNOSL* 2H 
ch 
THLnLegség* en 
LOTETOLEOT 

OOSTecoeT 
ONeOTSSELT 
En 
*th 
LneogéLgse 
Léggnet "th 
CQSETIETSE 
OGOETLEOQE 
660egt6TTh 
EH 
6LS9TENLTO 
En 
En 
LOCQIETTEL* 
LETQTSCO9g* En 


ELTSE6LQ002 
2h 
SLECHOTINS 
“eh 
egtoLetege 
SQETHATZOE 
QJOSTOSLEZE 
OQTOSHEDNE 
OSOOBERTLE*2h 
TQESEESLOEE 
OSLOLTECEN 
E0egennges 
O2TQLLELOS 
esetnegnsg 
LLAQT9NEOL* 
LISTS9TISL* 
CH 
BOLTHQQONE 
"EH 
“En 
OLTTSLEGET*EN 
ELQOREGINT 
LZO6TONGST 
9ENENIQOLT "En 
HLONOLE 
ZLOTETOBET 
neotnégnte* fn 
6LELZOzEG2 
TeeetgtL_te 
SOMETSEEME “En 
OOTLEO9TLE* EH 
£tgggoncon* 
60ST62S2TS 
netegéésss 
6gTOOTHEOS 
nS6EGNTSS9 
COLETESTIL* EH 


ESTSEQESET “2h 
TSEEOELS02 
en 
2h 
O92 * 
O9Z6S00ZEE* 2h 
LLSELAQSSE 
QLSSTOSTQE "2h 
* 2h 
9L29RSES0S 
EGSECESENS 
OL9ESO929S *en 
en 
2H 
GLIENZLECL* eH 
QTCTETOLLL* 
66TQ2n6R6Q "en 
Q9LEOSELIE* 2h 
TeTeRgeTHOT 
En 
HENSSCTOET “EH 
ELSQOTSTST*En 
ELOQEEOEST 
CTEQEQESLT 
ZEQOSTEERT 
En 
EH 
6060029612 
EL gegnneot 
€L90S92L2E “En 
eggeeesstn en 
en 
NTO9SLIBES *En 


“en 
OS9LNEGEET 
MOTLEQTLES* en 
6TQ6ELATLE* 2H 
ONONSSETTE “en 
TLELEQeNEE 
EQGOLTEQSE 
290TELTNGE 
2h 
NERIELET 
OSHLEZOE0S 
T6995 
QLECEENLLO "2h 
opetgeegeL*en 
en 
en 
Q9LOTIENOT 
IT’ 
EnLONSLTCT 
ZELNLETTET 
T29CETHTNT 
TO6SELNNIT 
ETSONLALLT 
ZLOTSS609e 
6LESTQTege* En 
En 
EnLLOCEOEE 
ZLESTIOLSE 
GOLLLT9ETH* EN 
LELTETANSH 
OO62LLTLEN* EN 


2h 
LOTHOMELET 
ENEESHLONT 
EQnéneggst 
ZESCECELT* 2m 
A906LLOTET 
OLLMELEETS* 2m 
SLQOLETS2*en 
en 
en 
OTHEELLQOE 
* en 
QEMELTHLEN 
6EETESEEOH 
STSSESLEOS 
EQ0EBEQONS 
TIESHTQORS 


LEQGCOTEEET* EN 
LREGLSQOST* En 
LLLLOTLE9T* EN 
62ECQ9SSLT 
GZ0Q00SERT* 
LHTOTTONO? 
£9L0S00T ee 
TOTTEELOES 
En 
QEQEETETOE “En 
“En 
LEEEENOEQE “En 
QLNELEENTH 
TOQROOSE Th 


TSOLENOENT 
“2h 
“2h 
LONLSGLOQT* 24 
OnQO9TLOG2* 2h 
L29ETOLOOE * 2h 
NSLSZETEZE 
TREGEOE SHE * 
LLOMNEQTOLE 
2m 
SogegTagsn 
JenonS206n 
ent 
EQLELOOLO9* 2h 
ETL9OZEOOL* 2h 
O2OENZEOTQ* 2h 
EH 


En 
ETEMEETENE 
TSESEOSTLE 
956500620 
CQETTETIEN 
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tetet 
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OLEQLILT2S* 
QLH6SOTSSS*TH 
BBRS6ELOES * TH 
"TH 
Ih 
QEELZQOESL* TH 
TH 
2SSTLLELT6* TH 
QS2Q66EQLE* TH 
2QOLOOEHOT 
2h 
L26EOTS9ET 
LOOTHLENET 
QLSTHEONAT 
SSOQHELEST 
LSLTESTHOT 
2QLETSESLT "2h 
6QSENLTLOT 
LES22S666T 
HEESENGETS 
LQO6SEE 
STOTEQ6O92 * 2h 
QEBOETTELE “2h 
SO29TLSQ62 "en 
6LEETONETE* 
Le6etégtnt 
QS96NSLTQH 
HTZOONEITS 
LO6SeQnSESS 
20SLZE99E9 
NOLQCEOEEL* 
99966n698L* eh 
O2EOLEOSLE* 
HZOOSQLHOT 
En 


egegetoLnt* th 
Th 
Th 
O9ELEHOLTH “TH 
Th 
Th 
922S6Q6T0S * Th 
SBBISETHES “TH 
9QLLESS99S TH 
* Th 
Th 
T6299L2999 
ELENSQEOEL* TH 
GeLSQSLELL* Th 
Th 
Th 
Th 
6SQ2LTSOOT 
SELLEGBTLOT 
en 
* eh 
ZELESITOET 
6ESTOTEGHT “en 
O6HTTIEGST 
69ETOLTEQT 
ch 
2EzELOEQON * eh 
€90SEESOee 
ENLE69TeS2 
THISSHLE9E 
Lozontggge* 
* eh 
6ETLEEQOEE 
Oz90EZONSE 
6SQTESTOLE* 2H 
LIO6LOTION* 2h 
E90S9TeSEH 
BLELOGHIDN 
22T9L000S * eh 
OOTLES29€S 
QnEQIETSLS 
SOTEZOES99 "eh 
TSRSOLEOTL* 
CETLESHEGL 
OTOEELE6BTQ* eh 
NZT6LS6099 * eh 


TH 
LESTOSHSSE* TH 
96S90SETLE* TH 
SeLSo6egge* Th 
O9TeTHEQOE 
2TQOELSSeE * Th 
* TH 
* TH 
“TH 
GeSeSETITH 
6629S99eHh Th 
LOLENLEOLY* Th 
6LOS9QTTOS * TH 
LOEQIENEES* TH 
69S001919S° Th 
SHGET TH 
Th 
2n206L9S99° Th 
QREOGTLOEL* TH 
Th 
TH 
S9LECESOOT 
NESQEBELOT 
Coenggeoet 
Q60QE6TELT 
60NTESLES 

TEBSES ‘ 
LELLHLOLET* en 
LeEQEEOLO?* eh 
6062egLET Ce * eh 
en 
QeLLogetse*eh 
OOLL9SO062 
LO666ETSOTE 
SS69S6neEE 
6TEQTETISE 
LogSeestgt “en 
ey 
* 
SST2TLEEOS 
QLOOOTSONS * 
* 2h 


LS29S6LELT* Th 
NTISELOLSET “Th 
O£TSS2QLO? "TH 
96SSt990¢2 "Th 
En€Eneente 
Onssnéggne* Th 
L6LLOE 6092 * TH 
SEELETELTE*Th 
TH 
TH 
6TTSQTHTQE "TH 
TLESOLZESH* TH 
TTEZOLLESN* Th 
ectntecoes 
QEIBTSQESS 
TESGLOLERS * TH 
S96L62T699* TH 
TEESTHECTL* TR 
ELOZELLESL* TH 
* Th 
Th 
L2OE0L90C6* Th 
TH 
TLOTESLHOT 
2h 
LEQGIEELET 
Zenen6tSnt “cn 
EHTQZ6OSST 
T2TREOLSOT 
LSLOOGOLLT 
“en 
G6SL69E 
eh 
eTESLENTES* en 
eh 
"eh 
6TSQETLEQS* CH 
THTEEGLEOE 
OQOOTTRELE 
EESLTTZOOR* 


TLEQZ660ET* Ih 
ESTQOSEBET* TH 
OS6ESSELHT 
TH 
OTENSIQSIT “TH 

Coegggt * th 
LO6QT2EGET* TH 
LEgogogtee* th 
TH 
Zegnggeese Th 
TH 
Th 
QETESHLZOE 
SEEQELTCHE 
*TH 
LHTOTEOLGE’ TH 
HLEOOELTTH* Th 
TH 
Th 
THELLTEG!S 
669995295 "Th 
6TLECSTOES *TH 
"TH 
* Th 
Th 
OOTHHeneLL* Th 
S6SQTOTELE* Th 
SESSIESIES TH 
N@NTOLO666* 
QTLSOO99OT * 
QnO699LETT 
QOS TeT 
O9OTS99ECT 
CELLEQHBET 
* ch 
£eg_gzogst *eh 
QERESLEQIT 
QHTEOTSORT 
TIOLTE66ZET 
Ln6S1Qt90e "eh 
eh 
2h 


GETLHLOLOT 
LT@HLOZEOT* TH 
CLQETHLECT 
Th 
SOSHTROBET "Th 
egcognesst* 
HLIT6EGHOT 
OLHEQLENLT TH 
TEOSEBISET 
HSQOTETLET 
* Th 
"Th 
Th 
th 
O9LEETS99E° Th 
H2LELSOTOE 
OTSTLSOOZE * TH 
QETOLTEONE* TH 
SQESSZETIE 
JOTILLENGE “TH 
6€2TELSEON* TH 
Th 
Th 
QOT2LSLE6H “TH 
TH 
Th 
£00S6L296S Th 
"TH 
QLTEQLETEL* TH 
TESLE96E9L “TH 
69L66996T9" Th 


Th 


TH 
Onegsn6n66 
QLOTOLT9OT 
OLHOQTELTT* eH 
* 
CETENLTECT 
TEQETQELAT 
TTSQ6LNLST 
"en 
*en 
QSELSLSOe* eh 


£ZOLLOTONG “On 
LT6£39E969" On 
ZEENESLENG 
O69n9L990T * TH 
Th 
NOBZEEQETT 
QLTQEETLET 
TH 
QOLQTZOENT* Th 
66SSZELTST 
HISTOOOTOT* Th 
Th 
Th 
TOSSTESZET 
* Th 
60nE960LTS* Th 
Longégsote 
TH 
* TH 
Th 
n6QOZOONES * Th 
E26STSSeTE* Th 
Th 
OST6ELESLE “Th 
JeQEGn666E * TH 
"TH 
LILQQZESSH* Th 
TH 
Tegseetets* th 
* Th 
Th 
TH 
TH 
TSOQOTROSL* TH 
L6L220S009° Tt 
S6T9O99ESE* TH 
Th 
TECQITHTLE* TH 
9enSES9EOT 
LSLS9EQOTT* en 
26QTOTTQTT 
OSES9TT9ZT 
ETEQOSRENT 
6ESLONLEST 
TESQETENQT 
6SE9RS9SLT 


QLOETHGSEL* On 
On 
6STLOBNE TE’ On 
On 
TSS0QQSZOT*T 

Th 
HEEQoResTI 
th 
“TH 
LOSQLERLET* TH 
QOETOSGSHT 
JeotgnLnst “Th 
gogTtetergt Th 
eLooggenLt 
* Th 
TIESHENIET 
* TH 
GeQegestee* th 
Th 
QTLO6000S2* Th 
Th 
Th 
nN29SenTOOE *TH 
ESSSOTIETE* TH 
TECCEELEEE * TH 
"TH 
2LOSLEGEQE 
Th 
QLLEESLHEN 
Th 
Th 
0S99£ 2996S Th 
QS2LLS0S6S° TH 
"TH 
TH 
OLQNIOLLOL* TH 
"TH 
QLESHLOELO* TH 
Th 
Tt 
T6ETHSOQOT “2h 
6TSHLOCETT 
QTZQE9QOET 
QLNSELOEZT * 2H 
nNE6BLLELET* 2h 
ZOLNGOELHT 


9689S TT0S9 "On 
626062899 °On 
ET6HSTERQEL 
LSOGOROTLL* ON 
TOSEZSO99* OF 
“On 
£EORSSESLE* On 
TH 
LOQOTTSSTT 
TEZOOORTET 
"TH 
“Th 
COMLEESONT 
TH 
“Th 
926Q9eTSLT 
L9OSLO6SET* TH 
* Th 
gsggnegeze "Th 
Th 
gen0ggetse* Th 
£ETQOTOLIS "TH 
2QQLZELEQe*TH 
ZOLQO9LOZE * TH 
HENLOOTTNE 
TH 
TTEQLLOOTH* TH 
QSTENSOLEN® TH 
Th 
Th 
TTLLO6ELES* Th 
LOTOBOLT9S* TH 
* TH 
TH 
QLLEEE HELO TH 
Sningtenel* th 
OOQLOTTZLL* TH 
TH 
TOO9T9TGLA’ TH 
TH 
nES6L29666° Tr 
* 
QTZegggett 
SSTSCO9TeT 


"On 
"On 
NOLQELTZEL* On 
LSESESOOTE*On 
QLSN9LTE06°On 
On 
OONSSOTZOT* TH 
Th 
EntS€ggtet 
TH 
LETLOSQIET* TH 
TH 
QSTSQSHEST “TH 
"Th 
Th 
TH 
Th 
Th 
nLn€o6egne’ Th 
TH 
HLHL660TR2* Th 
"TH 
OSOBTSLEE 
TEORTHEESE 
TH 
* TH 
9951992955 "TH 
Th 
ESQSOSETEO* TH 
Th 
22ESLHTLTL* Th 
TH 
OETTONSE9R* TH 
On6NeTILZE * TH 
QTOSSLLEQ6* TH 
QLOELEQSOT 
QLEZTOLETT 
25LE96E02T 2h 


ZQEESNTOED’ On 
On 
LOT2E2E9OL* OH 
EQTTHLELHL* On 
On 
TOQOSSTEEQ*On 
QtL6620698 "On 
ETTLITOSHE* On 
On 
ZOLZELTZIT* TH 
Th 
CEQSEGLEET* TH 
Seggeégtat 
QSLLONSOST* TH 
TEROZHLEST “TH 
HONBEESEOT 
LO96SSEELT “Th 
nSEOTSOTET 
* Th 
* Th 
Th 
OnEE9EE RSE 
TOO60SS4LE* TH 
Th 
Th 
EnTSHT662E TH 
6nO9SLE0SE *TH 
CENLLECELE* TH 
SSLHE6OLEE * Th 
£00999S2eh "Th 
ORELESLEth 
TH 
QS2TOQLE0S TH 
6TT9SLEQLS 
QnETOROITS* TH 
9050264959 ° Th 
Th 
SHTQOTESEL* Th 
Th 
* Th 
LEegHTOEOT 
OLLELS660T* eH 
eh 


| 
| 
c 
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TABLES MODIFIED BESSEL FUNCTIONS 


MPSEY AND BENSON: 


O2ETOLELE2 * On 
nEnSO6ONOE *On 
ZORLESTTCE 
Z2TSHLTZEEE 
Tengentgse 
QO9TSESELE 
6Q6LEH000n * On 
On 
LLOTeSE66n 
*On 
OLQ6LES06S * On 
On 
* On 
“On 
69960266EL *On 
TESQOTTEGL* On 
On 
LS9SL61626°On 
StrhS666E 96 
6TLSQTeHOT* TH 
T6QHLECOTT* Th 
SLOOOJESTT 
TH 
LESGLOZEET*TH 
LELOEQHOST* TH 
6EQESTESOT TH 
O0Ze99ESLT 
S20ESS99QT * TH 
QETTSZ2OR6T "TH 
T66OTETOT2 
th 
OOTSERIIES* Th 
TH 
e62026L99e "Th 
69THOCE Lge * TH 
* TH 
QS2QELIETE 
* TH 
* TH 
TH 
HEZOTTQLON* Th 
Th 
TH 
HTEEQEHOEH 


Strh2EQn66T ‘On 
esg6egnote 
6enteonnte on 
“On 
192 "On 
TES2ESESLE “On 
NQZLEQET6Y 
L9SSENLLOE “On 
O6HO9LOSZE 
getezentnt 
Tensgentgt ‘on 
On 
goszoztger "on 
‘on 
ZTTLSHESES *On 
TO62LE299S "On 
9260766965 
6LSEEELEEQ “On 
6£TS2QSOL9 “On 
LSTLELIEOL 
622nOSTISL* ON 
EE6TOS9TSEL* OR 
LOLOOBET * On 
0R9069N TER "On 
HEL 
9EQTZOOOOT * TH 
SSENEEESOT TH 
6TQESHECTT 
Th 
2TIQTLOTAT 
TH 
TH 
EENTILEQIT “TH 
6ST9NETOST 
TH 
OLSSLEELZ2* TH 
Q6ETOENTHE 
2TENTLEQSS* TH 
TH 
QSn6TEZEQge "TH 
SOTTEGeLOE * TH 


LEOOTTHENT 
OZLEQESTST “On 
ENLOLTSEST*On 
£S996S2g9T "On 
EENSLOSLLT 
Of909B2LET "On 
QENELEILET “On 
2679S 902 * OF 
n6TOSETOZE* On 
6SqLSLE2ES* On 
QUTLLZESH2*On 
S60THEOES2 * OF 
gesgonstLe*on 
ENLEESTSOE 
QTLT96SONE “On 
ELT QCOEESE * On 
2LLEESEORE 
* On 
“On 
OTTLE6OSLH 
ESETEONTES 
THESESLHES “On 
$96SL62629 
261696599 
EEHELITONL 
NTSS0666QL "On 
"On 
QEE99ERSER "On 
ELETELEQEE* On 
On 
QTOHETESOT* TH 
TH 
etonegegtt* th 
6ENESLESCT 
ESTLECECET* Ih 
LOSETHEONT “TH 
“Th 
QLETIESQST*TH 
TH 
600LLSEEQT 
69LESLE002 * TH 
6260THELTS 
0S6£S0692e" Th 


SOLESLOLOT 
6L9SLSQETT 
"On 
QLTLOTHS2T 
gégentectt 
QSOBRENSET “ON 
NOBTESOLNT "On 
GOTETOTSST 
6LSL9O9ZLT 
QQELETTZQT "On 
SSSHOQTZET 
egsgtegeoe*on 
9196996622 "On 
SQTSTESQES On 
"On 
19g0S70992 
* On 
On 
LONETOSETE 
OLZ2EOETEE 
TELENEOOSE 
G9SEOOOOLE * On 
LIQEOBOTEE 
LTIQORBENETH “OR 
S9C6TETLEN "On 
On 
2ESTHSOLTS “On 
* On 
TELEQGOLELS * On 
266 “On 
6120616599 °On 
EGLZTLERIL* On 
ST6EZQTHTQ “On 
"On 
LEZTELEETE* On 
66206£9L96° On 
6LHLLZSZOT* Th 
Q96NTS9QOT* TH 
SLOLQSTSTT TH 
nnssnLocet 
th 
QTHOTHELET* Th 
Qeonesssnt* th 
COSTIEENST “TH 
LELEZELEQT* TH 


TETTLESOEQ*6E 
6029TSENLE*6E 
* 6E 
ZEQEOSETLE*6E 
ELOSS9EZOT 
LLOZBRRBLOT 
tHLeZeLETt 
On 
“On 
SLYNELELET* OF 
on 
ES6ES9E9ST 
TEBLLSENIT * OF 
*on 
TONSOS9EGT “On 
On 
On 
eesegnngge 
“on 
gegogsenst 
ZETINENLTH 
“On 
*On 
“On 
TIQLeHEzeS * On 
0206929255 "On 
‘On 
On 
626LE66S9 On 
6QLLLIEEES *O 

QOEOLEONEL* On 
NSQTEEELLL* On 
SOSTHSLEZE "On 
On 
CEMELTLEOT* TH 
"TH 
Ih 
HISQSESEST 


TT609LH699 6E 
6E 
THLgeeg_ge 
66TOOZOLOT *On 
LQOTOTG2TT 
NOLZLEEQTT 
“on 
60ZTONHEET 
T2SLOLOLAT “On 
TTQSLETSST 
9SS9999E9T 
OBRISEIELT “On 
“On 
* OH 
ZOSELLENT2 "On 
OLT9S22922 "On 
ToesoLegte *on 
"On 
Z9ETIOE LEZ * On 
TZONQTQOSE 
ESenSTQOLE 
LESZS0026E *On 
“on 
Q69Q00SE9" "On 
£T9G9S206H 
ZTHOOEIBTS “On 
"On 
"ON 
SQTegnn0S9 
"on 
QTTLE9062L* OF 
LyTLLEO?LL* On 
9£6n999LTR* On 
69E6TZT99Q "CH 


On 


96LNEQEOIS *6E 
*6E 
6SHLNOLT29*6E 
*6E 
ZQLEETESEL*6E 
TOETELSNOL*6E 
TO6RELESOR 
6£ 
*6E 
6159869266 
29TLeegnOT “On 
S60SOQZOTT "On 
“On 
BLELESLENT “On 
TLQ6OHOTST 
EOOLHLEEST “ON 
LENEIBLEIT 
OLNNQOTELT 
2nOLEQELET OF 
OTSTT6EQET 
Loget2en60e* on 
SEELOSHELS 
On 
69STT2OSLe*On 
92920L690E * On 
6ZO00SLERZE 
Eesogntegt *on 
* On 
QOLQORENON 
gotoetegen “on 
LnSQLeLg90s ‘on 
EOOTOSTIES “On 
LLHELEZLIS * On 
SSTOSLEO09 "On 
"on 
On 
TOTSHSZETL* On 
gesencensl*on 
QENOIBLEEL* On 


QLHOLLOEIS * 6E 
0$90SS229g 
nN6OBLEQOOT * On 
L2O6SOTCTT 
QEEQELTINT “OR 
QSSRESITLT 
* OF 
SmnL6ZQOS2 * On 
9692RSREL2 
gntgS6eeTh On 
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"9h 
6LEZETOTZe 
EES 
"Oh 
"9" 
6QT62TeSLE "On 
E6TSERSLOE 
GenTQotset 
T9BEQITHIE 
0069n9nSBE 
"on 
HESTITELSH 
9S002625gh 
S9HOOeHSHS 
22602 TSGLS 
eH THEORETI" 
"SN 
TENLSE9TE9 
629 
TSONSEITEB "IN 
LENOLESLEG 
TOQEENTIOT 
6Z06900LTT* Lh 
EHOLTSEO?T "Lh 
QnQLTTeget 
“Lh 

SOEGOLSAT*LH 
OSOeHTHSST 
HTOOSERSIT 
S6Se9COLLT "Lh 
* Lh 
O96E9LETO?* Lh 
HETELLESTS "Lh 
QnTEgtgote* Ly 
ene ton6gne* Ly 
genoggenge* Lh 
Lh 
LETOQETEOE 
Tetetegnet * Ly 
2QO6LSERNE * Lh 


9S0660S20T "94 
TESGELELOT 
LEGENSLETT 
LROSTLEBTT 
“On 
QLEQHOZEET "On 
SLOELINONT "On 
L2QeSTE9ST 
EZHOLSTHLT “ON 
O6SEENQERT 
QOZTISTHET “ON 
L96900TS02 
THQOTEQOLE 
66Ln6e999e 
ORSERILEOE 
MHUT6LT 
EgegnETone on 
QEISTIQTGE 
GOLLETLHON 
getetetéen "on 
29LSSHTSSH “On 
"9H 
LSOLEET HHS 
QZOHOERLLS 
esgeLgnes9 "on 
ENE TOSLLEL* 
OSTTE THERE 
QORSLTELNE 
"Lh 
E9QLEESLOT 
ELQORELATI 
HLETQ29OET* Lh 
ZTOSOLNEET 
S€gesLégnt Lh 
EQSELETEST 
SS2S96TOLT* LH 


9106995195 
Sh 
N6SQLOSLE9* Sh 
9QTEOZ0099 St 
OSHLSIENES 
Sh 
Sh 
06L0T9$608* St 
TSOT9SS966 "Sh 
9SOTEOSOT 
L9O6Z2LOTT 
Oe9HTSLOTT 
LOMAQETECT 
QeegOLET 
*9n 
OOREECETOT 
ZTOTLIROLT 
9QTZESOORT 
eH TSSOLO6T 
LOZ 
QOEGLSLETe 
"9H 
Leesg6Sege "On 
6060RLE6ES 
QTE9SHEGLE "On 
QTOENHETON 
QSeggTtesn 
TIE 
SQSLEQHOTS 
2TOTHZEILS 
SET9999ETS "SH 
LTESQOTSE9"9N 
9690£0968 


LINTZO6ESE 
TLEIENSTLE "Sh 
LEGETSTOBE * SH 
Sh 
ETHOESHOEN” SH 
HEOQLOESLA 
QTLE909664 * Sty 
* Sh 
ETEGTSESSS* Sh 
*Sh 
SH 
Sh 
ELSELQSCTL* Sh 
St 
EQOTLTSeER* Sh 
ETEZSESLLE "SH 
HESELLSHZE* Sh 
TEOR6ONGOT 
HEZQELEHTI 
TOTELONLET 
696SQTSHET 
6T69Q9OZNT 
LLOEEQOOST 
QLEBLESRST “On 
“9h 
LTSO9OEQET 
QOTEES 
Olggét 
"On 
TIESSLEEQE "OH 
OLE TE629962 
TLIESESHTE 
6TEOEGLEEE 
Teeggtenst 
ZTELESOGLE 
GSLLETERSH 
8S0676605 


Sh 
STTONZELHS* SH 
LLE6600092 * Sh 
6SE6LETIQE* SH 
6QTESSZ00E Sh 
SEOQOETSTE 
QSEE90QOEE * Sh 
NOGETELLHE 
Sh 
O69Q9ETERE 
OLLOLOSCOn * Sh 
* Sh 
2TRESTELON* SH 
LEBTEIE TEN’ Sh 
SQOEOLLOTS*Sh 
LTMBHOIENS 
HSMQESETLS* SH 
0992626109 
Sh 
SH 
LLATT9SZOL* St 
6TZETIOONL’ SH 
2STLOSLELL* SH 
Be26e92998 "Sh 
TSS9SELETE* SH 
SH 
nSOOL2E9TOT 
OOTT6ESLOT 
OMHESOSETT 
6c0e90e9cT 
ELOGETELET ON 
TOLLHETILT 
6LECSRELET 
QEEQLEQQOS * On 
EgnLottnte 
OeetSLngne* on 
CESSHTONTE “Oh 


QIESSESLLT 
LenzLSoggt 
HETEQSENET “SH 
O€GeLTEHO? * Sty 
Sh 
6LLHOESHES* Sh 
SSQS6ENSES 
ZEnSHZEQne 
Sr 
SOST9SLTLE*Sh 
OOLANTE6ES 
GoesLtente 
* SH 
O£20S96E9E *SH 
Sh 
699LeLtzen’ Sh 
Sh 
SH 
60LLE60T6N* SH 
T9T9LEQITS*SH 
* Sh 
QS06NSTELS * SH 
ogttozt2o9 
"SH 
SONTLEQEOL* Sh 
SH 
8S8L0L8898"° St 
Sh 
Sh 
ZTEQLHOZOT 
2Te20SLETT 
QIHTSHTOZT "On 
OQTESEE9ZT "ON 
S26SLATHET ON 
ELOGOTETHT 
SOLESS66NT 
"On 
STLOGeQLOT 
LS2S629LLT°9N 
620n29088T 
OLHSTLTEET 


QLELIAIOET 
L29EQLLIET "SH 
69eQELTENT’ Sh 
HEESSEQENT 
626%6S69ST* Sh 
QT9ZOBTELT* Sty 
HEL EQEEORT "Sh 
TOSESTIGET 
Sh 
SSTOgOE ge? 
sn 
Sr 
Sh 
LetLeLsgle*sn 
9906962062 * St 
2QELTILHOE * Sh 
£6202TO0eE * Sh 
HILLEGO9EE 
SQEEQSOLSE “Sh 
25622S60LE 
*St 
OnS6TLLE0N * Sh 
St 
LEELSES9LH* Sh 
NOQTSEETOS 
£962LL26SS "Sh 
9009 *Sh 
St 
St 
QOEHTETERS 
Sh 
Sh 
Sh 
TLT60L0068 * Sh 
2ELQLLTEEE 
2Q66092T 66 * St 
LEEQTSOHOT 
LLESQTSOTT 
TLETERLOTT "On 
2000nSELST 
ZOQTGLEOET 
ON 
O6LLELLSHT 
SOLEQTeHST "On 


TQZ9O LHGOT 
T2TSOeTHET “SH 
ET2LEQQECT 
THENENESET 
ONGISTESHT 
OTRQOSOOEHT * St 
ETLHSEO9ST 
LTOLS2HEQT "Sh 
SELZTIELQT’ Sh 
QeEOLOOLET* Sh 
Sh 
QSTOTESOTS*Sh 
sn 
* St 
Sh 
SOLE0STSLE*Sh 
Sh 
LEQTELCEOE "Sh 
QQTeLEngtt 
* Shr 
LOQSELETSE 
LygLenoget * Sn 
OTOLETELON* SH 
Sh 
ZOZE9LOTSH* SH 
Sn 
nSTEEONSES *Sh 
ZESOTLOESS 
*Sh 
OOSEQLZETO "Sh 
2120620979 
£0ggL69089 
SH 
6S0T9Q29SL*Sh 
TLOLLLALEL* SH 
O0gezETT Hg 
QONTTTNLER 
Sh 
EELELHS Sh 
BLLEGLEHOT “On 
OOLQEHZOTT "9h 
“9h 
606TEQOECT 


GROTHSLAIE 

60RnLZESOT "Sh 
Sh 
HT6LESeTeT Sh 
QeLLEQQ92T 
Sh 
2LO6NTO6ET "Sth 
6LUQTHSSHT 
2EETTONCST 
TIETIT9EST 
TLLTCQTSLT*Sh 
ZOTLSHSERT 
Z6SOS9EZET 
ggegnngtoe*St 
Begg 
"Sh 
Sh 
TOSHELLOTE 
*Sh 
*Sh 
* Sh 
* Sh 
* Sh 
St 
* SH 
S200ELHeTS *Sh 
OOLOZEHEES * Sn 
OLTSTLOQES 
SH 
nS29E6H669 
NELELSELEL* Sh 
OCEHTLALLL* Sh 
Z0SLOQE6TR 
"Sn 
SH 
Sh 
€Q9O6SLTOT "9H 
QLLELEHLOT “On 


2TLEELOLES 
EEETEQNSLE 
Qn6ELEETOT’ 
OB9ERIS9OT 
ZEQIOLSTIT’SH 
* Sh 
TESOSELLET "SH 
ELECSSLEET "SH 
SH 

LHEOLTSEST*SH 
HLNQ6LLOOT 
Sh 
QTOTTTENGT "St 
Sh 
GOSeT9TEO? *Sh 
"SH 
THIELLEE 22° Sh 
eTORL Gente *Sh 
099906LSHe *Sh 
gte2ottnge 
nZECggegée *Sh 
TTZEEEZETE 
Sh 
ZELSSZOSHE 
TENCTLIEQE “Sh 
TQOO6EITRE * St 
QEEQLEZTON 
Sh 
SH 
HLELLGOTEH 
*Sh 
QTLSEGEENS 
ELSOLOLELS*Sh 
TLEELOZEO9 * Sh 
"Sh 
OLBHETSE99 * St 
OT9ELS9SOL*Sh 
Sh 
Sh 
QENELTELER 
neLegtge.e* sn 
69LETTITZE*Sh 
£566S9920T "9h 


° 
a 
w 
c 
a 
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FUNCTIONS 


BESSE 


TABLES MODIFIE 


NSON: 


AND 


DEMPSE 


Sh 
OSEOESSLZE 
* St 
St 
BTESESENEE 
6LEZ2OTETH* SH 
Sy 
Sh 
SH 
SEROTILOEH 
Sty 
* Sh 
QTLOORHSLS* Sh 
S290R6LE09 
St 
LET9TEOS99" SH 
"Sh 
O£96S662EL° Sh 
262LQnL69L* 
Sn 
Sh 
TOLOS2H26Q Sth 
Sh 
Sh 
TEQETHOEOT 
9SLOZQEQOT 
QOTZQOSEET "On 
OLOnE2SONT 
6ETIOSELNT 
LO6NIELSST "ON 
ELE6ESTSQT “On 
O6TLELETET 
HLELZEEETS 
"On 
"On 
nQE9TOSOSe 
S96S999n62 "9" 
QETZLECTIE 


S92926SS9T "SH 
LOGEETTELT* SH 
266nSTOTOT 
€90Sg626gT 
QT96E0TLOZ* St 
Sh 
EogeztLgee* Sn 
Sh 
8992200092" Sh 
SH 
SETOELSQ6e° SH 
EQ99TOLETE "Sh 
SQ99LLOLZE* SH 
QSOZTOEENE * St 
NQOLEOSLLE * Sty 
TZ09SSnS6E 
Sh 
6SEEnLENEN 
SH 
St 
* Sy 
ZETLSEE9ES * SH 
T2ELHLTESS 
Q6EEO9LORS 
"Sh 
HONTEGHOND 
St 
T9OSEQLOGL* Sh 
Sh 
Sh 
6Q0£099L06* St 
6TTOELLHSE 
QLELESHOOT "On 
ST9QTTLSOT "9H 
SO9TgeTIT 
SHOQTHILTT On 
TLLLISEEST 
LETHETEECT 
HENEQEETST 
6ST69S6Q9T 


TTQO6RZOE6 * 
C6LQHESTLE try 
G2ESQTHIOT 
OEZOETESOT "St 
ZeQenegort 
LHOOTOSSTT* St 
NOBLHELOST 
OORBLSTICT “SH 
OLETENQTET 
QLOGSSONAT 
OOZLETIOST "SH 
TULTTEMLST 
OTSSTQZELT* Sh 
2QLEOEZORT 
* Sh 
SOO9SEELET *Sh 
Sty 
SQTHQETITS "St 
SH 
* Sh 
Sh 
Sh 
* Sh 
QCTHSOSETE "St 
* Sh 
TIE 
SLTOLSO6BLE * Sh 
SNTSESOLEE 
TS9COSELTH 
EQSTLEQESH 
Q6T9S9990S Sr 
QTEEENTCES 
TOO6LRERSS * St 
2556005 * SH 
SH 
Sh 
QOLES9OLTL*SH 
QLSONNTHSL* Sh 
Sh 
ELHETEERLO* Sh 


ORELQZTELS * rh 
QOORO66E09 
1609666629 
66LS2T2LS9* 
€TL6GR9SE9 
LTZOCERSTL* 
HZLEESEENG 
nSSTT98988 
NOLETST9ZE* 
60TRSTELIE 
HEIOSHOTOT 
9OSLS9SSOT 
QTHENOLOTT 
6TST2LESTI 
Sh 
SEQQENSSET "Sh 
9eTSOT9OST 
2OLOLSSLST 
ZOLLOSQHOT 
OLEEQOS2LT 
SSTOLE6QRT 
TESEHETLOS 
"Sh 
NOLLOTSENS *Sh 
Sth 
Sh 
* Sh 
QLOSILESTE 
EQHOSENTEE "St 
QOSTSTLLAE "Sh 
69QT TINIE * Sh 
QELEISETON 
96g09TOZen 
SH 
SHSSSLLETS*Sh 


OLQTOE GREE * 
THHECE TEEN 
29Q0RSHTSH 
Q2ESLHIOEN 
HSOTSOSTTIS 
THEMROSEES * 
ggzgzengss* th 
* try 
69TL266L99 
6LN60S9ENL * 
OLQETAIZEL* rh 
6SLE0R2LTO 
HZELLOST6Q* rh 
TES * 
HOHOLZOELE * 
TETOGLOTOT 
SH 
QENTETTOTT’ SH 
HOQLELLEET 
"Sh 
LIELOTESHT "SH 
HQSOSETEST "SH 
6002LSS99T * St 
EQ9SEQT LET "Sh 
6TSHTS200e* Sh 
S99 QLLE02"Sh 
TIELETS "Sh 
Sn 
LEGESTH THE 
OnLasgegle*sn 
2992S9QT62 "Sh 
SETASQTIOE "SH 
Legéeggetet 
HECQRSTLEE 


Hh 
SOLQETOOge* 
th 
SLOOERSEOE * 
HOOSSOITOTE * 
69RETOEECE * 
QQSTESOENE 
ES606SELSE 
LLEQQEESLE * rh 
* rh 
* rh 
QTLLLSH Ten 
QO9LQELEEN 
Leesstogsy 
OSQORRS LL" 
SSSETENETS 
LTEHSHLERS * 
EQOLIEESTO* 
OnST2LEOLO 
SSTELTQ669 
SLLLLSQOEL* th 
995929629 
ESLTQORIEL* 
0255011699 rh 
* rh 
9626190266 * nth 
QTEELE9EOT St 
ELSQOSELTT 
TOOEQESQTT 
LLOSQLEEST * St 
"Sh 
"Sh 
LLEHEQETHT Sh 
ORST9ESENT 
6gtcenssst 
LLELELGEIT* SH 
TETTQOSOLT 
GEELELOLET* Sh 
GESOLNESET 
660NSTESO? *Sh 
OOOTSESS22°Sh 


OTO6E9T TO? * 
ETO6ELZ602 * 
trty 
LSeLEQS922* HH 
6EEGEORSES * Hh 
ensggensne 
QEELOLHSSE 
95965S669e 
* 
* 
QLEEOMLETE * 
* 
STHQELEGEE 
TIQOTESESE 
Q2CEELELGE * rh 
QHI6ETZOON rh 
QTCTOQTLTN 
* 
$S6S2062L * 
* th 
66LE2Q99ES * 
QSOBINEESS * 
S606LE2NGS * th 
$%S206L609 th 
BLESZE * 
6S2Q90LE69 
“Hh 
nOZOTOS9SL* 
th 
OT9ZE0ETO6 
OSTHSOSHOE* 
SQSOLZEZOT* Sh 
TIBIGESZTT "Sh 
ZELLSOLLTT’SH 
*Sh 
Sr 
Sh 
“Sh 
STHB6609NST SH 
TOLOLEQTIT Sh 


HSLLIOLOST * th 
QQTLQZEELT 
* 
HLOESEOSET *trh 
* 
th 
* 
* 
QQOOLE6LES * 
HSHETLOLNS* 
OLELESSELE* 
TEE * try 
* 
QLELZORSTE * 
TE66HELLENE * 
QOTINLTLSE * 
th 
69SLOO96En * 
th 
6E9TLEORLA 
909068986" 
* try 
SSLE9GLENS * 
* 
SQOLTOTTES * th 
LLEGHLEQTS* try 
* rh 
O6SQHLTCOL* 
SHeLLS6S9L* try 
O6SEEQTOOR 
nh 
th 
TOQORZOE TE try 
QLEEGIENSE 
QETSCELLEG * 
255996060T * St; 
Q6EEQOTHTT 
* St 


* rh 
TIESTQSTRT 
QLOSTLILAT 
OZTQLEG2ST * 
rh 
OTHEZONSIT 
O29LHOOSLT * 
* 
22OELSSEET 
HTIQTEQETO?* 
* 
HESTEZOQTS 
Ot 
* 
ESGELEQSHe* rh 
Q9QTOZESSE* trh 
96008S199e * 
STOBTONLLE* 
6666E trh 
HILSLS600E * rh 
EZLOLENETE* 
ELTSQ6S9zE* 
SHRENESHSE 
O9S6EnSE9E * 
HLEHTETSBE 
LHTSS2STOR’ 
THESSQTSSH 
Hh 
15590256 * Hh 
Q6620299TS * 
OSSOZEOBES * 
MONETENZIS * 
SQ082S0L9S 
TSHEOLLETS 
Q96ZEOLEES * 
* 
TLeHTENQeL 
6959088092 
SST2TEQNEL * 
S6HES26L99 
T9OTSZTLOG 
QTTHOIE T66* rh 
S2Q6T9QLOT "SH 


*trh 
2TETELEEET 
EETOSSOQET 
LESLIELENT * rh 
LOQOZELSST 
TS2QHZETOT 
QOGOHTTSLT 
2662SE * 
ZEQERSHERT 
TROCLEOLET * 
* trh 
TeELTQEETS 
* 
* rh 
rh 
rh 
* 
TO9TETLOOE * 
QOEQOTSETE 
* try 
NOOOLLEONE 
NOTOTENTIE * 
HEBLODLOLE 
rh 
TISEQEHEON 
OOLQSLENIN * trh 
TLL22QTSOS 
2QIDNLOLES * 
Q9EGEREENS * 
O9BCELEELS 
SHTELE066S * th 
HLLTTONSES* Hh 
6092TS62S9 
£99eT TATOO" 
HETOLELENL * rh 
QT9LTHEGLL 
6206TSLTTO 
259L9TS988 
9£99269926 


SOSTQESETT 
Q2SLSOSHET 
* 
ZOSGHETHET 
SOLENTHEET * 
“th 
NEBLOSIOST 
SOEQLSEGOT 
QLOCEETOLT 
ETTESETEST 
662L69S06T * hh 
6EHELSEQET * try 
OLQOBLE9O? * 
Lenégngnte* 
6eQtogttee * th 
th 
HE * 
th 
QO6ZOTES9e * Hh 
ELEMQETE LE * 
try 
QOLOLESQOE * 
Legtegntet * 
6SQEQSENEE trh 
* 
CEQLOLTELE 
EQSETTESEE * 
* 
* rh 
LETLSOOR * 
SOLLLHELON 
* rh 
SOSSLHnGOS * 
LSTHOENOES 
6ENTILSESS * th 
* Hh 
* th 
LLTEHES629° 
€5Q0982LS9 
LLLOETEITL* 
S2696ELENL * 
th 
TOLTL9ZLTQ* Ah 
th 
lL6S26Q nth 


a 
4 
Rim 
4 
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a 
an 
ao 
w 
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° 
c 
i 
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FUNCTIONS 


SSE 


DEMPSEY AND BENSON 


QETHNSSLES 
"En 
LETTOORILS 
£4 
EN 
"En 
etnegngetL 
QENLTOQEEL 
Q6E9EOTHER 
tn 
6TELAONZOT* Hh 
an 
HSOTISCOLT “th 
ZELSLLENTT 
O2TeggggIT 
6TRCOLLLET* Hh 
OLSS96SELT 
THELITOLET 
gntgeLegnt 
OOMESZEEST 
O9E9B0GEST 
eLQLOEESIT 
NELEQCOLET 
LZ2OTTSO02* hy 
606 THOL GO! * Hh 
HgTe * 
T60EeL 
eGQSLigtte 
th 
é9gtegtege "nn 
4h 
MOLOTSSEQ? 
L29969590t 
LETELTERTE “Hh 


ELEOQESLLE “En 
O92LOSTETH* En 
E6SOLSONEN*EN 
60S£096S9n 
EQeEEneegn 
SLHS6ES66H 
EQTSHENLTS 
HLLENETOES 
nSOREBNSSS 
gt 
nSegetgong* 
LSSEeTONIO 
eee En 
ETOLIGESEL* EN 
LIQEESCLOQ*ER 
6£69L90026° EH 
oh 
EN 
ognggeLeot* 
Th 
ELIQLOSOTT 
HLTLOBLATT 
6096HETELT 
LESTLytget 
STLOONZEET 
QLOLIEEQET 
THSQOESENT 
QEEQSHTENT 
9e99SBOLOT 
LLESHOSELT 
HETHTTEOQT 
QTEELEELOT 
Hh 


THe" 
zogeoo6ese 
EH 
609606992 "En 
€9L0SH696e" En 
QneOESELOE 
ST9SETOESE “En 
HOTLSTSS9E 
"En 
THOSTSETEE 
9ST9e26S0n 
ORTOELHSEn 
HITQSEOTSH 
HZ6THTELON’ EH 
GgLeestngn’ 
GET2SE9TOS 
nI9Z6QLETS 
TOLLTHOQES *En 
o€2ttzegss 
E2OLLSSQLS* EH 
OSOL9e9Te9 "En 
"tn 
TLyS6OEQTL* EN 
Sogoségnnl 
"tn 
EQONETTTEG* En 
EN 
62E TOTOgeE* 
LTQOTZEZ96" Eh 
S26T “EN 
QnZOOOLEOT * hh 
NS90629LOT * Hh 
QEQLOTLITT 
6296696S UT 
QEQLSELOST 


eggcensert 
tn 
STHeORHLST* 
009n0ege9T 
ETLONGONLT 
Engtecoogt 
LagnSLT9QT 
LQOITSSZET*EN 
negssortete 
2Lgotgnoce 
6290L66SEe° 
OLEZE0STIE* EN 
TSEQENIOLS 
96602966ge 
SEQEOITOOE* En 
6eQQ9TLIZE* En 
OOSLROTEEE*EH 
602neQTLSE 
S6TLETTEGE 
OOSERST 
ETSOQTOELN 
292LLELOEN 
TSOTTZ29QOS 
Lytgngtgns 
22ELELEIIS* EN 
QEETEETLES 
SLOCTESECL* 
ETTZESTLSL*EN 


ZOO6SLSTET*EN 
299QSOO9ET* En 
6OEEEESONT 
enLognesnt 
OLSLLOESST 
* 
ZOTINELTLT “En 
nensogLegt 
@LESTLOOET *En 
Teegetggét 
En 
EN 
TEELEIEASe 
En 
69LOTE99Ee 
LOTSTLTLOE 
LOGETLOTHE 
* 
629S9969LE 
2609869907" 
G9GLZES9EN 
HSeggonesn 
En 
ft 
6TOOBO9EOS 
TETSTOZeS*EH 
6EOLBETTNS 


0865600951 
er; 
LILSLILetQ* en 
209090098 
GRR 
LOQTHEORNE* en 
er; 
En 
LETHHEOGOT *En 
th 
LELLGLELeT 
GegenTIQET 
CTELLTLONT 
getennssst*tn 
LEEGINGOIT* En 
QSIETHEQOT 
En 
QGETLOSONGT 
EH 
THSTESE96T 
tn 
HTCOgoTgot 
LQLIESOBTE 
Toggtezene 
O906TLEORE * 
QITOSLONEE 


en 
OETLOGONEG* et 
69T9NGSOEL* 
GELOLOTHNL* en 
neg. *en 
QGLESLOOLEL* 
en 
Ltoeg099Ng* en 
LECQIENNLE* en 
* en 
OTSQEECELE* 
SLLSESTHOE 
BINLEST9EE 
LELEQZEZOT* En 
BENARSEJOT 
TESSOQSETI 
En 
TONLETHSET 
62QTLNQEET 
QLTQICONHET 
Langeget nt’ tn 
OSONTE TANT 
OBTSSLTIEST 
LEOETONGST* En 
SLT LEeHE9T* EN 
JEGESHeSLT 
LOnenLetet 
nOL9eeSLAT*th 
QOLSnTLooe 
En 
"En 
JonToétgte 
EN 
enO6TQTH9e 
66QCTOSELE* 
6ESZOLTEQE 


EGTEELESEQ 
LOGELEN9S9* 2h 
en 
LEENLEGEES * 2h 
en 
LEELONEgHL 
GOELELLOLL* 
QETLESOQEL* 
Ch 
* 
2h 
LOLEOSELQE* 2h 
* 
TREEHELQOT 
TT’ En 
EN 
QLOTLESHET “En 
LESLOQOEET* EN 
TETLOTREHT 
6TOTROLGNT 
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LSEQLLEZES 
"en 
Legenteteg*en 
2S62EHONEL* 
eh 
g9S92908n8 
ELSEQHSSOT 
LHZ6CO9LTT 
eQnStteeet 
6eLELELTET 
QTLYnT6THT 
En 
LENOT 
Jengn6LLLt 
660NSELTET*En 
6Te9gg690e 
OnnLogeTye* eh 
TSLZ6T9O9S 
2EELTEQTES 
O2LLENTHOE 
SLLESLT9SE* En 
On66QSnSee “En 
2Q9SETELTH 
S6029LETSH* En 
S2960TLEQh* Eh 
nOSEQOQOES *En 
EH 
SOLLTHT66L° 
HH 
QOS9TELTIT 
LELESEQTST Hh 
LSQT29E AT 
BLEESQELST 


6£g9699R9e "en 
LOLNSLOLGe 
QINTEOTEOE 
LLETSQOOEE 
6TSEOBOBLE* 2h 
Letégengon * 
OBLLASESER 
E9QS9ESLON 
69TSLS9OTOS 
OBINITHEES 
OSTS62ORLS * 2h 
*en 
* eh 
en 
QESLLEOOLL* 2h 
9S0Seggleg*en 
* 
6S9NNELLSE* 
6nOQLOOEOT 
TLELAHEOTT 
TLLETONETT’ En 
TEnSeLgget 
Legngtgget 
TELQEOTTOT 
HLESESLELT 
ZTEHESHLET 
eegtgtteoe 
néoegengte* 
602TZOESES * 
6n69STTeZe 
Lgesgesent 
LeHTLHeLLe 
En 
HOZETESELY 
6ETIBEGESS 
OT2SQLLTTO* 
O99QELSTSQ 
900ST00926 * 
OngS6ELOOT 
HTESSH96OT 


629620602 * 
60690L19S2 "2h 
LySgtLenle*en 
OTITOELEG2* 2H 
QOBOBETLEE * 
EQTEELETOE 
EQHSS9ELQE 
EQOSEGESTH “2h 
en 
QZEQQZOETS *2h 
QSQ9SELOSS 
LSS9OT2SEQ* eh 
Lenetgnegg*en 
QEEQLENEEL 
en 
2h 
SSHZLOSTQE 
S9TSSHOSOT 
QEOEGHLETT 
LegnTTSeet* eH 
QTLEZEETET 
QEQOTISCET 
Songeosnze* 
En 
EH 
SOLLTQEEQe "En 
EH 
EQETIOERSE 
LISTLOSHEN 
En 
* En 
6THOQESTED "En 
g0gso9gse9 
SESHTTTSHL 


O9NQTAQOLT* 2h 
TTggtoLegt 
JEQOOEESET 
nee2L0060e * 
LELEGZEEES* eH 
eh 
S926 
6LSQINLATE 
ELSOPLELEE 
ZTELEOSTIE 
S2TLELSLEE 
ONTSSH9STH 
Of “en 
eh 
2SnSQLHETS *en 
OneSggeTss*en 
LOLELOOZES * eh 
TTSELOZERE 
En 
6LIQQIEETT 
6z0€egLeet 
OONTLLEZET* En 
S6929LS2nT 
TTENQeLEST 
LE6LEORSIT 
QO6S669RLT 
HE9LOGOLET 
L9egogngoe* 
TOROLLTSe2 "Eh 
6S2n66ZEne 
"En 
Q9TégnEnge 
HEZBHLSLOE 
ztggnegeee 
T9SELQZOIE 
ZOZEQITOBE * En 
QTETEST 
QSLHLELGES 
* £4 


STTQ2SRE 
9SOBEHLLIT* 
SENSENEELT 
NSQELELTET 
606S66TS0e* en 
OT9TSOSET2* 
O662TTENES 
2h 
eh 
2LSOLSEgRe 
SQTT9OEQOE 
ELQOSOOTEE 
QTE9LLLHSE* et 
GOnTESEORE 
TETOZEGLON* 2h 
HTEELESLEN 
* 
O9BTEQREDS 
TOSHZEEONS 
SETHEEQORS 
E6gesgéteg 
2h 
S9SnS2902L 
ELQOSOLALL* 
2h 
NENELESNIE 
T6ngeegeor 
S6QESOQTIT 
SSn66en0eT En 
SHETTOLEET 
ONeNTOQEET 
62HOQéLOST 
nSnLangegt 
THE 
*En 
QOETOEtHOe * 
go90tggoce* En 
"En 
En 
tn 
OLEn6S9TOE* En 
“En 
ETETONQERE “En 
TOTTTNSHTH 
En 
O6SLELEQES 


69S6HTLEET 
QONEBLLTAT* 
ZEQLQENTST “2h 
eh 
QENTTOTELT 2h 
EEOLOLELET* 2h 
NQOCTELTTC* eh 
2h 
2ESETLSeHe* eh 
S99ELEQ6Se "en 
ELNSEOTETE 
EQ6OLHOeHE * eh 
* 2h 
2h 
SS9EOTHESH 
£QLLOOSSEn * 
6SSLOSTTCS * eh 
QELTERSESS 
QLLELZOTO9*2h 
QLOLELEEEQ* eh 
en 
en 
eh 
OSESSECEEE * eh 
6LOGHIRSTT 
TISHSEQHET 
QSOTESSHET 
TEO*ZZOSHT* Eh 
LEZEETLEQT “EN 
nLnlggoegt * 
QTORSHS9ET 
En 
LETEIELLHE* En 
TQOHESHETE 
TLHQTEELOE “En 
EH 
Teevgstten’ 
OBLTEEELON 


O92ELST9eT 
NO6O2SLHET* 2h 
2G60S96ENT * 
OQOOEEHQEST 
HOENLENHOT 
OSHTLEELET 
QLZO9gOTO?* eH 
Zonseengne* eh 
TLE6SSS9E92 "en 
en 
O2Q96E920E 
ZOONTHENZE 
TLOTT29LHE 
€eSeoLgeLt *eh 
6290TEQ6EE 
"eh 
LELZETEESH 
6£9996SE6n 
6nLegegées 2h 
6LS260069S *en 
Q9ESEBTTIO* Cn 
SnSS69959 
EOETZELSOL* 2h 
cn 
eh 
ESESEELLLE* 
TOS069SHNE 
TIOHELOTOT 
OLLESLHEOT * EH 
EQOOSTELTT 
HTIGLSHOLET 
En 
LOQQETZEST 
EQOCELLTLT° EH 
TELHEGESET “En 
2TZOLSto0? 
QT6eSSSEEe 
OTHnLenese* 
99T0002S62 "En 
En 
* Eh 
OSOTSSSEEN* 
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NCTIONS 


SSEL 


MODIFIED 


TABLES 


AND BENSON: 


CE6EQTHTED 
O622ELTOEL* TH 
Th 
Th 
tn 
QTSZEHELNE 
BNLELEGLOT *2h 
eh 
Gonesnhncet 
ZLOTZOQOET *2h 
EQRROLLEET* Cn 
NBEBESQORT 
E2O60STEET 
TLSQTeE90e * eh 
LETTEENOZS* 

TEHLSSES* eh 
Ennonteége* er 
LOLTLEELOE 
TESEESHECE 
LOgg9SS2SE “eh 
LET9QIELLE 
BNITEZOHON 
BTEINIVEEN 
en 
6LEEELTCES* eh 
gtegLenols*en 
OL99S29TT9 "eh 
QEESSOLEOL* 2h 
GHEQSTISSL* 
2S2996n0TR 
En 
HLn6E9QLOT 
QLLTEEESTT 
ntgetagnet 
LOLLTNONET 
O9EECETHAT 


QELETTZOSH* TH 
Th 
6622654605 *Th 
6TT9N60ENS * TH 
OLEELZEGLS * Th 
TH 
OLTEGTQES9*TH 
TLee969°TH 
ELEZOSSTHL* Th 
EQSeLT66QL* TH 
Th 
ETSQETEIER’ TH 
TH 
NONLZQIQOT * 
SQSENTESTT 
LEEGIEBTET 
LEQOEETOST * eh 
2QESEIETLT 
LSQOQEOEQT 
LEQQOSSSET 
998nS96g0e 
nesnosttee*en 
TEQLELLEES 
QOLELEESS2* 
en 
ELESLEQHEE 
GeTSgtzonn’ en 
60TSELLTLA* Ch 
O£OHGOLSOS *eh 
N9ZOOTSTRS 
"en 
QORELELTLL* 
henLentgeg*en 
9861908688 "en 
STQELE09SE * 2h 
HEQLTGLZOT 
TOQEORBHOTT 


6LHTETHLQe 
OTEO9QESOE * TH 
6nQseesnze* 
TH 
“Th 
6S0S26Q69E * TH 
TH 
gtonozornh 
th 
OTQORTZEEH* Th 
ESLENQeTES* TH 
€9TSg90209° 
ELTSSOTTHO* TH 
zegesogegg 
2ZOQTLENLL* Th 
Th 
TH 
SLOQEQOOOT 
en 
TSQESQEQET 
TEHTENLLAT 
O6ZQOBLLST 
"2h 
OSSgeS00gT 
LZTESONZET* 2h 
OL20€S9S0e 
OSHELOTSES*2h 
6ETCEQHTSE* eh 
er 
GESeSHTQOE 
"2h 
LOCTSLLS9N 
S9gnoené66n "eh 
“Cn 
LOLHOEHLS 
6ELENEQOTL* 


TH 
énSnezéste*th 
Th 
6ELOSGHENS 
9ST0S299S2" 
O299ESLALe*TH 
Th 
EZETESZOTE "TH 
OLTLHLLEZt* TH 
LS669QS0SE*TH 
T6OLELLELE*TH 
* TH 
QEGTCH “TH 
ooesgtoenn 
TEEHQRELLA 
SSLT9LOGOS * TH 
* TH 
TH 
S6TOTETETO "Th 
TH 
T60EE9LSE9* TH 
609SENETHL’ TH 
n6ESOQOOEL* Th 
TLHLLO6LEQ* Th 
TSLSHTSLSE*TH 
LSZO6ZT2OT* en 
LOETHSEQOT *2h 
6TLOOGZOTT* en 
HCETECENTHT 
CLENLACTST 
HLLATSETOT 
ezeHonengt eh 
SQ0SZHE9ET 
HELLOETS22 
EgeLnggone 
QLELLHTOTE 
AMHETESQEE 
L6660eLe9t 
LEETAIIERE 
eh 


6£09SQ00ST 
E6EETCLEST 
* Th 
LINGEQZELT "TH 
QEETEQZOET * TH 
Th 
th 
TH 
Q609QOETHE* TH 
TH 
ELENOLOELE* Th 
LS9EHLEQOE *TH 
LOgOSTELZE 
EQELENSENE “TH 
TH 
QEETOESHEE* TH 
Th 
H2TRELNSOS * Th 
Th 
TH 
Th 
Th 
th 
LEENELLIGL* TH 
"TH 
BINGE Th 
LLECEZEESE* TH 
NIQBTOLTOT 
LOQOELSQOT * 
EMELLLESTI 
TOSQEEQEST 
OFBLSEOTHTI 
Eggeangost 
"2h 
Z2STOOOZLT *2h 
Lognengest 
£QOSESS9ET 
OLNgeenLse*en 
LON9ENSSLE* eh 
TETSEZOS62* 2h 
LSOBESHAEE 


Of GONSSSTT 
th 
O26HLOGECT * Th 
nEZOSLLLET* TH 
2EETIETONI TH 
ZONTEEOSST 
SLIQOLSHLT 
“TH 
GETZTOL9ET 
TH 
Th 
QSLesgsste* th 
6QTTENZOSe* Th 
"TH 
TH 
TETOTTEOOE *TH 
ELENENZETE “TH 
Th 
OO6TELEO9E *TH 
QTESeS6ERE 
STegegngon* Th 
HESTTegney 
666g0nSe9" 
QELIGLTHES * TH 
*TH 
TH 
Th 
LLTOC606TL* TH 
TH 
th 
TOO96OHTLE* TH 
TH 
TLIAQTST66° Th 
6nZ6LELSOT 
QLEONTECTT* en 
“en 
THnSg69get “en 
HLTQ6TAINT “eh 
QLETTOCELT 
QSNNEESTET 
* 2h 
OLTETQE0S2 


TEZLES9GZE 
LTEMINONGE * On 
TH 
6LEOHHSOTT "Th 
TOSLTETLIT* Th 
TH 
6LEONGLEET* Th 
6ntgegegnt “Th 
TEGLSOELST 
“TH 
ELEQSOTLLT 
QSTONGORRT* TH 
TH 
TH 
TH 
Th 
th 
TeegstooLe*th 
* TH 
TH 
HOSHEOEHHE 
OLSS62g99E 
ELSQTEZOEE * TH 
Th 
Th 
TH 
QTZ6EEQONS * TH 
THETIQOEES * TH 
JELNOTLLOS * TH 
NOL 

Th 
9960996999 Tn 
TH 
STHQSTTOGL* Th 
TSEOLHTLEQ* Th 
TH 
HITSTLEOOT “eh 
NZOLSSLLOT ‘2h 
€OZLOZOSTT* 
OSTSZ0geeT 
HILQEETTET 
NOBLOCEBST 
ZOSSLHLZRT 
eQzonenset 


69T9LOLELL* OH 
‘On 
Onetegngeé* on 
On 
*T 
626S2620TT* TH 
LILLQ2EQTI 
eQtéegétet *th 
TORLSGHTET* TH 
OLTEEQHEET TH 
QTQTOSELNT 
TH 
nLneSQsggT 
TH 
SQTESEGLET 
LITEQLZEET* TH 
th 
Th 
ENQTSSLEES "Th 
TH 
th 
S6SE20S9ge "Th 
* TH 
ELELLTIQESE* TH 
Coetnsennt 
€SLQ09e99E * TH 
Th 
OLETET9EOH 
TH 
TeQTQeZES* TH 
nENTQOOLIS * TH 
68LS$900%09 Th 
* TH 
TLOLOTESE9* TH 
TLOLOTOTEL* 
CHONEHEELL* TH 
th 
BLESTIEIER 
9S90LEQOOT 
NEQELLOLOT 
LELELELZET* eH 
LESELIOTET 
9S6T900ONT 
60S22996ST 


29H "On 
QeLEOECEZL*ON 
ZQOTEGESOL “On 
99S9L82TTR* On 
Th 
ZOSEEQEQOT Th 
"TH 
TLOOQEQTET 
LHETHHOLET 
TONIELESHT* TH 
SegeLnenst* th 
Th 

6TQICELSET "Th 
QHLIBOELOE * TH 
OSTSHTQOee 
TSZOTHSHES* Th 
TH 
TES6ETENGE 
OT2SennTee* th 
Th 
GQTREQOQTE * TH 
TH 
TH 
TSTLETLEGE "TH 
62H TOSLON* Th 
TH 
th 
O6EELEOTEN’ TH 
Th 
*HECELTEES* Th 
LSLOEQETED* Th 
TH 
OR90SHELTL* TH 
Th 
TH 
OTHLTELOLE "TH 
O660LENGCE* TH 
ETLLELQOBE* TH 
OT9S22LSOT*2h 
ESTTEEHOCT 
LESELLOGET* 


060955529 "On 
QIST2STIOL 
On 
* On 

"On 
TEESEEENLE On 
OF 
HTSEETOSOT* TH 
QBLTQZETIT’ Ih 
SEn6ECOQTT Th 
COESEQTSET 

TT9SQOGONT "TH 
TEESOQEENT' TH 
SOgNETSQST “Th 
TH 
LOEQOISQLT "TH 
OLETHOSEQT “TH 
* TH 
TOQSOSLETS* TH 
nSnLSnoLee* Th 
LLETEOe TH 
oLntegegse* th 
TH 

LELIELLLOE * Th 
TH 
agegetogne * Th 
TSTOQOZOLE "TH 
OLS6QOTETH TH 
TH 
Th 
6LELLSSS0S* Th 
ZESSALEGES TH 
LEOLETEOTO* Th 
25£SE96069° TH 
E2TEELLEED* Th 
Th 
g9Ret ESReL*T 

Sengoz6ong* 
TH 
S0S09LE69S6°TH 
LIMLOZT2OT 
ELOLQEEGOT * 
ETOLLS2H2T 


T9QES “On 
"ON 
on 
OSHELSZOTL* On 
On 
On 
OTHENZONER 
ZEQRECELNE 
TeeESELOOT 
TH 
EQTOOTRETT Th 
QSOOOT9ETT 
HEOTONAIET * TH 
ogsagesnet* Th 
OOTRIOLENT 
NECOIONTST* TH 
QLLEEQGOOT "Th 
SSLTMTSOLT 
TEQLEGEORT 
TLELTOTZET* TH 
TH 
LOGeSOLOTE* TH 
TH 
CELQELTOLE* TH 
TH 
egeeLotett ‘th 
ZLONTERTEE TH 
TH 
G29QSLASLE* TH 
* TH 

TH 
ELOTSAITEH 
TTELSRETS * Th 
6LEZLOZING* Th 
ZLEQETETES 
Th 
£219S£S099°TH 
EQQTOZONOL* TH 
696LT2SOSL* Th 
*T 

TH 
TH 
TOLSSLETLE* TH 
ETEEEQQLOT 
NELLTS9OTT 
2OTELETQTIT 
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TEOTIE ROSE 
62STETIOLE* On 
QO9TSSSTEE 
OLEQEELETH* On 
On 
nS62TTTZ9n 
TTEOLTHOTS * OF 
LETZLZO9NS * OF 
*on 
6899"B90T9 "On 
* On 
29S90LEER9* On 
QROLSHOEeL* On 
2LEQOETSOL* OF 
EQTOLENLOE 
QtELEgL096 "On 
Th 
6QnOLSLLOT*Th 
Ty 
TH 
SEENETEENT* TH 
ZOQTESSEST 
622Q0SLT9T* Th 
HIZ2QTSTLT* Th 
TH 
*Th 
* TH 
HON6CETLIS* Th 
TH 
TH 
Th 
6QOELELSLE* Th 
th 
QSSTCE6OTE*TH 
THeggseoee Th 
TH 
LEQGHLEZLE* TH 
ELOMHONIEE * TH 
Qeség_eten ih 
ST9SQELOLA* TH 
GO6ETEI9OS * TH 
GEGESRERES * Th 
* TH 
TH 


LeTggssgte 
"On 
HLTOLSQS92 * On 
602960L0ge * Oh 
“On 
ETESQIOETE* On 
ZOLONGIOEE * OF 
OLSHSELENE “ON 
LETHIQOEIE * On 
OLOLHOOOEE * OF 
On 
ORBEIEDSEN 
* On 
LTLQEELATS * On 
QHZOTT HANS On 
O9T6EORSLS * OH 
6069060609 * On 
On 
ESEQTEQTRD "On 
SQESENSTEL*On 
On 
* OF 
SOLN099SS8 "On 
On 
LEQOTTLES6* On 
STSHELOLOT*TH 
eegengontt 
TH 
Th 
Q62CEELSET* TH 
TT96EQ6ENT TH 
ZOSTENGEST* Th 
TH 
62OnEQITLT* TH 
€99S6TTeQT* Th 
* Th 
6£€QETEOSO? * TH 
OOLCOOGLTS * TH 
6STLOLE0ES* Th 
Téegstesne’ th 
Th 
Th 
Th 
TEOZOLOSTE * TH 
QOE9SHSTEE 
QBLIEHESSE “TH 
9eTS60SHLE * TH 
96260n 


HIQT6ELGIT On 
TLITEQSLLT 
"On 
TO99EQ9LET 
SETLET9IGO? * On 
SOTLTTHeLe* On 
SHE’ On 
9QLSSH06S2 "On 
* Oh 
OSTESTTSOE “On 
ELEZETSONE “On 
2LELEELESE *On 
6ZHLSTZORE 
2ESEENQTON 
SOO9TLLAZH 
Ln 
8SSE99020S 
2TONSTOTES 
2THESOLECY 
LOELLEZS99* On 
HSQOETTHOL* OF 
TOSTETESHL* On 
eLgsoenstg*on 
LS2S069ngg "On 
"On 
OeTSHESZE6 
enSSeStsot 
ESHOETHITT’ TH 
S6GR9LOQTT "Th 
SQOLTSTSET*TH 
TH 
6QECTOTENT 
TESTLETEST 
ZOeTLOBLIT* TH 
TH 
* Th 
6269E6LETS* Th 
TH 
Th 
TS6SSELHSS* TH 


é9gnetttet 
6S069200ET* On 
HHEQOOTLET ‘On 
6ELEGLSHAT “On 
LTLALLHZST* On 
LO6S62gO9T 
6£0S6S969T “On 
“On 
TILIGIBBST 
TSQHOTEBET “On 
gensstote on 
HZO9TOHHES 
On 
SEEQECETSZ “On 


*On 
SSOLLTES2E 
ETEQQELENE 
Q60S0S2E9E 
GeLLOZEEBE 
HIEHTEQSOn 
* OF 
QLOTLE9ELH 
EHETLEELOS 
OTLEENIIES “On 
* On 

On 
98e990E S69 On 
ELECTELGHE “On 
TH 
HSLOLIAGOT* TH 
Th 
QO99EESETT 
OTeLTQLICT 
Th 
SE0SSnSent Th 
LOQ6EQTIST* Th 
EQO9TLEOST * Th 
H9OOSHTOLT TH 
QeQLEnSOgt 


EE 
THQLETEHOT 
NOETSELSTT 
* On 
OLT9TCLSET "On 
LLSTQOTENT 
SOLE9EGOST 
Q9GBCTCEST “On 
S66Q296L9T *On 
E992STELLT 
EQTTOOOLET 
E99HENELET “On 
* Oh 
"On 
O6eHSOTSHS 
SStLetggse on 
LQOCELALEE 
QLEELEOSOE 
*On 
T6ZOL TEESE * 
6LLOQTHOGE * On 
"on 
“On 
L6SL9LEZ0S On 
ZTS9SIETES *On 
* On 
692@RSESSES 
TIQOEHEOLD "On 
LELETQOLII* On 
QOTIEITIOL* On 
QCOTLE9TEL*On 
LEZESSHEEG 
SenggoTege 
909SENGONE 
QT9S098966 On 
2E6ELEQSOT “TH 
TH 
LLOOTEHEET* TH 
TH 


QTEQTEELOQ’ EE 
TH * EE 
190606296" 6£ 
HSEOSLATOT “On 
SS0SO0SE90T "On 
690NHELETT On 
ALLENTEBTT “On 
OSTLECTZET 
QHTCHELEET 
S9OLSSE9HT “On 
TQS9ZTOSST 
QQSOTESEQT 
HLTSQESELT 
SEBLEETCET 
* On 
6OHELLONTS 
HNOTETEQQES "On 
On 
OH 
LEQLETTTGe "On 
LLOQELLETE* On 
SSeTEOHOSE 
ZTONETSTEE 
“On 
COLELEILER* OF 


"On 


6LEETTELTS 
OTSHONEL' 
6TLENOBELS * On 
6SSLS9n6n9 Oh 
ELnLEeSLeg "on 
HOTEQTELEL* On 
HENLEELOLL* On 
On 
On 
tHENCOLSOLE 
“Th 


EQLOOTOGOT TH 


*6E 
6E 
EE 
OMMEQLETTL*6E 
*6E 
296EQTREGL* 
"EE 
6S6LT9S9LE*6E 
*6E 
ESOOOLZELE* EE 
6Q00ELSZOT “On 
STEHT IQOT ‘OR 
GENELIELTT 
9ge29eSTO2T 
SHBLLEQIET “OF 
2QSTHESEET On 
ESLLEEGONT * On 
SETOOTOQNT "On 
QT2LLOLIST 
ZEQZLEESOT “On 
ON 
EQSQTITHGT "On 
OH 
HEQLT9TSO? “OF 
* On 
On 
266LECOSSE 
OW 
On 
€QO9LESOOE *On 
TETCEOGLTE 
CEEZELQSEE 
STETORLHSE “On 
“On 
O9TT9STETH 
on 
QLLTI9ELE9N * OF 
egtgegeséh "On 
QTEBLOGHES “On 
9991510555 
66L696£L9S 
On 
* On 
"On 
SOENQEILEL* On 
On 
ESELEZSOLQ "On 


2TTZEQOTeS 
*6E 
LHQSETELLS *6E 
LEQHOLLLOO* 6E 
SSNLLEQEED 
QOLERISEOL* 6E 
EE 
OBELSLTLEL*6E 
OSLS6E262g *6E 
EE 
6920675026 *6E 
ETELELOOLE*6E 
O2LLLEZZOT 
O9OLE9LLOT* On 
LELEEBSETT 
CESTOOLETT 
ZTHSEQ2deT 
QESLEQTEET* OM 
“On 
NESLLLEQST “On 
QQ9ESLEELT 
OTOE9ESEQT “On 
EQERELLEET 
* On 
TSEOLHLONE “On 
Z60geSenSe*on 
On 
“On 
* On 
O9R9TOLOTE*On 
TROBLEGHEE 
SOLLENLESE 
LLSOLLEELE * On 
ZENHEZESEE “On 
LENELOORTH* On 
ESOTLHOeHH On 
SONENESLON* On 
OT2EO9SHEN 
QLTON29LSS "On 
°On 
S6E960029 "On 
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general, for given argument, this formula was used calculate two 
the functions, with orders differing unity, and the recursion formula 


was then used compute the other orders required. 

For half-integral orders, the series equation (2) exact, and difficulty 
arose the calculation the corresponding tables. For integral orders the 
terms the series eventually diverge; and order cope with this, three 
separate ranges the argument were considered follows: 


this range the computer calculated and and built the other 
orders using equation (3). The asymptotic series was summed the point 
where the terms began diverge, and the series was then rounded off 
adding the sum one half the first divergent term. Within its range this 
approximation gave almost 12-figure accuracy the results. 


Over this range method due Dingle (1958a, 1959) was 
calculate more accurate converging factor than the crude value 0.5 used 
the first range Generally, this method, the asymptotic series 
allowed run for number terms beyond the first divergent one, and 
cut off after terms, where depends upon the values and 
The converging factor which applied the last term calculated, may 
written 


where the are defined the form 


r—8 
r=t 


The numerical factors and are given Dingle 1959) 
terms 3.* Values and were found from the relationship 

(6) 

(7) 

where accordance with the derivation equation (5), restricted 
small values; 0.5. Since and must integers, was chosen 
the half-integer nearest was then obtained from equation (6). Because 
the series for the converging factor itself asymptotic, was necessary 
insert into the computer program test which was arranged cut off 
correctly its least term. For this range the best results were obtained 
calculating and and using the recursion formula obtain the 


*In the calculation the tables, the expressions equations (4) and (5) were used 


i 
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other orders. Above the method gave results accurate least 
figures; but from this point down 3.6 the accuracy fell off steadily. 
the region the error the tables may amount unit the 
10th figure. 


This the only range which method independent equation (2) was 
used calculate the Bessel functions. The method adequately described 
Aldis (1899) and will not repeated here; makes use series which occur 
the calculation the modified Bessel functions the first kind. this 
range and were calculated basis for the application equation 
(3). The accuracy the method limited essentially the fact that the two 
basic are derived the differences other series each which could 
only calculated the computer limited number figures. Once again 
the region 3.6 the absolute value the error may approach 
one unit the 10th place; going lower values however, the accuracy 
improves until, the region 0.1, 12-figure accuracy obtained. 

order estimate the limits the various ranges the argument, and 
decide the magnitude the errors within the ranges, sets values 
and were computed the various methods described above, for 
12; these were compared with table more accurate values 
and over the same range, given Aldis (1899). These computa- 
tions served also compare Dingle’s method calculating the converging 
factor for the series with simpler method due Burnett (1930); 
this work forms the basis separate note (Dempsey and Benson, unpub- 
lished). 
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ANALYSIS SOME STATISTICAL PROPERTIES 
AURORAL RADAR REFLECTIONS AND THEIR RELATIONSHIPS 
THE DETECTION CAPABILITIES THE 


ABSTRACT 


statistical model auroral echo occurrence has been made from analysis 
observations 48.5 Mc/s obtained over number years continuous 
operation. The probability density distribution auroral target cross sections 
has been examined experimentally, and the resulting curve fitted simple 
mathematical relations. Both inverse power law and exponential law have 
been derived, the forms 


and 


= — 


These models have been interpreted terms distributed and localized targets, 
and used analyze the echo occurrence indices and the effect which variation 
radar parameters will have upon them. Both forms the target law are useful 
although considered that exponential form yields better agreement with 
observations over wider range the variables. 


INTRODUCTION 

previous paper (McNamara 1958) has given the details auroral 
radar designed for continuous unattended monitoring radar echoes from 
aurora. Two complementary forms record are obtained from the radar and 
are scaled give measures echo occurrence. This paper examines the 
parameters involved the design equipment for auroral detection, and 
their relationship the performance, echoes obtained, range distribution, 
scaling procedure, and indices derived from the data. 


RADAR EQUATION FOR AURORAL ECHOES 
For the jth single elementary target, small dimensions, the received 


(1) 


where refer the transmitting and receiving parameters respectively and 
function proportional the magnitude the Poynting vector 
the incident wave such time that the reflected signal from 
the scatterer returns the receiver the instant under 
sideration, 
power gain the antenna over isotropic radiator, 
radio-frequency wavelength, 
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pattern propagation factor the field strength, 
antenna pattern function expressing the relative free space field 
strength the direction specified elevation and azimuth angles 

interference function, between the direct and ground-reflected rays, 
giving the field strength the specified elevation relative the 
maximum field strength free space, 

equivalent backscattering area the elementary target, 
range the target, 
the aspect function the target cross section, giving the factor 
which the effective target area reduced when the target located 
the direction 


From consideration the characteristics auroral echoes, the actual 
targets are thought represented number scatterers which scatter 
non-coherent manner relative each other. The received power the 
arithmetic sum the powers from individual scatterers, and averaged over 
short period time (the integration time the eye, the camera viewing 
the cathode ray tube, other recording device). Observations (Currie, For- 
syth, and Vawter 1953; McNamara and Currie 1954; Unwin and Gadsden 
1957) also indicate that the auroral echoing region relatively thin layer, 
over which the elementary targets will considered distributed. 

For non-isotropically scattering targets, the returned signal modified 
‘aspect attenuation’ function which dependent upon the angle inci- 
dence the radiation asymmetrical targets and hence upon the geometry 
the radar case. The aspect function will denoted where 
the angle intersection with the major axis the target, which assumed 
have axial symmetry about the earth’s magnetic field, and relative 
the radar. constant height reflection, the aspect function can 
expressed terms range and azimuth from the radar the function 

Radar observations Mc/s show that echoes can received angles 
deviating much 15° from orthogonality the magnetic field. Hence 
the length field-aligned concentrations generally will not greater than 
several tens meters such beam width observed the scattering 
pattern the target. This dimension much less than Fresnel zone and 
the range dependence the scattered power from such elementary targets will 

The horizontal distribution targets may now classified under the two 
conditions: 

Case (7): The scattering region completely fills the azimuthal beam width. 

Case The scattering region more localized space, occupying only 
fraction the azimuthal beam width. 


both cases the extent range the scattering region will assumed 
appreciable, i.e., greater than the pulse lengths commonly employed. 
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The scatterers are considered independent and the received power, averaged 
over many pulse periods, the sum the average powers from each contri- 
buting scatterer. 

For case the density function now introduced. This 
the number scatterers per unit area cross section the 
scatterers are uniformly distributed horizontally, then the density function 
separates the integration yielding 


0 
the average cross section per unit area the ionosphere. Hence the signal- 
to-noise ratio 

(2) 64r°kT, B de (9, (0,¢)do dR. 


Now “beam width” may defined 


Within the next integral, non-zero only over short interval roughly 
equivalent the pulse length. Over this small interval the other factors 
the integrand are substantially constant and may moved outside the 
integral. There remains 

vi) 
where the definition terms the conventional values 
peak power transmitted, 
pulse duration, 
velocity propagation the electromagnetic energy. 


total noise power referred the receiver input terminals, 
Boltzmann’s constant, 
effective absolute temperature the source, including the effects 
system noise, atmospheric noise, and cosmic noise, 
noise bandwidth the receiver. 
Collecting all the parts 
Turning now case and starting with equation (1), the ensemble 
scatterers may replaced equivalent scatterers time and azimuthal 
average cross section located the centers the active region. 


Also 
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Since shall interested following the long-term variations 
through P,, shall introduce the approximation that the propagation 
patterns and aspect factor are broad and flat the azimuthal direction and 
cut off sharply beyond some critical azimuth. Then for fixed broad beam 
radar, equation (6) reduces the approximate form 

This formulation will also applicable rotating narrow beam. 

Other factors being constant, will noted that the signal-to-noise ratio 
proportional where the factor arises from spatial integration 
the distributed target the pulse length. For optimum signal-to-noise 
the receiver, the bandwidth should approximately equal the inverse 
the pulse length and the signal-to-noise ratio becomes proportional 
Hence peak power may traded for narrow bandwidth and large high-power 
radar may replaced compact, reliable, low-power equipment. 
example, for completely distributed target, 1-kilowatt pulse radar 
can approach the detection sensitivity 100-kilowatt pulse radar. 

For any given radar the received power sensitivity determined the 
antenna system through fourth-degree propagation pattern. Since practice 
all propagation patterns, both for horizontal and vertical polarization, have 
null zero degrees elevation, the detection echoes low elevation 
sharply dependent and hence antenna height. antenna height 
made too great, however, nulls the pattern may appear important echo 
elevations. 


DEFINITIONS AND INDICES ECHO OCCURRENCE 

The most complete and detailed description echo occurrence conveniently 
obtainable from range-time recording with fixed beam radar obtained 
method whereby the record searched units which are 
one equivalent pulse length range one unit time, and occur- 
rence defined auroral echo detectable this cell. These can also 
grouped under either the range intervals time obtain statistics varia- 
tion occurrence with range percentage time occurrence echoes 
respectively. The latter quantity, the occurrence any auroral signal any 
range unit time interval, advantageous that applies uniformly 
any equivalent pulse continuous wave system and much more quickly 
read from the film record. Moreover, this index can also obtained directly 
and immediately from the chart record produced the signal integrator 
(McNamara 1958). The integrator also yields another index which measure 
intensity the disturbance. For daily measures the area under the trace 
the time variation signal amplitude, integrated over all ranges, may 
scaled. the daily basis, plots the experimental data percentage occur- 
rence and integrated area, both derived from the integrator charts, show 
fairly well-defined statistical relation between the two indices, although these 
indices, occasion, may expected differ considerably. The relationship 


McNAMARA: AURORAL RADAR REFLECTIONS 429 


demonstrated Fig. where 423 daily points from the period April 
1956 May 1957 have been used. The data well fitted the parabola 


(solid line) Fig. The dashed lines give the values the first-order standard 
deviations. 
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Fic. The statistical relationship between 423 daily indices echo intensity and duration 
occurrence. The solid line the line best fit and the dashed lines are the first-order 
standard deviations. 


DERIVATION ECHO OCCURRENCE FORMULAS 


discussed Section observations indicate that the auroral echoing 
region relatively thin layer over which will consider the long-term 
unit-area target cross-section probability density function given 


(number time units which target 
with cross section between and 


(total number time units 
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will implicitly assumed here that the probability density function 
independent the geographical position over the ionosphere. This least 
reasonable over the limited region viewed radar. Geographical variations 
the frequency occurrence are accounted for multiplicative factor 
which represents the average rate occurrence all targets per 
unit area the ionosphere. Therefore, the number targets observable 
per unit area over extended period time proportional 


o 


(9) 


where the limiting cross section for threshold detection the radar. 

the geometry thin auroral echoing layer height over curved 
earth, may shown that the occurrence targets range interval 


eRe / 


Ri 


and for case 


COS 


JR; 
where radius the earth, 
angle subtended the center the earth the arc from the 
radar the subtarget point, 
and (h/p) for ranges which echoes occur. 

From this point onward, shall only carry out the development case (77). 
Case (7) then obtainable from any the case (77) equations deleting the 

The minimum detectable signal taken constant given 


From equation (7) 
But 


or 


(13a) since the integral goes zero the upper limit 


and therefore 
(14) 
eR, 


| 
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predict the range distribution, small range segments are taken. the 
range segments are all equal and sufficiently small,then the mean value 
theorem the range distribution well represented the integrand 


The term which most uncertain due the aspect factor con- 
tained within it. principle should possible plot the distribution 
setting and estimating values from other sources, and 
using the difference between this curve and the observed range distribution 
derive the aspect function 

From radar observations, the aspect factor appears tolerant small 
departures from normality but falls rapidly with increasing departures pro- 
ducing rather sharp cutoff echo occurrence some minimum range. 
Theoretical derivations aspect function have been made several 
authors the bases definite scattering models, and these could sub- 
stituted the above relations. However, for determination certain statistical 


properties total echo occurrence assumptions need made physical 
model, and simple linear approximation can used represent the attenua- 
tion. For illustrative purposes, will represented here the ramp-like 
function 


A(R) 


| R> Ri 


The cutoff range and the ‘full-signal’ range are suitably chosen with 
reference the observed range distribution examination the aspect 
angle (for example, Figs. McNamara 1958). Although this results 
inaccuracies range distribution near the cutoff, the integral total 
occurrence not greatly error are well selected. 

The function expressing the true geographical variation mean 
auroral occurrence not well established present. necessary know 
the zenithal occurrence each point. Again, shall make simplification 


(16) 


EXPERIMENTAL AND THEORETICAL PROBABILITY DENSITY 
FUNCTIONS FOR TARGET CROSS SECTIONS 

The law which governs the distribution target size required specify 
the nature the function used the preceding section. this end, 
calibration the over-all system was made with random 
noise generator the receiver input since this type signal most repre- 
sentative auroral signals. The recording characteristic was thus established 
nearly linear with noise power input, and hence linear with echo cross 
section Echo data were derived from 469 days consecutive operation, 
January 1957 April 1958. The function scaled was actually the inte- 
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Fics. and 2B. Experimentally derived target cross-section distributions fitted power 
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grated probability density, defined equation This was obtained 
determining the time that the signal level was above series selected 
levels the chart. 

Mathematical description the cross-section distribution was obtained 
fitting simple expressions the data. Figures and show the results 
assuming inverse power law and exponential law. Reasonably good 
fits are obtained over the dynamic recording range for both forms. the 
following the consequences both approximations will examined. Further 
discussion the applicability the laws will then given. 

(a) The inverse power law: 


Hence 
K’ 
L 
and 
(6) The exponential law: 
where mean cross section. Hence 
and 


The slope the best-fitting straight line Fig. yields the value 
From several selected periods, the quantity m—1 found lie the 
region 0.5 0.9, and representative value 0.67 will adopted. This 
result, 1.67, consistent with earlier work the author (McNamara 
1954) Saskatoon where the value 1.5 was obtained using much 
different method analysis data obtained from high-power short-pulse 
radar Mc/s. the present instance, additional checks were also made 
shorter periods operation with simultaneous recordings 50-km gate length 
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covering the interval and the standard 800-km gate length 
covering 300-1100 km. The distributions obtained over the period December 
1958 January 1959 with short and long gates are shown Fig. The 
close correspondence the slopes the two curves supports the validity 
the application the long-gate statistics the general analysis. further 
check the exponent was made using logarithmic receiver recording 
30-db range signal, yielding the value 1.5. 


---- GATE 


SCALE, ARBITRARY UNITS 


Fic. Comparison cross-section distributions obtained from simultaneous operation 
narrow and wide gates centered range The slopes m—1 the lines are 
and 0.53 respectively. 


Equation may thus taken giving the variation echo occur- 
rence with range, and giving measure the functional dependence 
occurrence rate the geometry, equipment parameters, and antenna pat- 
terns. Plots equations and for two different combinations antenna 
heights are shown Fig. corresponding case (77). Modification equa- 
tions and 22, deleting the multiplicative factor and changing the 
results equations appropriate case (7). These are plotted 
Fig. For all sample calculations the values N(R) 
were used. each figure, the curves correspond the 
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Theoretical echo range distributions for two radar systems differing only the 
heights their transmitting-receiving antenna configurations above the ground. Calculations 
are based the assumptions completely distributed targets, and for power law (A) and 
exponential law (B) target cross sections. 


inverse power law distribution target cross sections, and the curves 
correspond the exponential law. 

For the inverse power law, the shape the echo distribution independent 
the radar sensitivity factors, the vertical scale varying the 0.67 power 
the sensitivity. For the exponential law, the shape well the vertical 
scale varies with the sensitivity, and also with the mean the target cross 
section distribution. 

indicated Tables and both forms yield distributions and occur- 
rences not unlike the observational data, provided appropriate choice 
made for the exponential law. Additional data from (A/2) system 
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TABLE 
Echo occurrences derived from Fig. for the case completely distributed targets (case 


Relative Area under curve* Peak value first lobe* 
mean target 
cross (2d) (A/2) Ratio (A/2) Ratio 


Power law for cross section 


100 2.9 290 100 2.6 260 

Exponential law for cross section 
100 2.2 220 100 1.5 150 
100 240 1.5 360 150 1.3 200 


*Normalized to the value 100 for the reference configuration (A) (2A) (A/2), (B) (2A) (A/2) with om = 10. 


TABLE 


Relative Area under curve* Peak value first lobe* 
cross (2d) (A/2) Ratio Ratio (2d) 
configuration configuration configuration configuration 
Power law for cross section 
100 2.4 240 100 2.4 240 
Exponential law for cross section 
100 1.8 180 100 1.4 140 
100 210 1.6 330 150 1.2 180 


*Normalized to the value 100 for the reference configuration (A) (2A) (A/2), (B) (2A) (A/2) with om = 10. 


with 20-db variability sensitivity suggests that the exponential law can 
fitted more accurately over larger dynamic range. Mathematically, the 
exponential law may more convenient since all its moments are finite and 
can normalized. the other hand, the inverse power formulation may 
often convenient for comparison system parameters since assumptions 
nor measurements need made concerning the magnitude the cross-section 
scale. 

should noted that although the equations contain radio wavelength 
term they not give the frequency dependence the echo occurrence. 
The term arises from the conversion effective area receiving antenna 
into antenna gain; the true wavelength dependence the echoes hidden 
the cross-section term, and not considered this paper. 


INTEGRATED ECHO OCCURRENCE 

The formulas the preceding section were developed terms the detailed 
echo occurrence defined the range and time division discussed Section 
Also interest, the description occurrence terms time alone, any 
echo echoes present given time satisfying the definition occur- 
rence. The formulas may used guide deriving these indices but 
cannot confidently applied because the unknown correlation echo 
occurrence space and time. fact, the geographical correlation distance 


Echo occurrences derived from Fig. for the case localized targets (case 
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target occurrence often large and the probability echo occurrence 
adjacent range increments not independent. the echoes each range 
increment were statistically independent time, the range-integrated echo 
probability would proportional the area under the curve. 

This may seen from the following argument. Assume that echoes 
adjacent range increments occur independently time and let the probability 
that range increment will have echo and the probability that all 
increments will simultaneously fail have echo Our 
definition echo occurrence purely time basis that least one 
increment does contain echo. This probability which, 
for small probabilities, may expanded 


This simply the area under the distribution curve equations and 22, 
for which examples are plotted Figs. and The areas under these curves 
are tabulated Tables and and may taken approximation 
the relative amount echo which will observed. 


NUMERICAL EXAMPLES 


Tables and II, corresponding cases and distributed and 
localized targets, are given the areas under the curves and the heights the 
peak the main lobe. The tabulated values have been normalized 100 for 
the reference configuration (A/2) with 10. should noted that 
the computations were made the basis uniform occurrence targets 
all over the ionosphere, 

Experimentally, for two such radar systems operating simultaneously 
Ottawa over period years, the occurrence ratio was found 2.5 
Although allowance has been made for small system gain differences, 
nor for the naturally increasing occurrence aurora with range toward the 
auroral zone, the agreement with the tabular values satisfactory. 


CONCLUSION 


Expressions have been derived relating the radar parameters the auroral 
echo statistics observed. For this done, was necessary adopt practical 
definitions ‘echo occurrences’ and interpret the appropriate activity 
indices. The sensitivity the indices demonstrated equations and 
and Figs. and and indicates the necessity for caution the comparison 
indices from different systems and the need for accuracy determination 
and stabilization the radar parameters. 
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LIQUID-SOLID INTERFACE SHAPE OBSERVED 
SILICON CRYSTALS GROWN THE CZOCHRALSKI 


ABSTRACT 


study has been made the macroscopic structure the inter- 
face which exists during growth silicon crystals the Czochralski method. 
two crystals the interface was seen contain (111) facet. The development 
such facets discussed with reference current crystal growth theories. 


INTRODUCTION 


The shape and structure the liquid-solid interface during growth 
silicon single crystals great importance and has been discussed number 
authors (Cressell and Powell 1957; Rosi 1958; and Tiller 1958). For example, 
Cressell (1957) and Rosi (1958) have observed that interface which 
planar slightly convex towards the liquid leads crystals with minimum 
dislocation density. this laboratory interfaces have been observed which are 
parallel, over macroscopic distances, the (111) planes lying normal the 
growth direction. 

The interface shape may observed number ways. may 
viewed directly rapidly separating the growing crystal from the melt during 
growth (Elbaum and Chalmers 1955). The possibility that liquid adhering 
the crystal will solidify after decantation and the fact that the interface shape 
observed one position only are the limitations this method. second 
method revealing interface shape requires the addition volatile impurities 
melt under vacuum from which crystal being grown the Czochralski 
technique. For example, silicon crystal grown from melt containing boron 
impurity will P-type, but when excess antimony added the 
melt the growing crystal becomes .V-type. Under vacuum, antimony evaporates 
rapidly compared with boron and the growing silicon crystal reverts P-type 
(Bradshaw and Mlavsky 1956). The grown crystal may sectioned longitu- 
dinally, i.e. parallel the growth axis, ground and etched reveal the 
and junctions. The interface shape the time the addition the 
antimony shown the shape the junction. this way the interface 
shape number times during the growth the crystal may studied. 
Rate-grown junctions (Tanenbaum 1955) may also used give 
continuous record the interface shape during crystal growth. However, this 
method imposes stringent conditions the growth rate the crystal and 
the impurities added the melt. 

The most useful and convenient method interface examination makes use 
the variations crystal composition which occur during growth. Unless 
special precautions are taken during growth, Czochralski-type silicon crystals 

received October 1959. 
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exhibit variations resistivity and oxygen content both radially and longi- 
tudinally. The longitudinal variations are closely related the rotation rate 
and many instances their period equals the growth per revolution the 
crystal. crystal cut longitudinally and the exposed face ground and 
etched these variations striations are revealed corrugations the surface. 
The striations give revolution-by-revolution macroscopic picture the 
interface shape during crystal growth. Striations with separation 0.001 
have been resolved. variation this method Camp (1954) has used 
copper-plating technique delineate resistivity striations germanium and 
has successfully resolved striations with separation just less than mm. 


EXPERIMENTAL 


Six crystals were examined; these were grown the Czochralski method 
with the solid—liquid interface convex towards the liquid. Each crystal was cut 
longitudinally and then the newly exposed faces were ground and etched 
described the Appendix. The corrugations produced the flat surface 
the etch had separation equal the growth per revolution the crystal and 
curvature expected for interface which was convex towards the liquid. 

Two crystals examined had two distinct breaks each the striations and 
the sections the striations between the breaks were straight lines (Figs. 
and Typical straight line sections were 0.45 long striations with 
radius curvature 2.5 cm. Near the top the crystal, shown 
Fig. la, the reduced radius curvature the interface restricted the develop- 
ment these ‘flats’ the interface and caused increase the change 
direction the breaks the striations observed maximum 73°. These 
straight line sections the striations were observed perpendicular the 
direction which was near the growth axis the crystal. Each the 
two crystals was cut again bisect the straight line striations and reveal 
second face normal the original section. These second faces were also 
ground and etched order reveal the striations. The two flat faces each 
crystal were masked with ‘Halo’ wax and the remaining crystal surface was 
etched hot caustic soda which structure-sensitive etch for silicon. 
optical goniometer was used show that the (111) planes approximately 
normal the axis growth would give traces the two sectioned faces 
parallel the straight line sections the striations. 

The growth axis one crystal was 4.5° from the [111] direction and this 
angle shown Fig. la. The growth and rotation rates this crystal 
were cm/sec and 4.6 r.p.m., respectively, which gave striation 
separation cm. The crystal was P-type with resistivity 
approximately 0.4 ohm-cm. The exact growth conditions the second crystal* 
are not known. This second crystal was with resistivity approxi- 


*Kindly donated D.R.B. Northern Electric Co. Ltd. 


Fic. la. The form the striations seen after etching section silicon crystal which 
had been cut parallel the growth axis. The lines and indicate the extent the 
(111) flats. 

Fic. 1b. Photograph the striations taken the approximate position the inset 
marked near the top left corner Fig. la. 
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mately ohm-cm and the growth axis was from the [111] direction. The 
striation separation varied between 0.001 and 0.01 cm. 

this second crystal was observed that where the diameter was increasing 
rapidly, i.e. the sections and (Fig. the striations abruptly 
changed direction within approximately 0.1 from the edge the section 
(Fig. 2). These were probably traces other {111} planes but unfortunately 
poor masking the sample led their loss during etching and before measure- 
ments were made. 


Crystal Growth 
Direction 


enlarged drawing the crystal-section edge the region Fig. la. The 


abrupt change striation direction near the surface indicated. 


DISCUSSION 

With crystals grown the Czochralski method the {111} planes often 
give rise facets the crystal (Billig 1955). Crystals pulled the (111) 
direction have three such facets spaced symmetrically around the circum- 
ference the crystal; the fourth {111} plane which normal the growth 
direction lies approximately parallel the interface. The {111} 
planes are the most closely packed the silicon structure. Growth direction 
normal these planes slow because the low accommodation coefficient 
(Chalmers 1954) whereas growth lateral extension the planesis 
preferred (Tiller 1958; Billig 1955). The preferred directions lateral growth 
the (111) plane are the [211], [121], and [112] directions (Billig 1955) which 
give rise the three symmetrically placed spikes commonly observed when 
crystal pulled the [111] direction. 

the (111) set planes lying perpendicular the growth axis [111] 
crystal which has been detected macroscopic scale the etching experi- 
ments and which normally does not give rise any visual characteristics 
crystal. two crystals the interface was found limited over macroscopic 
distances the densely packed (111) planes which are approximately normal 
the growth axis. 

The temperature gradient existing the interface during crystal growth can 
estimated only approximately. Calculations made Billig (1955, 1956) are 
inapplicable because the considerable radiation received the ingot from the 
adjacent melt and susceptor ignored; this leads values temperature 
gradient which are much too high. With growth from crucible the temperature 
gradient decreases crystal growth proceeds. Observations made with 
pyrometer, pertinent times during crystal growth, indicated temperature 
drop 150° the first centimeter from the interface. 
assumed that 


4 
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(a) the gradient actually existing the interface 150°/cm, 

the curved portions the interface are coincident with the 

freezing point isotherm, and 

(c) the interface has approximately constant radius 

curvature, 
then the two crystals considered there approximately 1.5° super- 
cooling the center the (111) flats. 

Nucleation (111) planes will occur this point. With interface slightly 
convex towards the liquid these (111) planes will grow outwards stepwise 
manner (Tiller 1958; Elbaum and Chalmers 1955) the limiting diameter 
the crystal. the interface appreciably curved, Fig. la, the completion 
(111) layer will take several revolutions the growing crystal. However, 
the interface concave towards the liquid there the possibility nuclea- 
tion several points around the interface periphery. With nucleation the 
periphery the (111) planes grow inwards stepwise manner completion. 
the interface periphery there liquid—gas and interface. 
Twinning and heterogeneous nucleation are likely occur here, possibly 
because the lattice spacing not quite normal the last few atomic layers 
near the interface and also because the presence minute 
particles oxide other impurity, for example, silicon carbide carbon 
susceptor and heater are used. The growth strays and twins into the main 
body the crystal thus encouraged the concave interface and this has 
been discussed Chalmers (1953) and Bolling al. (1956). 

conclusion, lattice misfits and twinning are more likely occur (111) 
planes nucleate the crystal periphery and grow inwards rather than 
nucleation occurs near the crystal center with growth outwards completion 
the case with interface which slightly convex the melt. Results 
have been given which suggest that the (111) planes normal the crystal 
axis play important role the growth mechanism. equilibrium 
the interface will composed {111} planes they have the lowest accom- 
modation coefficient and equilibrium temperature. Consequently, interface 
limited the (111) plane normal the growth axis most likely occur 
crystals which are pulled slowly from melt which the temperature gradients 
are small. This picture crystal growth, which especially applies the two 
instances cited, agreement with the observations Cressell and Powell 
(1957), Rosi (1958), Bolling al. (1956), and Billig (1956) that the least 
twinning, lineage, and the lowest density dislocations occur when the 
interface planar slightly convex towards the liquid during crystal growth. 


APPENDIX 


The etching silicon mixed nitric acid etch 
basically two-step process. The oxidizes the silicon and the dis- 
solves the oxide re-exposing the silicon surface the action the 
(Robbins and Schwartz 1959). The oxidation the structure-sensitive step 
and suitable control the content leads maximum information 
regarding the structure and composition the silicon. 
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For the delineation striations the sample surfaces studied were first 
ground and then chemically polished mixture parts HNO; and 
part (48%). The sample was then placed and 
approximately 0.03 (one drop) HNO; (70%) was added the solution 
every seconds for total time from minutes. the 
specimen during the etching process greatly enhanced the delineation the 
striations. The rate addition drops and time etch were not critical and 
the most suitable conditions for silicon given resistivity were experiment- 
ally determined. slow oxidation step appeared essential and the 
sample was immersed comparable HNO;-HF mixture fast reaction took 
place with evolution considerable heat and very poor, any, delineation 


striations. 
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EMISSION SYSTEM THE MOLECULE! 


ABSTRACT 


The “methyl iodide flame lying the region 4100 6300 have 
been photographed using 21-ft concave grating spectrograph. The bands are 
shown arise from transition the molecule. Rotational 
and vibrational analysis the bands has been carried out and the molecular 
constants obtained. 


INTRODUCTION 


number red-degraded bands, now known the iodide flame 
was first observed Vaidya (1937) emission from flame 
methyl iodide mixed with alcohol coal gas. The bands were analyzed 
into two systems, and and were tentatively assigned the 
molecule. Coleman, Gaydon, and Vaidya (1948) obtained the bands 
emission from oxyhydrogen flame which iodine was added, thus support- 
ing the assignment These authors observed approximately bands 
the region 4100 6300 and showed that all the bands could arranged 
into one system instead two. Vibrational constants were given but 
rotational analysis was attempted since the bands were observed only with 
low resolution. 

Durie and Ramsay (1958) observed six bands the same system absorp- 
tion during the flash photolysis mixtures iodine vapor and oxygen. The 
bands were observed with high resolution but rotational analysis the 
absorption bands was attempted since the same bands, extending much 
higher values, were also observed emission from oxyhydrogen flame 
which methyl iodide was added. Most the bands reported Coleman, 
Gaydon, and Vaidya were rephotographed under high resolution Durie 
and Ramsay. The rotational and vibrational analysis these bands forms the 
main subject the present paper and confirms the assignment the bands 
the molecule. 


EXPERIMENTAL 

The source used excite the methyl iodide flame bands consisted 
diffusion flame hydrogen, saturated with methyl iodide vapor room 
temperature, burning atmosphere oxygen. The bands were photo- 
graphed using the second order 21-ft concave grating spectrograph and 
Eastman Kodak 103a0, and plates. The exposure times 
varied from hours. iron are was used provide reference spectrum. 
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reproduction the spectrum given Fig. enlargement showing 
the rotational structure one the bands and the rotational assignments 
given Fig. may seen from Fig. that many the bands have discrete 
and well-resolved rotational structures while others are completely diffuse 
show distinct broadening the rotational lines. The presence diffuseness 
certain bands due predissociation the excited state. Using the 
vibrational assignments Coleman, Gaydon, and Vaidya find that bands 
with and are completely diffuse, bands with have rotational 
lines which are slightly broader than the slit width, while bands with 
and have rotational lines which are sharp. These different types rotational 
structure may seen (Fig. comparing the rotational structures the 
(4,0), (3,0), and (2,0) bands. 

Wavelength measurements were carried out using comparator equipped 
with photoelectric scanning device the type described Tomkins and 
Fred (1951). Standard wavelengths were taken from the wavelength 
tables (Harrison 1939) and were combined with the vacuum corrections 
Edlén (1953). The rotational lines eight bands were measured and the 
vacuum wavenumbers and assignments are given Table The wavelengths 
the sharp band heads are given Table the band head measurements 
Coleman, Gaydon, and Vaidya are also given. Both sets measurements are 
compared with the wavelengths the band heads calculated from our final 
molecular constants. 

ANALYSIS 

The electronic transition responsible for the bands expected 
the type analogy with the known transitions for ClO and 
BrO (Durie and Ramsay 1958). For Hund’s case coupling, each band 
splits into two well-defined subbands, viz. and and each 
subband consists predominantly single and branch and very weak 
branch. For the discrete bands show single and branches but 
branches were detected. Moreover, for each band, only one the subbands 
was observed, presumably due large spin splitting the upper state. This 
not surprising view the large doublet splitting the ground state the 
iodine atom, viz. 7598 has been assumed that the observed subbands 
are the type although should pointed out that was not 
possible verify this assumption observation the and branches 


the lowest values. 
The various molecular constants were determined the following standard 


procedures (Herzberg 1950). 

The rotational constants and were obtained from the combination 
and and were obtained using data from more than one band. These 
values are considered more accurate than those obtained for the other 
vibrational levels and are given Table III. The other and values were 
refined the manner described Sections and 
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TABLE 


Vacuum wavenumbers and rotational assignments for the subbands 


(0,2) (2,2) (2,0) (3,0) 
J R(J) P(J) R(J) P(J) R(J) P(J) R(J) P(J) 
2.5 
3.5 21130 .25 
4.5 129.10 
5.5 127.88 
6.5 20130 21134 .09 126.49 22470 .76 
7.5 129.27 133 .65 125.00 470.24 
8.5 20137 .65 127.74 133.09 123.38 469.64 
9.5 137 .03 126.10 132.46 121.64 468.87 22951 .26 
10.5 136.46 124.22 131.63 119.77 467 .95 948.86 22937 .25 
11.5 135.56 122.40 130.72 117.77 466.89 947.77 935.10 
12.5 134.67 120.34 129 .67 115.63 165.71 946.46 932.69 
13.5 133 .59 118.26 128.48 113.38 464.39 944.99 930.2 
14.5 132.45 115.98 127.18 462.92 943.49 927 .52 
15.5 131.14 113.62 125.74 461.32 941.83 924.85 
16.5 129 .83 111.15 124.16 459.57 939 .96 921.86 
17.5 128.32 108.57 122.49 457 .70 918.80 
18.5 126.74 105.88 120.69 455.67 935.84 915.64 
19.5 124.98 103 .02 118.72 453.52 933 .55 912.27 
20.5 123.15 100.07 116.66 451.21 931.15 908.81 
21.5 121.17 097 .00 114.46 448.80 928 .57 905.17 
22.5 119.08 093 .83 112.15 446.24 925.83 901.40 
23.5 116.91 090 .53 109 .68 443.51 922.98 897 .46 
24.5 114.57 087.15 107.11 440.66 919.99 893 .36 
25.5 112.13 083 .61 104.40 137.66 916.87 889.21 
26.5 109 .60 079 .99 101.55 434.55 913.58 884 .87 
27.5 106.90 076.22 098.60 431.29 910.11 880 .34 
28.5 104.11 072.34 095.51 427.89 906 .57 875.70 
29.5 101.23 068 .32 092.29 $24.36 902.85 70.93 
30.5 098 . 23 064.19 088 .95 420.67 898 .97 866.00 
31.5 095.09 059.98 085.48 416.86 894.94 860.92 
32.5 091.78 055.66 081.87 412.91 890.78 855.71 
33.5 O88 .42 051.17 078.15 108.81 886.49 850 .32 
34.5 084 .97 046.59 074.30 404.58 882.03 844.85 
35.5 O81 .32 041.88 070.31 400.21 877.44 
36.5 077 .62 037 .05 066.19 395.71 872.71 
37.5 073.79 032.15 061.96 391.06 867.81 
38.5 069 .83 027 .12 057 .60 386.27 862.77 
39.5 065.73 021.95 053 .09 381.34 857 .60 
40.5 O61 .56 016.69 048.48 376.29 852.2 808.76 
41.5 057 .27 011.28 043.70 20998 . 60 371.08 846.81 802.25 
42.5 052 .83 005.76 038.86 992.69 365.74 841.20 795.57 
43.5 048.29 000 .16 033 .86 986.61 360. 26 835.43 788.76 
44.5 043 .64 19994 .38 028.73 354.65 $29.50 781.77 
45.5 038.86 988 .56 023.46 348.90 823.45 774.09 
46.5 033.98 982.58 018.09 343.01 817.26 767.45 
47.5 028 .97 976.49 012.58 336.97 810.91 760.04 
48.5 023.84 970.24 006 .93 3 804.42 752 .52 
49.5 018.57 963 .97 001.15 324.51 270 .92* 797 .79 744.82 
50.5 013.26 957 .50 20995 . 26 318.05 263.41* 791.00 737 00 
51.5 007 .77 950.99 989 .26 311.47 784.05 729.02 
52.5 002.24 944.30 983 .09 304.73 776.97 720.90 
5 19996.47 937 976.83 297 .90 769.74 712.63 
990 .62 970.41 290 762.38 704 .23 
6 984.70 963 .89 283.75 754.85 695.67 
5 978.63 916.45 957.21 276 .49* 7A7 686.95 
5 972 .46 909 .16 950.44 269 .05* 739.37 678.12 
5 966.18 901.86 943 .52 261 731.40 669.12 
959.74 894.42 936.48 723 .26 659 97 
5 953 .22 886.60 929 .32 715.01 650.66 
5 946.58 879.01 922.00 706.59 641.23 
62.5 939.81 371.21 914.56 698.04 631.63 
63.5 932.95 863 .26 907 .04 689.34 621.90 
925.95 855.24 899 .36 680.47 612.03 
45 918.85 847 .02 891.55 671.45 601.99 
66.5 911.70 838.68 883 .49 662.30 191.83 
67.5 904.21 830 .37* 875.50 653 .00 581.50 
68.5 896.81 822 .00* 867 .40 643 570.99 
69.5 889 .25 859 .03 784 633.90 560.41 
70.5 881.56 850.57 624.17 549 67 
71.5 873.83 841.94 614.27 538.77 
72.5 865.92 833 .29 604 .22 527.64 
73.5 857 .66 824.41 594.00 516.50 
74.5 849.55 815.45 583 .63 505 .08 
75.5 841.56 806 .22 193.56 
76.5 832.75 797 .14* 481.97 
77.5 824.57* 786 .37* 170.34 
78.5 815.43* 540.72 458.75 
79.5 529.48 446.16 
80.5 518.28 


*Denotes overlapped line. 
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TABLE (Concluded) 


o 


* Denotes overlapped line. 


387 .90 
381.78* 
375 .82* 


.81* 

3. 63* 
.25 
350. 83* 


344.53* 
337 .95* 


628.03 
621.31 
614.51 
607 .60 
600.60 
593.46 
586.28 
578.96 
571.54 
564 .06* 
556 .40* 
548 .48* 
540. 85* 
532 .83* 
524.79* 
516 .47* 
508 23* 
499 .74* 
491 .49* 
482 .85* 
474.16* 
465 .44* 
456 .36* 


(2,9) (0,6) (0,4) (0,3) 
R(J) P(J) R(J) P(J) R(J) P(J) R(J) P(J) 
5 
5 
5 18830 .38 
5 829 .09 
16727 .65 18836 .69 827.77 
.5 16736 .02 726.36 836.39 826.32 19472 .37* 
9.5 735.76 724.97 17570.51 17559 .21 835.74 824.79 470 .75* 
735.37 723 .52 569.72 557 .67 835.14 823.24 469 .02* 
5 734.91 721.97 569.18 556.06 834.54 $21.35 467 .18* 
734.31 720.31 568.55 554.31 833.80 819.48 55. 24* 
3.5 733.70 718.60 567 .82 552.39 832.87 817 .5: 463 .20* 
.6 732.95 716.74 566.94 550.52 831.85 815.35 461.01* 
5 732.07 714.86 566.17 548.45 830.78 813.24 
5.5 731.23 712.88 565.09 546.38 829.59 810.85 
5 730.18 710.80 564.06 544.16 828.25 808.51 
5 729.11 708 .62 562.73 541.86 826.83 806 .03 
9.5 727 .92 706 .39 561.45 539 .42 825.34 803 .38 
.5 726.68 704.05 560.05 536.94 823 .69 800.64 
725.32 701.62 558.56 534.38 821.97 797 2. 
699 .10 556.97 531.73 820.16 794.90 4139.80 
5 696 .47 555.28 528.96 818.20 791.85 436.63 
5 693 .82 553 .56 526.05 816.14 788.73 33.35 
5 691.06 551.67 523 .06 814.02 785.47 
5 688.18 549.71. 520.06 811.77 782.09 
685 .23 5A7 67 516.90 809 .39 778.65 
5 682.16 545.47 513.66 806 .94 775.08 
9.5 679 .05 543 .23 510.31 804.33 771.4: 
.o 540.98 506.88 801.64 767 .66 
5 5 503 .35 798.84 763.74 
499.71 795.96 759.79 
5 496 .00 792.97 755 . 398.89 
4.5 699 .79 492.20 789 .86 751 394.50 
5.5 697.14 188.30 786.63 747. 389.99 
694.44 484.28 783.32 742.7 
5 691.76 ; 480.16 779 .86 738.29 
688 .89 646.95 476.03 776.36 733.61 
9.5 685.96 643.01 515.55 471.7 772.74 728.97 
682 .92 638.89 512.29 467 .37 768.99 724. 
679.81 634.71 508 .83 462.89 765.12 719. 
2.5 676 .67 630.42 505.40 158.33 761.20 714. 
3.5 673.28 626.10 501.79 153.68 757 13 708 .97 
5 670.01 621.68* 498.13 448.94 752.97 3.6 
666.53 617 .19* 494.41 444.09 748.69 3 
5 662.90 612. 66* 490.54 $39.13 744.33 : 
659 .41* 607 .84* 486.60 434.11 739.81 
.5 655 .60* 602 .97* 482.60 428 .97 735.23 
651.65* 598 .11* 478.40 423.78 730.57 
647 94 593 .17* 474.18 418.37 725.75 
5 ti44 .00* 588 .22* 469.91 413.06 720.85 
5 583 .00* 35.46 407 .47 715.83 
> 35. 577 .75* 401.96 710.67 5 
5 631 .44* 396.26 705.52 
622 .74* 384.64 694.71 132. 
5 378.66 689.21 125.9 
5 613 .48* 372.60 683.56 619.2 
5 609 .21* 366.44 677 .82 612.38 
604. 26* 360.18 671.96 605. 
5 2 353 83 665.98 598. - 
5 416.08 347 47 659.91 591.3% 226.74 
589 .42* 410.638 340.98 653.75 584.07 219.16 
5 405.12 $34.40 647.46 576.74 211.45 
399.45 41.09 569.2 203.64 
393 .72 634.60 561.72 
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1 
1 
1 
1 
| 
2 
2 
4 
4 
4 
5 
Vv 
6 
6 
6 
67. 
68. 254.58 
70. 530.54* 240.10 
225.17 
201.93 
471.81* 
454.15* 
80.5 
81.5 
82.5 
83.5 
84.5 
85.5 
86.5 
87.5 
88.5 
89.5 
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TABLE 
Band heads and band origins for 


Coleman, Coleman, 
Gaydon, and Our Gaydon, and Our 
(v’, vo, Vaidya measurements Vaidya measurements 
(5,0) 4189.0 3.0 
(4,1) 2.4 
(3,1) 0.6 0.10 
(1,0) 
2,1) 0.5 0.01 
(3,2) 4624.9 0.0 
(4,3) 3.5 
(1,1) 1.3 
(2,2) 21133.11 4731.2 4730.29 1.0 0.01 
(3,3) 4763.3 
(0,1) 4806.1 0.4 0.02 
(1,2) 4844.5 1.4 
(3,4) —0.4 0.04 
(1,3) 0.4 
(2,5 5208.8 0.6 0.33 
5248.7 0.3 
(1,5 5345 
(2,6 5384.45 1.0 —0.01 
(3,7) 
(0,5) 5495 1.8 0.01 
(1,6) 5533 1.8 
(4,9) 3.4 
(0,6) 5692 2.3 0.01 
(1,7) 5730 2.5 
(2,8) 5767 —0.06 
(0,7) 5900 2.4 
(1,8 5939 4.0 
(2,9) 16734.46 5976 2.6 0.02 
(1,9) 6152.5 
(2,10) 6193 0.7 
(3,11) 6231.5 0.3 


*Co'eman, Gaydon, and Vaidya quote 4286. 2 Jk which appare ‘aie’ is aieas to a misprint. 
tOverlapped. 


The band origins were determined from the equation 


plotting against J*, using approximate 
values and extrapolating The values obtained are 
given Table More accurate values for the differences and 
These values were then combined with the and values determined 
Section obtain and values for the other vibrational levels. The 
rotational constants obtained this manner are given Table 
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TABLE III 
Rotational constants for 
Upper state 


Bands (measured) (calculated) 
(0,2) (0,3) (0,4) (0,6) 0.27491 0.33 0.27498 
(2,0) (2,2) (2,9) 0.26974 0.35 


Lower state 


Bands (measured) (calculated) 
2,0) (3,0) 0.33890 0.33891 
(0,2) (2,2) 0.36 0.33346 
(0,3) 0.33068 0.36 0.33070 
(0,4) 0.32794 0.36 0.32793 
(0,6) 0.36 0.3223: 
0.31377 0.37 0.31377 


(2,9) 


provide check the molecular constants derived Sections and 
the quantities AG, (B,,—B,,), and were determined from the 
appropriate combination relations between corresponding lines pairs 
bands having one state common (Herzberg 1950, 188). The evaluation 
these quantities was carried out two stages manner similar that 
discussed Section First, approximate values were assumed 
and the data used determine accurate values AG. Second, these values 
were combined with the experimental data determine more accurate 
values for (B,,—B,,) and The values AG, and 
(D,,—D,,) were found good agreement with those derived method 


The variation with for the upper and lower states shown Fig. 
the lower state marked curvature observed, but the variation very 
regular and was found that the values could fitted the quadratic 


equation 
(2) B, = 


within +0.00002 Final values for B,, and were determined 
the method least squares and are given Table IV. The values 
calculated from these constants are given Table IIT. 

For the upper state apparent that one more the values per- 
turbed. This perhaps not surprising view the predissociations observed 
near and Since the magnitudes the perturbations are not 
known, linear dependence was assumed and values for and 
were determined the method least squares. The constants are given 
Table and the calculated values for are given Table may 
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0.00273 
(cm-') 
0.2766 


0.2764 
0.2762 


0.2760;— 


Fic. Variation with for the ground and excited states 


noted that while the theoretical value calculated assuming Morse 
potential function (Herzberg 1950, 108) and our final molecular constants 
(Table IV), agrees well with the experimental values for the lower state 
factory for the upper state 0.00309 0.0027; 


TABLE 
Molecular constants for 


A 2T13/2 


Weve 


From the graphical procedures used Section was found that the 
values for the lower state showed very slight increase with increasing The 
magnitude this effect was comparable the accuracy the measurements, 
hence accurate value for could determined experimentally. Theoreti- 
cally, the value derived from our final molecular constants and the 


Oo 
0.3400 
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bining this value with the values given Table III obtain 
For the upper state the values cannot considered 
highly significant owing the perturbations present. The theoretical value 
for Combining this value with the values given 
Table III obtain (see Table IV). The experimental 
values and are good agreement with the values calculated from the 

Values for were calculated from the equation (A2) 
where the reduced mass the molecule atomic mass units 
126.9453, 16.0000). The values are given Table IV. 


The vibrational constants for the upper and lower states and the electronic 
term value for the excited state were determined fitting the band origins 
given Table the equation 


Only two vibrational constants could derived for the excited state owing 
the limited data available; for the lower state was found that cubic 


expansion fitted the data within the experimental error. The constants 
obtained the method least squares and are given Table IV. 


DISCUSSION 


From the rotational analysis the bands clear that the carrier must 
diatomic linear polyatomic molecule. Since the bands have been observed 
absorption during the flash photolysis mixtures iodine and oxygen, the 
carrier can only contain one both these elements. The only possibility, 
consistent with the rotational constants, and interesting note that 
the ground-state bond length obtained the basis this conclusion 
consistent with the ground-state bond lengths related molecules. Thus 
and 0.331 (Herzberg 1950). should pointed out 
that the alternative carrier proposed Vaidya (1937), viz. CI, cannot 
eliminated the basis the rotational analysis alone. This possibility 
excluded, however, the flash photolysis experiments and the observation 
Coleman, Gaydon, and Vaidya the bands emission from oxyhydrogen 
flame which iodine added. 

Our results confirm the vibrational assignments Coleman, Gaydon, and 
Vaidya. The wavelengths the band heads calculated from our final molecular 
constants are found satisfactory agreement with the measured values 
(see Table The vibrational constants the ground and excited states are 
now known with greater precision, but not possible improve significantly 
the earlier estimates for the dissociation energy since considerable extrapola- 
tions are needed for both states. The values quoted Coleman, Gaydon, and 
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Vaidya (1948) and Durie and Ramsay (1958) are kcal/mole and 
42+5 kcal/mole respectively. 

interesting note that several diffuse bands were observed emission 
from the flame source used even under conditions high-resolving power. This 
result marked contrast the general absence such bands emission 
from discharge tube sources. The appearance diffuse bands from flame 
sources may attributed the higher pressures prevalent flames and 
the fact that conditions hydrogen diffusion flames approximate thermal 
equilibrium (Durie 1952; Gaydon and Wolfhard 1953). 
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ANALYSIS SPECTROGRAM THE RED AURORA 
FEBRUARY 10/11, 1958, THE WAVELENGTH RANGE 


7300-8700 
VALLANCE JONES 


ABSTRACT 


high-dispersion photographic infrared spectrum covering the 7300-8700 
region was obtained from the great red aurora February 10/11, 1958. This 
spectrum shows the and 8446 lines have been enhanced least 
five times relative normal aurora. The most striking feature the spectrum 
was the appearance with high intensity the forbidden multiplet 
The measured wavelengths and relative intensities the features the 
spectrum are listed. 


INTRODUCTION 


The great type-A red aurora February 10/11, 1958, was observed simul- 
taneously over very wide geographical area. Spectra the have been 
obtained Alaska Belon and Clark (1959), Yerkes Observatory 
Wallace (1959), and near Moscow Mironov, Prokudina, and Shefov (1959). 
Saskatoon, spectrogram was obtained the region 
7300-8700 addition, set nine low-dispersion (330 spectra 
(over the range 3750-7000 were obtained with patrol 
spectrograph. common with the results obtained other investigators, 
these spectra show unusual enhancements some atomic lines oxygen and 
nitrogen. particular the infrared spectrum shows the forbidden 
doublet with very high relative intensity. The presence this feature 
was tentatively assumed Omholt (1957) low-dispersion spectra obtained 
from the higher-altitude portions auroral forms. Dufay (1959) obtained 
further evidence its appearance just resolving the doublet photo- 
electric spectrum low-latitude aurora. Wallace (1959) provided very good 
evidence for the appearance the doublet high-dispersion spectrum 
obtained Yerkes Observatory December 1958; his spectrum the 
lines appear superposed the almost coincident and lines the 
underlying airglow band. While the spectra Dufay and 
Wallace leave little doubt the identification the multiplet, the spectrum 
obtained Saskatoon puts the matter bevond question since the lines the 
doublet are completely separated and the night airglow lines are absent. 


EXPERIMENTAL 

The infrared auroral spectrum was obtained with dispersion 
the first order the auroral spectrograph Saskatoon. The 
spectral slit width was 5.4 Kodak spectroscopic plate, hypersensitized 

Contribution from the Department Physics, University Saskatchewan, Saskatoon, 
The research reported this paper has been sponsored the Geophysics 
Research Directorate the Air Force Cambridge Research Center, Air Research and Develop- 
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bath, was used. Wratten filter eliminated the 
second- and third-order spectra. The plate was exposed the red aurora 
February 10/11 for hours minutes, between 1820 and 0420 hours Mountain 
Standard Time. The plate remained the spectrograph and subsequently 
received further 5-hour 53-minute exposure aurora over period lasting 
until February 23. Some this later aurora was also type-A red aurora; from 
past experience with this spectrograph believed that the spectrum obtained 
almost entirely attributed the very intense display February 10/11. 

print the spectrogram reproduced Fig. together with night 


TABLE 


Relative peak intensity 


Normal Normal 


Observed aurora aurora 
Feature wavelength wavelength aurora (Meinel) (Saskatoon) 
(2F) 7319.7 7319.4 17.0 
13-3, 7330.1 7330.0 11.8 
7341.6 
7350.1 
7384.2 2.7 1.5 
7466.6 8.9 
7480.1 
Ne, (4-2) 
9.0 6.7 9.2 
No, (3-1) 7619.6 


7625.8 
(atmos.) 


7713.2 
7738.5 
7826.9 8.9 8.8 
Meinel, (2-0) 7860.8 8.2 8.6 
7875.9 9.7 
N3, 2-0? 7913.6 
(19) 7994.8 7995.1 3.1 1.3 2:6 
8054.7 4.6 4.6 4.6 
Meinel, (3-1) 8090.8 
NI, (2) (13-24) 8185.9 8186.6 1.4 
NI, (2) 8216.4 8216.3 4.2 1.9 2.8 
Unidentified feature 8241.9 
Unidentified 8411.8 3.0 
Unidentified 
OH, 6-2 8429.6 8430.7 


6 
Origin 
P-max. 


Wavelengths and relative intensities auroral features 7300-8700 


° 


— 
if 
a 


_ 
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Oll 
2F 19 


4-2 


Microphotometer tracing the February 10/11 red aurora spectrum. The positions 
the principal molecular bands and atomic multiplets are indicated. The atomic multiplets 
are identified the numbers adopted Moore (1945) and Chamberlain and Oliver 
(1953). 


airglow spectrum obtained with 45-hour exposure the same plate. For 
comparison, spectrum normal aurora obtained Saskatoon with 
spectrograph dispersion 132 A/mm and spectral slit width 9.3 repro- 
duced the same figure. The fourth spectrum Fig. that obtained from 
hollow-cathode discharge, with spectrograph still higher dispersion 
(20 A/mm). 

The wavelengths the features the auroral spectrogram were measured 
setting interpolation based the comparison spectrum super- 
posed the auroral spectrum. Unfortunately, the wavelengths the principal 
atomic features obtained this way showed systematic difference 1.4+0.6A 
from the accepted laboratory values. This was attributed the fact that the 
comparison spectrum was not applied until the plate was about removed 
from the spectrograph days after the principal auroral exposure; can only 
supposed that some slight mechanical shock thermal effect led the 
observed shift. The final wavelengths, corrected for this displacement, are 
listed Table The wavelengths should accurate within 

The relative peak intensities some the features listed Table were 
calculated with the help intensity calibration photographed the same 
plate. This calibration was produced exposing the spectrograph low- 
brightness source (Shepherd 1954), with rotating step sector the slit. The 
low-brightness source was calibrated against black-body standard. The 
relative intensities obtained are set out column Table 
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These relative intensities may compared with those listed Meinel 
(1951) whose values are listed column Table The spectral slit width 
employed Meinel seems have been about compared with 
about for the present study. Consequently, his spectrum, lines 
and sharp band heads may appear weaker factor 2.5 com- 
parison with unresolved bands and continuum. 

further comparison may made with the relative brightnesses listed 
column Table These values were obtained from normal auroral display 
observed from Saskatoon during the night April 5/6, 1957. The spectrum 
was photographed with 5-in. f/0.8 spectrograph (Petrie and Small 1952) 
plate. The plate was developed together with similar plate which 
intensity calibration had been applied described above. The spectral slit 
width for this exposure was about 9.3 For this spectrum, the relative 
intensities the features not differ significantly from the values reported 
Meinel. 


GENERAL FEATURES SPECTRUM 

comparison with the spectrum normal aurora the most striking points 
are: 

The and 8446 multiplets belonging the transition 
array are from times more intense relative the bands than 
normal low-level aurora. 

The OI, multiplet (19)* definitely present; its relative 
intensity may slightly enhanced comparison with normal aurora. 

The NI, multiplet (2) present with somewhat enhanced 
intensity relative the normal auroral spectrum. 

(iv) The relative intensities the first positive bands and the Meinel 
bands not seem much different. 

(v) The multiplet lines are present with extraordinary 
intensity. The measured wavelengths agree within 0.3 with the theoretical 
values. The intensity ratio between the lines was found 1.4. This 
fair agreement with the theoretical values Seaton and Osterbrock (1957), 
who predict the range values ratio 1.24 1.31. 

The atmospheric band appears particularly weak com- 
parison with the bands. This possibly due the increased dispersion 
the spectrograph. However, the spectral slit width the 9-in. spectrograph 
about half that the lower-dispersion spectrograph while the relative intensity 
the band about times less than for the normal auroral spectrum. 

The sharp maximum 7687.1 attributed Chamberlain, Fan, and 
Meinel (1954) the head band appears clearly. Strangely, 
its intensity relative that the neighboring 2-0 first positive band not 
very different from its value normal aurora; anything its relative intensity 
greater. This observation might appear raise some doubt the correct- 
ness the identification the 7687-A peak; however, must borne 


*The multiplet numbers are those adopted Moore (1945); see also Chamberlain and 
Oliver (1953). 
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mind that the sharp head the band would less affected the reduc- 
tion spectral slit width than the headless branches the band. Moreover, 
the rapid drop emulsion sensitivity beyond 8600 makes the intensity 
calibration liable considerable error the region the band. 

comparison between the auroral and the laboratory spectra Fig. shows 
that all the measured intensity maxima within the and bands correspond 
heads lines occurring the laboratory spectrum. The maximum 
7687.1 the auroral spectrum does not, course, appear the laboratory 
spectrum. 

There are few unidentified features listed Table The most intense 
the line-like emission 8411.8 This not night airglow line 
may seen comparing the auroral and night airglow spectra Fig. 


CONCLUSIONS 


The most striking feature the spectrum the February 10/11 red aurora 
this spectral region was the great enhancement the forbidden 
multiplet This observation may important one under- 
standing the mechanism type-A red aurora. The enhancement this 
multiplet very beautifully paralleled the enhancement the related 
multiplet observed Wallace for the same aurora. The 
intensification these transitions may very well height effect. 

The other notable peculiarity this spectrum the generally enhanced 
relative brightnesses its atomic lines comparison with those the mole- 
cular bands. This effect, which has been noted other observers for the same 
aurora, probably also height effect. 
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ROTATIONAL AND VIBRATIONAL INTENSITY DISTRIBUTION 
THE FIRST NEGATIVE BANDS 
SUNLIT AURORAL 


VALLANCE JONES AND HUNTEN 


ABSTRACT 


Spectra sunlit auroral rays were obtained from Saskatoon during the auroras 
September 3/4 and 4/5, 1958. The resolution these spectra was sufficiently 
high enable measurements made the relative intensities the lines 
the 0-0 first negative band well the relative intensities bands 
shows they are consistent with excitation process which ions thermal 
equilibrium with the atmosphere 2200° fluoresce under the influence 
solar radiation. The vibrational intensity distribution also consistent with 
fluorescent excitation from state thermal equilibrium about 2050° 
shown that the results are not consistent with fluorescent excitation process 
which the rotational and vibrational degrees freedom the ions come 
into radiative equilibrium with the solar radiation. Earlier conclusions that 
radiative equilibrium did hold for vibration are shown error result 
the high rotational temperature and the low dispersion used. concluded 
that the destruction ions asa result dissociative recombination proceeds 
sufficiently fast prevent any significant approach radiative equilibrium. This 
investigation provides strong indication that the kinetic temperature sunlit 
auroral ray (perhaps the region) the neighborhood 
2000° This may somewhat higher than the temperature the normal 


atmosphere this height. 


INTRODUCTION 


One the principal objectives the study aurora the determination 
the properties the upper atmosphere. The emission from the most common 
forms aurora concentrated between and 120 km. While the summits 
some auroral rays may occasionally reach high 600 the work 
Stérmer (1955) has shown that the forms occurring the greatest heights are 
the sunlit auroral rays. The highest measured points for these rays lie near 
1100 while the lower ends are generally found between 200 and 400 km. 
Consequently, spectroscopic observations sunlit aurora should enable 
information obtained about conditions the 400 600 height range. 

Spectra sunlit rays have been obtained Stérmer (1939). These spectra 
which were obtained from parts the rays heights between 400 and 650 
showed the first negative bands particularly strong and exhibit 
unusually great development the vibrational sequences. 

These results were interpreted Bates classic paper which 
was shown that the experimental data were consistent with the excitation 
sunlight resonance fluorescence process which special equili- 
brium distribution was set among the vibrational levels the ion. This 
distribution reflects the tendency the vibrational degree freedom come 
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into radiative equilibrium with the solar black-body radiation. For such 
steady state can arise for the case when collisions are unimportant since pure 
vibrational and rotational transitions are forbidden. Bates also calculated the 
population rates for fluorescent excitation the upper state for the case 
where the collision rate high enough maintain thermal equilibrium the 
ground state. The results the latter calculation did not seem fit 
data unless atmospheric temperature 5000° was assumed; this tempera- 
ture appeared too high even for height 650 km. However, view 
the high rotational temperatures reported for sunlit rays the present work, 
seems that data require reinterpretation. 

Swings (1949) pointed out that the rotational structure sunlit aurora 
should exhibit irregularities the same kind were found the spectra 
comets Swings (1941) and McKellar (1942). These irregularities were shown 
Swings and arise from the presence Fraunhofer lines the 
exciting solar continuum. The observation such irregularities the spectra 
sunlit aurora would provide definite evidence for the operation the 
fluorescent mechanism the excitation sunlit aurora. 

the nights September 3/4 and 4/5, 1958, high- and medium-dispersion 
spectra were obtained Saskatoon sunlit auroral rays. The analysis 
these spectra, which this paper devoted, shows definitely that the sunlit 
rays were excited resonance fluorescence and provides strong evidence 
support the view that the rotational degree freedom the ion remained 
thermal equilibrium with the atmosphere temperature approximately 
2200° The observed vibrational intensities the sequence 
bands also appear close that which would result from fluorescent 
excitation ions the distribution for thermal equilibrium approximately 
2050° This differs significantly from the distribution which would result 
from equilibrium with the solar radiation, which corresponds approximately 


OBSERVATIONS 


Displays aurora exhibiting sunlit auroral rays were observed from 
Saskatoon September 3/4 and 4/5, 1958. These rays showed the greyish- 
blue color described (1955) and times exhibited the characteristic 
intensity minimum where the rays crossed the boundary between the sunlit 
and dark regions the atmosphere. The spectrographs were directed towards 
the sunlit, upper parts rays. The sunlit portions the rays were recognized, 
general, their unusual color. The characteristics the spectra obtained 
show that this procedure was successful. 

Three separate plates were obtained from the sunlit aurora. Two these 
plates (one each evening) were obtained using the 5-in. spectrograph 
described Petrie and Small (1952). This instrument has linear dispersion 
A/mm. These two exposures were made Kodak plates and 
covered the regions and respectively. third, higher- 
dispersion spectrum was obtained using 9-in. spectrograph having 
dispersion 20A/mm. This exposure covered the sunlit parts the display 
both September 3/4 and 4/5; the spectral region covered was 
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reproduction the 3914-A band from the high-dispersion plate 
shown Fig. reproduction the same band obtained from normal 
aurora also shown. Microphotometer tracings the band have been pub- 
lished preliminary report our results (Hunten, Koenig, and Vallance 
Jones 1959). reproduction one the lower-dispersion plates shown 
Fig. together with spectrum aurora obtained after dark the same night 
(September 4/5, 1958). These two spectra were obtained using step slit 
having segments with widths ranging from 0.2 2.0 geometrical 
progression. Such slit provides means determining the characteristic 
curve the emulsion from region where the spectrum continuous. 


RESULTS 


The rotational structure the band for the sunlit auroral spectrum 
Fig. clearly corresponds much higher rotational temperature than does 
the normal spectrum. The relative intensities the lines the branch were 
derived from the microphotometer tracing with the help intensity 
calibration applied another part the auroral plate; this calibration was 
made photographing the spectrum standard low-brightness source 
(Shepherd 1954) with rotating step sector the slit the spectrograph. The 
relative intensities the branch lines obtained this way must then 
corrected subtracting the contributions the higher lines the branch 
which coincide almost exactly with the branch lines. The correction was 
made assuming that the branch line would have the same intensity 
the corresponding branch line from the same level. The results are given 
Table 

first approximation the corrected branch line intensities were then 
used derive rotational temperature according the usual graphical method 
(Section plotting log against K’(K’+1). The value for the 
state was used the relation between the slope the plot and the 
rotational temperature. The temperature obtained was 2100° This result 
must regarded preliminary one only since account has been taken 
the nature the excitation mechanism. more correct treatment given 
Section 

clear from Fig. that irregular fluctuations are present the line intensi- 
ties the branch; such irregularities are expected for fluorescent 
excitation the band system sunlight. 

was found Stérmer (1939) the relative intensities the vibrational 
bands the first negative system are also abnormal the sunlit auroral 
spectrum. This particularly evident for the 4278-A (Av —1) sequence 
bands may seen comparing the two spectra reproduced Fig. 
Figure shows intensity plot this region derived from the spectrum 
Fig. This intensity plot was obtained from microphotometer tracing 
across the fourth step the spectrum (counting from the strongest step). 
characteristic curve for the plate was derived scanning across the spectrum 
3775 and 4260 where the auroral spectrum quasi-continuous. Near 
the edge its field the spectrograph suffers from severe vignetting which was 


Normal 


Aurora 


Aurora 


The high-dispersion spectrum sunlit aurora compared with that normal 
aurora later the night; the rotational temperature the former much The 
second positive bands are abnormally weak the sunlit they are not 
enhanced resonance. 


PLATE 


SUNLIT 
AURORA 


NORMAL 


AURORA 


SEPT. 1958 


Neg. 


ow-dispersion spectra sunlit and normal aurora, taken with step slit. 


The enhancement the negative system with respect the rest the spectrum again 
noticeable. 


4 0 A 


VALLANCE JONES AND HUNTEN: FIRST NEGATIVE N; BANDS 461 


taken into account calibrating the sensitivity the instrument from 
4200-4300 with the help low-brightness source. Figure also shows 
plot the band sequence obtained from the non-sunlit auroral spectrum 
Fig. order keep the peak density the branch the 0-1 band 
within the range accurate microphotometry the intensity plot for the non- 
sunlit aurora was derived from the seventh step the spectrum. The conse- 
quent differences between the spectral slit widths for the two traces Fig. 


= 
100 
4 


Relative Intensit 
oO 


b 


Wavelength 


Fic. (a) Intensity tracing the 4278-A group the sunlit spectrum Fig. The 
dots show the synthetic spectrum for the band rotational temperature 
similar tracing for the normal auroral spectrum Fig. 


will not seriously affect the measurements the intensity ratios the bands 
since the shapes the bands are very similar. The vibrational sequence appears 
very much more strongly developed for the sunlit auroral spectrum. The 
significance this result considered further the next section. 


VIBRATIONAL INTENSITY DISTRIBUTION 

Bates has shown that the relative intensities the 4278-A (0-1), 
4236-A (1-2) and 4200-A (2-3) first negative bands are fundamental impor- 
tance the theoretical interpretation the spectrum. The relative intensities 
required, however, cannot taken directly from the curve Fig. because 
the high rotational temperature leads overlapping the bands. This 
not very important for normal aurora but major significance the case 
the sunlit aurora which has been shown above have rotational structure 
corresponding temperature over 2000° similar difficulty more 
acute form was encountered Clark and Belon (1959) the 
somewhat lower-dispersion spectra the sequence the bands. 

The correction the intensity the band was made calculating 
synthetic spectrum for the band. The wavelengths the rotational lines 
were taken from the data Childs (1982). value 2200° was adopted 
approximation the rotational temperature. The synthetic spectrum was 
then derived for spectral slit width summing the contributions 
all lines within range each wavelength which the spectrum was 
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required. The synthetic spectrum was calculated 4-A intervals for the range 
4236 The 4-A spectral slit width corresponds the value for the 
segment the slit which produced the spectrum from which the intensity 
plot Fig. was derived. This procedure supposes that the resolution function 
the spectrograph was rectangular; fact for this slit width was found 
distinctly trapezoidal (as determined from the profile the 4047-A 
mercury line). consequence expected that the synthetic spectrum 
would tend exhibit somewhat structure than the observed 
spectrum. 

The synthetic spectrum obtained shown superimposed the experi- 
mental curve Fig. The two curves fit well within the limits experimental 
error set the grain the emulsion. 

The relative intensities the 0-1 and bands may now calculated 
the assumption that the ratio will equal that the peak intensities the 
branches after the contribution the branch the 0-1 band has been 
subtracted. The result obtained 100:26. The result which 
would have been obtained correction had been made for the branch 
contribution the peak the band 100:39. The ratio derived from the 
intensity trace for the non-sunlit aurora (where the correction negligible) 
100:11. This agrees with the usual value the ratio for normal aurora (Hunten 
1955). For the sunlit aurora, the relative intensity the third band (2-3) was 
estimated extrapolating the background under its head about 10; the full 
set intensity ratios thus 100:26:10. 

discuss the excitation mechanism need the relative population rates 
rather than the The relation between them well known 
(Bates and need not repeated. longer necessary use 
factors instead relative transition probabilities the 
calculation since table the latter has been given Wallace and Nicholls 
(1955). The factors multiplying the relative intensities convert them 
relative population rates are then 0.93 for 1-2 and 1.06 for 2-3. thus 
find the third very uncertain. Reference Bates’ 
Table (1949a) gives vibrational about 2050° This 
not far from the rotational temperature 2200° found 
Section and certainly lower than the 4200° given Bates for the case 
radiative equilibrium with sunlight. 


THEORY THE ROTATIONAL ENERGY DISTRIBUTION 


The work Bates the vibrational intensity distribution showed 
that could apparently explained assuming that the molecules are 
equilibrium with the solar radiation. were not for the Fraunhofer lines 
the sunlight, the molecules would have vibrational temperature equal the 
temperature the sun’s surface; Bates’ result corresponds vibrational 
temperature about 4200° suggesting that the Fraunhofer lines have 
appreciable effect. Since the molecule homonuclear, has rotation— 
vibration nor pure rotational spectra; thus the only transitions governing the 
equilibrium are electronic. This effect illustrated results 
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(1942) comet spectra, which the molecule gives very high tempera- 
ture while the gives temperature about equal that the cometary 
gas (Swings 1943). 

The new observations sunlit auroral rays high dispersion show what 
appears the analogous effect the rotational energy distribution; 
although the temperature lower, about 2100° this might explained 
the presence Fraunhofer lines the exciting sunlight. The spectrum 
expected this mechanism was calculated and found indeed fair 
agreement with the observations. However, was thought best calculate 
another spectrum the assumption that the ground state had Boltzmann 
distribution corresponding the temperature 2100°; our surprise, this 
gave much better agreement. The methods used will now described; the two 

complete calculation for radiative equilibrium should take into account 
the absorption and emission all the bands the first negative and other 
systems, allowing for the intensity the sunlight each individual line. This 
would yield new version Bates’ result for the vibrational distribution, 
well the rotational distribution each rotational level ground and 
excited states. However, not practical for several reasons: far too much 
work would required; probably the wavelengths all the lines have not 
been measured; and the theory would almost impossibly complicated. 
Consideration was therefore restricted the first negative system and the 
first two vibrational levels each state; the band (3914 the most 
important, and the 0-1 band (4278 next. The 1-1 can ignored because 
the transition probability small, and the because the solar continuum 
much weaker the ultraviolet (3582 A); moreover, these bands the level 
which does not directly concern us. Only two bands were thus 
the calculated spectrum may expected exaggerate the effects individual 
Fraunhofer lines and allowance must made for this the comparison with 
observation. 

Let first consider the simple case which only the levels and 
are included. Figure shows the rotational levels numbered with 
the quantum number which convenient use when the spin splitting 
each level can neglected. Some the possible transitions are indicated; 
will seen that there are two independent sets levels since neither the 
absorbing nor the emitting transitions join them. The odd levels the one 
state are equilibrium with the even levels the other. The equations 
describing the equilibrium may set with the aid the principle 
detailed balance: the probability transition and its inverse proportional 
the statistical weight the final level. use and represent transition 
probabilities for emission and absorption, for the population state, and 
for the flux radiation the wavelength line. Then 


and 


y 

) 
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The first these expresses the fact that equilibrium the number transi- 
tions per unit time must the same for both absorption and emission between 
upper level numbered and lower level numbered K;; this transi- 
tion the branch and the flux called The second equation expresses 
the same fact for the transition having the same lower level. Dividing the 
second equation the first gives 

The quantities and represent radiation fluxes multiplied the statistical 
weight the upper level concerned; their use will allow later equations 
written more compactly. 


(a) Rotational levels and transitions for part the band. 
same, but including both and bands. (The vibrational numbers 
the state should and reading upwards.) 


When include the level the situation becomes much more com- 
plicated; shown Fig. Because the assumption equilibrium 
possible consider the transitions between the dotted lines themselves. 
But within these lines there now loop, that equilibrium 
without having each absorption transition individually equal rate the 
corresponding emission. Instead, must equate the number transitions 
reaching each level the number leaving; the four equations thus found, 
one redundant. Quantities relating the 0-9 band will given the sub- 
script and will used for the 0-1 band. The relative transition probabilities 
these bands will called and The result will first for 
the specific levels shown Fig. 4(6) and then generalized. 

constant, 


statistical weights appear the left because they are included the 
quantities and For constant, 


for constant, 
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the first two equations are substituted the third eliminate the 
populations the lower levels, the result 


more convenient form obtained dividing top and bottom 
using s(A) defined (1) and writing for the same time may 
replace the particular values general ones: 


interesting consider the special case which the illuminating flux has 
black-body spectrum; the ratio depends the wave-number difference 
the and lines with the same and this just the difference 
between the energy levels K+1 and K—1. From the definition s(K) (1) 
follows that 


thus this case the rotational temperature equal the temperature the 
black body would expected from more general reasoning. The solar 
spectrum the blue and ultraviolet deviates considerably from that 
black body; the detailed calculations show that the Fraunhofer lines happen 
lie such way simulate much lower temperature near 1200° 

Equation (2) not suitable for calculation for two reasons: first, its com- 
plexity, and second, because gives only the ratio adjacent odd even 
levels and large error could accumulate over the whole width the band 
unless unreasonable number extra digits were carried. (2) examined 
seen produce weighted mean and with the quantities 
square brackets acting weights. The value the mean will not greatly 
changed all the weights are taken equal then may write 
approximately 


The numbers follow from the values and given Wallace 
and Nicholls (1955) and the definition just before equation (2). The 
radiation fluxes and may found fractions the continuum from the 
Utrecht Atlas (Minnaert a/. 1940) and will call their ratio U(K); then 


i 


where have again used the Wien approximation the black-body law. 
When these are substituted into (3) the Boltzmann factors are common; and 
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when whole series ratios given (3) multiplied together, all the 
Boltzmann factors combine: 


(even levels). 


similar expression holds for the ratio odd levels Build-up 
error has now been eliminated from the Boltzmann factor but not from the 
term involving thus, one must careful estimating the U’s from the 
Atlas and making the calculations. can taken equal 
The case much simpler. The question what the 
molecules are thermal equilibrium with will considered later; for the 
present, the term will taken imply that the ground state has 
Boltzmann distribution the rotational levels, instead distribution deter- 
mined the radiative equilibrium just discussed. This distribution will 


longer represents the temperature the sun, but rather that the 
ground-state distribution. absorption the band alone considered, the 
distribution the upper state 


This will altered very little both the are replaced what nearly 
their average, then 


Throughout this development, normalizing constants have been omitted since 
only the distribution among levels interest. (6), and can taken 
the fractional intensities and from the Utrecht Atlas with little error: 
the variation the continuum cancels first order since the and branches 
opposite directions. second order, however, there appreciable 
variation the mean wavelength corresponding and lines; this what 
causes the head the branch and the increasing spacing the lines the 
branch. approximate calculation this effect showed that the measured 
temperature would about 50° too low 2100° and the color 
temperature 7150° This small correction hardly significant view the 
other errors measurement and interpretation, but included our final 
result; was not included our preliminary announcement (Hunten, Koenig, 
and Vallance Jones 1959). 

The band can included adding the last term (6) the amount 
where the primed fluxes refer this band. The fact that the first 
vibrational level less populated than the zeroth balanced the greater 
solar flux the longer wavelength; thus, only the ratio the transition 


a 
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probabilities enters long the vibration the molecule radiative 
equilibrium. Actually, the vibration somewhat less than this, but use the 
full factor will help compensate the omission other bands. Replacing the 
suggested above, finally get from (6) 


the next section the abbreviation used for the quantity the end (7). 


NUMERICAL CALCULATIONS AND RESULTS 


The biggest part the labor calculation the reading the fractional 
intensity sunlight the position each line the and 0-1 bands. 
already mentioned, the solar spectrum used was that the Utrecht Atlas, 
which gives intensity tracings large scale for the light from the center 
the sun’s disk. would better have spectrum integrated sunlight, 
but this much more difficult obtain and not available. The finite resolu- 
tion the spectrograph used obtain the Atlas fills the narrow lines 
however, the lines for the whole sun must somewhat filled 
also because the rotation and the resulting Doppler shift; thus, these two 
errors tend cancel. The wavelengths the lines are given Childs 
(1932); they were marked lightly the tracings the Atlas and the fractional 
intensities estimated closely possible. Because the Fraunhofer lines are 
often narrow and the rotational lines may fall their steep sides, some 
these results are rather uncertain, perhaps even sometimes (Those 
who have Utrecht Atlas available may find reproduction one its pages 
just this wavelength region the paper McKellar (1942).) The radial 
component the earth’s orbital motion and its rotation cause Doppler shift 
the whole spectrum which may have appreciable further effect; was 
not taken into account and thus provides another source error. With comet- 
ary spectra considerable change radial velocity required affect the 
spectrum seriously; thus, this error should small. Each branch each band 
was followed little beyond 40; for those lines with measured spin 
splitting both intensities were read and then averaged. 

The tracings the Utrecht Atlas were supposed adjusted that the 
continuum was the 100% level. has since been discovered Chalonge 
and co-workers that the near ultraviolet the continuum well above 100% 
places, primarily because the wings Balmer lines (Canavaggia and 
Chalonge 1946; Canavaggia, Chalonge, Egger-Moreau, and 
1950; brief summary Minnaert 1953). The best estimate the corrections 
Michard (1950) are follows: the continuum taken having color 
temperature 7150° 114% from 3926 3892 and 124% from 
3893 3793 The discontinuity occurs the middle the band and 
must taken into account, especially the case radiative equilibrium. 

Calculations for radiative equilibrium were made means (4) except 
that the weights used were 0.70 and 0.30 because numerical error. was 
not thought worth while repeat the calculation, especially since the result 
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neglects all the other bands the system and slightly too large weight for 
the 0-1 band would help compensate for this. B’’ was taken 1.922 
(Childs 1932; Herzberg 1950). Separate calculations were made for odd and 
even levels and the results are shown separately Fig. means lines 
joining the plotted points; the observations are shown the dots. The two 
sets could combined assuming them populated according the 
ratio their statistical weights, but there did not seem any advantage 
this. Since the ordinate scale logarithmic, the various sets points can 
moved and down without change shape; this procedure was used fit 
the corresponding ones together well possible. The vertical displacement 
between the odd and even lines arbitrary. While there agreement some 
places, the general fit between theory and observation poor. 

Equation (7) was used calculate the results for thermal equilibrium. The 
temperature assumed was 2100° found Section from the observations; 
another calculation using 1700° was also made, but this temperature was 
obviously too low. The weights used for the two bands were the ones (4); 
this result given Fig. along with the observations. The 2:1 intensity 
alternation has been removed doubling the measured intensities the weak 
lines; this gives smoother plot which easier interpret. The agreement 
between calculation and observation considered good and show that 
the distribution rotational levels the ground state close Boltzmann 


INTENSITY 


LOG 


(a) The synthetic spectrum the 0-0 band for thermal equilibrium, indicated 

lines joining the calculated points. The observations are shown the dots. Since the 

ordinate scale logarithmic, the spectrum can displaced vertically any desired amount. 
observed intensities the weak lines have been doubled. 

(6) Similar plots for radiative equilibrium. arbitrary displacement has been 

used separate the lines odd and even 


(a) 
e 
e 
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one. The biggest disagreements are found the ends the branch where the 
lines are faintest and the measurements least certain; addition, the first few 
lines ride above intense series lines which must subtracted. The 
possible errors the calculations have already been discussed; the principal 
ones are the errors the fractional intensities the sunlight and the omission 
the weaker bands the system. 


TABLE 


Observed intensities the lines the sunlit spectrum, corrected for the presence the 
branch. Also included are the calculated Fraunhofer correction factors 
defined the text after equation (7) 


Observed Observed 
169 199 676 188 
163 108 496 153 
533 126 254 149 
574 132 200 178 
826 201 241 180 
195 110 234 215 
(1-1 head) 214 


Once has been established that excitation from ground state with 
Boltzmann energy distribution, more accurate temperature for this distribu- 
tion can found modifying the standard plot A’) vs. 
take account the Fraunhofer lines. well known (Herzberg 1950), 
this follows from equation similar (5); replaced take 
into account slight dependence transition probability A’. Since the 
excitation now described (7), the ordinate must changed 
where the average residual intensity defined after (7). Such plot shown 
Fig. along with what thought the most probable straight line; 
gives temperature 2150° Adding the correction 50° for the slope 
the solar continuum, find 2200° the final result; could probably 
error +200°. The line was drawn favor the points from 
and 35; the points the ends are uncertain for the reasons already dis- 
cussed, and the ones between and appear raised the presence 
the band. The values used are given Table along with the measured 


line intensities. 


| 
‘ 
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(log scale) 


500 1000 


Fic. log against for finding the rotational temperature. Inten- 
sities the weak lines have been doubled. 


For later discussion will convenient have estimate the rotational 
corresponding the distribution for radiative equilibrium. The 
representation temperature not very good, but plot the synthetic 
spectrum according (7) gave 1200° for the strongest part the band. This 
much lower than the solar temperature because number 
strong, wide absorption lines the region the branch the band; 
prominent among these //-zeta. 


DISCUSSION 


this section wish discuss the relation between the rotational and 
vibrational temperatures the ions and the kinetic temperature the 
atmosphere. Since the term used somewhat specialized 
sense, will best review this first. may assumed that the mole- 
cules are thermal equilibrium, that strictly correct speak 
rotational and vibrational temperatures which are equal the kinetic tempera- 
ture But observe ions, not neutral molecules, and the observed 
rotational and vibrational may differ from such case 
not strictly correct use the term temperature, but convenient 
use parameter which specifies the distribution among the states. 
will use the term represent the rotational and the 
vibrational the ground state the ion The observed 
intensities actually give the distributions the excited state but since 
the excitation mechanism known the ground-state distributions 
deduced. This was done for the rotation Section and for the vibration 
Section 


° 
° 
2 o 
° 
° 
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The processes that control and will now considered. First have 
the original ionization process 


(8) (fast) 


will usually electron even the primary particles are fast protons 
because the large number fast secondaries (Omholt 1959; Bates, 
McDowell, and Omholt 1957). The measured cross sections indicate that most 
the ions are produced their ground electronic state (Omholt 1959). 
factors for this process are given Bates (1949a) and from 
guarantee that the vibrational distribution resulting from excitation ioniza- 
tion can represented temperature, but this case good representation 
possible.) The result shown Fig. For low about 1400° 


6000 
=x 
° 
Ww 
N 
tu 
z 
2000 
> 


Effective vibrational temperature the ground state the ion resulting 
from electron impact temperature 


reflecting the transition probabilities from the lowest vibrational level; for the 
range 2000° which interests most, about 250° higher than 
first approximation, unchanged ionization because the 
restriction the selection rule and the small change internuclear distance. 
However, the actual effect easily found since the RT’s are the same ratio 

(8) slow proton the most probable reaction charge exchange 
which the electron after the collision bound Experiments bearing 
this process have been made most recently Fan (1956), Carleton (1997), 
and Roesler, Fan, and Chamberlain (1958). They seem show that there 
bombarding energy kev. The effect when the initial temperature 
high not known, but there would probably still some increase. Thus one 
would expect find greater than charge-exchange excitation were 
affecting appreciably, but would still the same 
significant difference was observed for either the sunlit the 


aurora, 


a 
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The life ion ended recombination 


(9) 


which has large coefficient; the value has been measured 
(Faire and Champion 1959). lower value may more appropriate for our 
purpose, since the molecules the laboratory afterglow may well vibration- 
ally excited; adopt cm*/sec which may much order 
magnitude too small. the presence aurora, the electron density probably 
rather higher than normal; assume 

The mean life ion thus only about seconds; smaller electron 
density will give longer life, but unlikely exceed 200 seconds and may 
short seconds. this short time and are not able depart 
very far from their original values; thus, the cycle (8) and (9) effective 
keeping the ions close thermal equilibrium and far the most effective 
process high altitudes. 

the ion scatters solar radiation, its vibrational and rotational distributions 
are altered; the final result for the latter discussed Sections and and 
for the former Bates have seen, the rotational distribution, 
for the case radiative equilibrium, can represented roughly 
1200° The corresponding may found comparing Bates’ 
Tables and 13, and about 4200° Both temperatures differ from the 
sun’s color temperature because the effect Fraunhofer lines. The speed 
which radiative equilibrium approached depends first all the rate 
scattering, which was calculated Bates Lytle and Hunten (1960) 
have revised the result agree with more recent measurements, and find 0.11 
sec for the whole negative system; thus, the mean time between 
scatterings seconds. the estimate the lifetime ion given above 
correct, can scatter only 5.5 quanta, the average, before recombining; 
even the longest reasonable life 200 seconds permits only quanta. 
Because the selection rule (AK +1) and the large number rotational 
levels, the rotational distribution can affected only slightly this time; 
neither restriction operates the which may rise rather rapidly towards 
4200° The observed value will time average and will lie between 
and 4200°, probably closer the lower limit. near 2200° 
appears be, may slighly less than its original value 0.96 

Let now consider the effect collisions with other gas molecules. 
known (Massey and Burhop 1952) that rotational energy easily 
but vibrational energy transfer essentially impossible the time scale 
are concerned with. For the atmosphere there are, however, two special 
processes which effectively cause rapid transfer both types energy: 


(10) (charge exchange), 
and 
(11) (atom exchange). 


Charge exchange creates new ion for which and are both equal 
The effect atom exchange not obvious, but seems likely cause 
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equipartition among rotational, vibrational, and translational degrees 
freedom; thus, was large before the reaction, will increased after 
it. This seems the only mechanism which can become greater 
than 

The coefficients for (10) and (11) may the order gas-kinetic, 
they are compete effectively with (8) and (9) the mean 
time between reactions must less than about seconds, that the density 
estimated with the aid Nicolet’s model the thermosphere (1959): 250 
300 for and not far above 100 for Rotational energy can 
transferred atomic oxygen, presumably again with about gas-kinetic cross 
section requiring before density atoms/cm? which found the 
400- 500-km region. 

Finally, the possibility direct excitation electrons should considered. 
Excitation accompanies only few per cent ionizations, but each ionization 
followed the scattering five six photons. Similarly, excitation the 
ion electron impact during its short life improbable. Photon excitation 
always much more important than electron excitation, with the latter serving 
only produce the ions that they may scatter light. 

conclude that above about 300 the only important reactions are (8) 
and (9); approximately equal and somewhat higher. Below 
300 km, exchange reactions become more important and and should 
not differ appreciably from Our aurora was probably well above 300 km, 
perhaps 400 500 km, and gave 2200° and about 2050° 
Thus appears that was not far below 2200° Whether this should 
taken the temperature the normal atmosphere another question, since 
considerable energy was being deposited reactions such (8) and the 
action solar radiation raising VT. That found less than 
instead greater not significant, since neither accurate better than 
+200°. 

Some accepted conclusions about sunlit and aurora must 
reviewed the light our results. Section was shown that cannot 
found from the spectrum without first making substantial correction for 
the overlapping bands which occurs when high. Since this high 
was not suspected previously, many erroneous have been derived. 

The measurements (1939, 1955) blue sunlit rays were used 
Bates (1949a) deduce that corresponded with the value expected from 
radiative equilibrium closely would expected. Study 
microphotometer tracing shows that the appearance the 4278-A group 
very similar that Fig. though much lower dispersion. gives 
the intensity ratio 100:59; this compared with our 
uncorrected result 100:39. highly probable that the true intensity ratio 
close our 100:26 and the close 2050° The lower resolution 
spectrum would expected raise the uncorrected intensity 
the band observed. 

curious type auroral spectrum has sometimes been observed from low 
latitudes (Rayleigh 1922; Barbier 1947), and more recently from the latitude 
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Moscow (Mironov, Prokudina, and Shefov 1959). Its characteristics are the 
very high apparent shown the first negative bands and the almost 
complete absence the second positive bands. Both these facts strongly 
suggest that the aurora sunlit, idea which means new (Seaton 
1956; notice particular his statement that Barbier’s aurora may have been 
sunlit despite earlier statements the contrary). Sunlit aurora should 
observable for long distances because its great height, and its spectrum 
observed from far away less likely contaminated light from ordinary 
aurora than when observed from below (Chamberlain and Meinel 1954). 
the other hand, seems likely more contaminated other high- 
altitude emissions than the spectrum which was produced guiding 
visible rays. All the spectra mentioned were produced 
long unguided exposures, though some cases the sunlit phase must have 
been rather short. would seem that rather than 
the important characteristic this aurora, although the 
latter name has the advantage established usage. 

All the spectra appear show even greater development 
the negative system than this particularly evident the 
4709 sequence which the 1-3 band brighter than the 0-2 and the 2-4 
bright. Unfortunately, our spectrum did not cover this region, but one the 
plates Mironov, Prokudina, and Shefov (1959, Fig. covers very well. 
Fortunately, the resolution this spectrum just good enough show some 
lines certain multiplets whose other lines must blended with the 
heads the two anomalously intense bands. rough allowance for these lines 
brings the relative intensities into reasonable agreement with Stérmer’s. The 
and bands are the most likely free blends; they appear very 
similar ours. Their Fig. shows band whose appearance strongly 
suggests high RT. 

suggested that some the exposure the atomic lines may have 
accumulated during the night before the visible aurora appeared, and that the 
latter then provided the sunlit bands. rough calculation the shadow- 
height during the visible aurora, using the times and pointing data given, 
shows that varied between 410 and 330 km. This consistent with the 
hypothesis that the display was sunlit. 


CONCLUSIONS 


The spectra taken the sunlit auroral rays occurring September 
and 1958, have been described. Even inspection was clear that the 
rotational and vibrational the first negative bands were 
unusually high. preliminary value the former was 2100° the latter was 
found about 2050° comparison with calculations Bates. Detailed 
calculations the rotational structure resulting from fluorescent excitation 
sunlight were made either two assumptions: first, that the rotational 
distribution was determined the fluorescence process 
and second, that the ground state the ion had rotational 
distribution the Boltzmann type equilibrium’’). The second case 
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gave better agreement with the observed spectrum. was then possible 
derive improved rotational temperature method which took account 
the solar Fraunhofer lines, and the result was 2200+200° 

thus appears that both the rotational and vibrational degrees freedom 
the ion were near thermal equilibrium with the atmosphere. discussion 
the possible mechanisms thermal relaxation showed that the most important 
high altitudes seems the short lifetime the ion against recombination, 
after which replaced new ion. The new ion rotational temperature 
slightly less than the kinetic, and vibrational slightly higher. 
photons are scattered the vibrational rises rather quickly 
but the rotation affected only very slowly. After scattering only few 
photons (as few five six reasonable assumptions) the ion recombines. 
The kinetic temperature the auroral rays cannot have been much below 
2200° but the rays may possibly have been hotter than the normal 
atmosphere. 

This work reverses the accepted conclusion that the vibration the ions 
radiative equilibrium sunlit auroral rays. Earlier spectra did not have 
high enough dispersion show. the large rotational temperature which causes 
neighboring bands overlap. appears that aurora requires 
the same correction, and also others for blending with atomic lines. With these 
corrections this type aurora gives results very similar those found this 


paper. 
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DAWN ENHANCEMENT AURORAL 


ABSTRACT 


The 0-0 band 3914A was observed once per minute with scanning 
spectrometer sunlight illuminated the upper atmosphere the morning 
September 1958. reasonable assume that the density ions was 
maintained constant during this time, and that the observed rise intensity 
was due resonance scattering sunlight. The ions were found distri- 
buted between 180 and with peak 250; density here was 
There was some evidence for another layer Rotational 
temperatures measured while the main layer was fully illuminated agreed well 
and gave average 1060+100° 


INTRODUCTION 


Measurement rotational temperature bands aurora has been 
carried out for many years. This experience has shown that difficult 
find measurable aurora much above 100 km, most temperatures being below 
400° Sunlit auroral rays are much higher 1955) and are bright 
enough give good spectrum the band (3914 short time. 
accompanying paper (Vallance Jones and Hunten 1960) gives details the 
measurements high sunlit rays September and 1958, with high- 
dispersion spectrograph. During the same display was possible observe 
continuously the dawn enhancement September with scanning spectro- 
meter spectrum/minute. There reason suppose that the density 
was maintained constant during the half-hour required for this enhancement; 
its vertical distribution can then found well-established methods used 
for twilight analysis. The ions were found occupy layer about 100 
thick; measurement intensity gave peak density 
250 km. The measured rotational temperature 1060° presumably applies 
this mean height and should essentially equal the kinetic temperature 
the atmosphere. 


OBSERVATIONS 

The spectrometer was used with resolution normal for rotational- 
temperature work with this instrument (Shepherd and Hunten 1955). 
usual make scans lasting seconds aurora, but this occasion the 
amplifier for the high-speed recorder was not working and was necessary 
use slower one Fortunately, the aurora was unusually 
quiet, especially during the dawn enhancement; thus, the slow scans were 
disadvantage and fact gave improved sensitivity. The typical spectrum 
shown Fig. illustrates the low dispersion and appreciable noise level; but 
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3875 A 


Fic. Sample spectrum taken shadow height 150 km. The estimated solar back- 
ground shown also. 


the results are consistent from scan scan and appear reliable. The region 
covered was from 3875 somewhat past the head 3914 while this 
normally includes the whole the band, did not this time because 
the high temperature which was not recognized until later; lines with 
greater than were missed. now obvious that the filling the zero gap 
the returning branch indicates the high temperature normal 
aurora shows the gap clearly the resolution used (Shepherd and Hunten 
1955). 

Towards morning the aurora was observed rather faint but with con- 
stant and fairly uniform intensity all over the sky. Scattered clouds were 
present but had observable effect the intensity, presumably because they 
were well illuminated auroral light. The sky was fairly clear the north; 
the spectrometer was left pointing this direction zenith angle 55°. 
Since the sun rises almost exactly the east this time year, the shadow 
height for observations the north the same the zenith. The intensity 
function shadow height shown Fig. with indication time 


200 300 
SHADOW HEIGHT (KM) 


Band intensity geometrical shadow height Mountain Standard 


(105th meridian) also shown. Each point average about five individual spectra. 
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well. The intensity remained kilorayleighs (kR) for least minutes, 
until the shadow reached the 280-km level; then steadily rose 
between 170 and 120 km. The further rise starting 120 seems real, 
but the measurements this region are unreliable because large continuous 
background had subtracted; was equal the band intensity km. 
The assumption that the enhancement was entirely due resonance seen 
reasonable, but cannot proved. Proof might have been obtained 
measuring some other emission which does not come from resonance transi- 
tion, but this was not feasible. 

Rotational temperatures were measured for the region 200 140 when 
the intensity was nearly constant. The unresolved profiles were corrected for 
convergence the lines described Shepherd and Hunten (1955). The 
plots vs. gave the expected straight line for between 
and 25; below the returning branch raises the intensity and this was 
observed. The effect illustrated Hunten, Koenig, and Vallance Jones 
(1959) for higher temperature obtained during the same display. The spectra 
were averaged four groups giving temperatures 1100, 1125, 1000, 
and 1025° the average these 1060° the shadow fell still more, 
the temperature seemed rise still higher, but this measurement even more 
unreliable than the intensity one low shadow-heights, because even small 
error estimating the continuous background could give large error the 
temperature. The measurements given for the greater shadow-heights are not 
appreciably affected this difficulty. 

DENSITY AND DISTRIBUTION 

find the vertical distribution ions must differentiate the curve 
Fig. this can only done roughly but least the heights which the 
derivative goes zero are reasonably well defined. The symmetrical curve 
chosen shown Fig. for any twilight investigation, also need 
calculate the transmission the lower atmosphere function height 
above and below the geometrical shadow, since this gives the form the shadow 
imposed the incident sunlight. The calculation was done essentially 
described Hunten (1954) using extinction 0.040 per atmos- 
phere (Allen 1955); this includes small contribution from ozone which was 
lumped with the rest since small. Dust the lowest atmosphere was 
represented coefficient 0.042 per atmosphere falling off with scale 
height 1.0 km. Numerical differencing the final table gave 
which call the geometrical shadow-height used Section 
plotted Fig. the center mass above the geometrical shadow 
and the second moment has recently been shown (Hunten 1959) 
that the brightness derivative plotted against km) Fig. may 
regarded the ion distribution broadened “slit From the 
relative breadths these two curves may seen that the broadening 
negligible considering the uncertainties the observations, and this con- 
firmed trying the sharpening procedure used for twilight sodium observations. 
thus gives the ion distribution directly. the lower laver really 
present, its height would about 120 km, common height 


aurora, 
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the transmission function plotted against the height above the geometrical shadow; this 
height shown the right. ion density (units 


find the total amount use the observation that resonance 
the same layer seen the zenith would fainter the van Rhijn factor 
0.616 calculated for height 250 km; this gives zenith brightness kR. 
According the Appendix, illuminated sunlight scatters 0.068 quan- 
tum/ion sec; thus have (38/.068) ions/cm? (column). The 
half-thickness the layer shown Fig. about km; thus, the average 
density about and the peak This gives the ion- 
density scale Fig. 


DISCUSSION 


should repeated that most this analysis rests the assumption that 
the ion density remained constant during the minutes taken the 
shadow sweep over the 200- 300-km region. While this plausible, there 
way being sure without other measurements which were not made. 
The height distribution found similar that high auroral rays, though 
not suggest that the observations were ray structures, rather, 
sort afterglow remaining from night unusually strong and high aurora. 
The high sunlit rays observed the evenings September and were 
presumably somewhat different phenomenon; this confirmed 
higher temperature about 2300° which suggests higher altitude 
around 500 km. 

Our rotational temperature does not depend the above 
assumption and seems reliable. This temperature reasonable one for 
height 250 km, and probably agrees with the kinetic temperature there 
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suggested the analysis Vallance Jones and Hunten (1960). The ions 
apparently occupied layer between 180 and 320 with peak 250. Their 
abundance was ions/cm? (column) and the peak density 
These numbers were deduced from brightness for 
the dawn enhancement observed zenith angle 55°. 

another opportunity comes observe high sunlit aurora with scanning 
spectrometer, the scan should possible extended 3840 include the 
lines higher rotation. Since this type aurora not very bright but 
rather quiet, scans can used and perhaps wider slits, 
instead preferred for normal aurora. The better resolution required 
only for the region near the origin which useless high-temperature aurora 
because the presence high lines the branch. 

all possible, twilight series the kind discussed here should 
accompanied measurements bright line band which shows 
resonance enhancement. This would allow check our assumption that the 
density remained constant during the twilight. bright band the 
second positive system would probably best. 
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APPENDIX 
SCATTERING POWER FOR THE BAND 

The fundamental calculation this quantity was made Bates (1949). 
The ions were taken having rotational temperature 273° and the 
exciting intensities each line were found from the Utrecht Atlas the solar 
spectrum. For the case which collisions could neglected, found 0.13 
quantum scattered per ion per second for the whole band system, and 
this the 3914-A band. oscillator strength 0.04 was calculated and 
rather larger value was calculated later Shull (1950, 1951), but 
recent measurement Dalby and Bennett (1959) agrees better with Bates’ 
result: Combining this with Bates’ figure given above, 
find scattering power for the band 0.068 quantum, ion sec. This not 
strictly correct for our much higher rotational temperature, but should not 


greatly error. 


ROTATIONAL SPECTRUM THE MOLECULAR 
BEAM ELECTRIC RESONANCE METHOD! 


GREEN? AND LEW 


ABSTRACT 


measured means the molecular beam electric resonance method. The follow- 
ing rotational constants have been determined (all frequencies 


=a, 


The quadrupole interaction constants measured the state are 
found 


v=0 —7.932+0.003 Mc/sec, 
—7.835+0.010 Mc/sec. 

The equilibrium internuclear distance obtained directly from 

The electric dipole moment the state 


INTRODUCTION 


When the rotational spectrum RbF was obtained this laboratory (Lew 
al. 1958), the remaining alkali halides for which the rotational constants 
were unknown were LiF, NaF, KF, and LiCl: these had all proved unsuitable 
for microwave absorption methods (Honig a/. 1954). Since was expected 
have the lowest rotational frequencies among these molecules, was the 
next one chosen for investigation. From the equilibrium internuclear distances 
17,400+175 Mc/sec. 

Potassium fluoride has been studied twice previously beam methods: 
each time observations the radio-frequency transitions between Stark 
sublevels particular rotational states. Grabner and Hughes (1950) obtained 
values for the quadrupole interaction energy the vibrational states 
inertia the zeroth vibrational state, and also the vibrational constant 
from intensity measurements several vibrational states. Schlier (1957) 
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obtained more accurate values for the quadrupole interaction using multipole 
version the electric resonance apparatus refocus and states 
and also obtained values for the spin-rotation interaction constants 
for both the and nuclei. 

the present experiment, through the observation 
transitions, have obtained the molecular rotational constants and 
and the electric dipole moment and also obtained quadrupole interaction 
constants agreement with those Grabner and Hughes and Schlier. 


THEORY 

The rotational spectrum obtained this experiment analyzed means 
Dunham’s expressions for the energy levels vibrating rotator. The 
frequencies the rotational transitions the first four vibrational states 
the absence hyperfine structure (h.f.s.) are given terms the Dunham 
Dunham (1932) and are restated Lew al. (1958). order 
use these expressions, the h.f.s. the state must taken into 
account. The major contribution h.f.s. due the quadrupole 
moment the nucleus which has spin 3/2; fluorine, with spin 1/2, 
has quadrupole moment. The h.f.s. pattern due the quadrupole 
was determined with the aid the table Casimir’s function given 
Townes and Schawlow (1955, Appendix 517) and with the aid the 
experimental results Schlier (1957) Grabner and Hughes (1950). First- 
order separations were found 1.5 Mc/sec; second-order energies 
were case greater than 0.1 kc/sec and were quite negligible (Townes 
and Schawlow 1955, Appendix II, 517). Further h.f.s. arises from the 
spin-rotation interactions for both nuclei. Schlier’s values for these 
interaction constants are 


+10.55+0.07 
0.380+0.097 kc/sec. 


Neglecting the potassium interaction, find for the h.f.s. the 
state the pattern shown Fig. The magnitude the splitting has 
been calculated from eq. Nierenberg, Rabi, and Slotnick (1948). This 
shows that the maximum splitting occurs line and approximately 
12.7 was expected that such splitting would observed and 
although splitting did appear some traces was not reproducible con- 
sistently and could not attributed unambiguously the interaction. 

The electric dipole moment may determined observing the change 
transition frequencies application static voltage the field. The 
most suitable line for this measurement line (Fig. 1). From eq. 10-31 
Townes and Schawlow (1955, 261), the 1/2 level neither 
split nor shifted weak field. Thus the only change the transition fre- 
quency arises from shift the level. This given (eq. 10-19, 
Townes and Schawlow) 
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Using Grabner and Hughes’ estimate seen that fields 
produce shifts less than 0.5 Mc/sec and satisfy the 
condition that 

eqQ. 


APPARATUS 

The apparatus has been described Lew (1958). The previously 
used FP-54 d-c. amplifier was replaced Applied Physics Corporation 
vibrating reed electrometer (model 31) using input resistor 
and taking the electrometer output via integrating circuit Leeds and 
Northrup 20-millivolt pen recorder. This system gave over-all sensitivity 
inch per ampere although was usually employed 1/3 1/10 
this sensitivity. The surface ionization detector was specially treated 
tungsten wire* which the potassium contamination had been largely 
removed. The background ion current from the special wire, though still 
appreciable, was much less than that from commercial tungsten wires. Under 
our experimental conditions the best signal-to-noise ratio was obtained with 
oxidized tungsten wire temperature about 1150° The absolute 
detection efficiency this temperature not known. After oxidation the 
efficiency was found decrease slowly over the course several hours. Hence 
capillary supply air the filament was provided permit periodic 
re-oxidation without loss vacuum. 
wire was kindly supplied Mr. Kilroy, General Electric Company, Cleveland, 

10. 
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attempt was made improve the signal-to-noise ratio chopping 
the beam cycles/sec and using tuned a-c. amplifier, integrator, and 
recording system. This had very poor sensitivity except high-filament 
temperatures where filament noise was prohibitive. The poor efficiency was 
attributed the filament time constant being greater than the chopping 
period and the inverse dependence this time constant the filament 


temperature. 


ELECTRONIC EQUIPMENT 

Two modifications were made the arrangement described Lew al. 
(1958). The first was the replacement the Varian X-12 klystron Sperry 
2K39 for measurements the transitions. This change was necessary 
because the Varian klystron had large bandwidth the higher frequency 
end its range. Since the Sperry klystron had frequency range 7,500 
10,300 was necessary feed the output into crystal rectifier 
and use the second harmonic for the inducement transitions. The second 
modification was the conversion the frequency multiplier chain give 
output 540 Mc/sec instead 270 Mc/sec. The sequence multiplications 
then The harmonic and associated 
markers could then used for klystron stabilization and frequency measure- 
ment rather than the much weaker 62nd harmonic 270 Mc/sec. 

The klystron control remained the same that used Lew but 
since recording system was used, the tuning condenser the frequency 
modulation receiver was driven motor rather than being 
set manually previously. 


EXPERIMENT 

The production beam molecules presented difficulties: 
thermocouple produced ion current the detector ampere. 
Commercial was used; load several grams filled the oven 
the bottom the slit and was sufficient for month’s work. 

The refocussed beam was found times greater than the predicted 
fraction part the main beam; part this accounted for states 
other than being refocussed. Under the best conditions the refocussed 
beam was several times the magnitude the residual beam, the residual beam 
being defined that arriving the detector with the stop-wire position 
and the field but the field off, the refocussed beam being defined 
the increase detected beam switching the field. 

The time constant the detection system introduced delay between 
transition and its recording. This resulted difference between apparent 
transition frequencies measured opposite directions scan. the rate 
scan was sometimes high per sec, this difference could amount 
kc/sec more. Consequently each quoted transition frequency the 
result averaging several scans each direction and then taking the mean 


the two averages. 
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Errors attached quoted frequencies cover the range measured values 
and are such give self-consistency for quantities such calculated 
from several frequency differences. 


RESULTS AND DISCUSSION 


Observed Lines, Quadrupole Interactions 

Rotational transitions belonging the three vibrational states 
were observed. These are listed Table Typical recordings several 
these are shown Fig. These are all for the line 3/2) 


ALL FREQ, Mc/sec 


LINE WIDTHS -080 - +100 


16573-943 16573-949 


-020 
16713-086 


For given vibrational state the successive peaks shown are repeated observations 
the same line, the arrows indicating the direction increasing frequency scan. The fre- 
quencies shown are not corrected for departure about parts 107 the Mc/sec 
standard from its nominal value. The correction made the values tabulated Table 


for the three observed vibrational states. can seen how the signal-to-noise 
ratio decreases rapidly the measurements the higher vibrational states. 
The line width increases the same time, probably because greater microwave 
power was required for the weaker lines. The .006-Mc/sec splitting shown 
one the lines was initially ascribed the interaction but 
attempts reproduce consistently were not successful. similar but larger 
splitting was observed line (Fig. 1). Since the splittings were not repro- 
ducible, they would not appear due interaction. Their origin 
not known. correction made for the effect finite field caused 
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contact potential shown below that this amounts most 0.4 kc/sec 
and therefore within the accuracy our 


TABLE 


F, of 


Line level v=0 v=1 


The quadrupole interaction constants measured the rotational 
state are given Table II. The values obtained Schlier and Grabner 
and Hughes are shown for comparison. seen that good agreement exists 
the cases and where our results are obtained from all three 


TABLE 


—7.744+0.040 
—7.658+0. 


lines, and (Fig. 1). the case where our result obtained 
from the and agreement not good but our result still 
overlaps that Schlier within the large uncertainty quoted. 


Rotational Constants 

The positions the rotational levels are determined the absence 
h.f.s. the use Fig. the case where there considerable 
uncertainty our value egQ, Schlier’s value was used conjunction with 
our measurement the stronger line The difference between the value 
obtained and that obtained using our value only kc/sec. The rota- 
tional transition frequencies unperturbed h.f.s. are then shown 


Table III. 
TABLE III 


Calculated rotational 
quencies the absence h.f.s. 


v=0 
v=1 


The rotational constants calculated from these frequencies are shown 
Table IV. For these calculations, use was made the values 
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(1953) from ultraviolet absorption spectra. assuming probable errors 
Mc/sec. Thus the quoted error arises almost entirely 
from the errors the vibrational constants. 


TABLE 
Rotational constants (Mc/sec) 


Yo = 

Yu = 69.991+0.010 
0.204+0.004 


Internuclear Distance 
using the relation may obtain the value the moment 


inertia 


order obtain the equilibrium internuclear distance r,, usual 
assume that the reduced mass obtained from the atomic weights the 
atoms concerned, may substituted the equation The value 
obtain this way 


where the error arises almost entirely from those Planck’s constant and 
the atomic mass unit 1955). Although this calculation neglects 
several electronic effects the molecule the above value probably the 
best one obtainable the absence exact knowledge the molecular 
wave functions. The electronic effects neglected are: 

Ionic Nature the Molecule 

The reduced mass ionic molecule cannot obtained from the neutral 
atomic masses the usual way because the displacement the valency 
electron. calculate the extreme case this effect, where the valency 
electron considered the fluorine nucleus, find that the change 
with the experimental error 

(it) Electron Slip during Molecular Rotation 

This may occur both valency electrons and electrons bound specific 
nuclei. Valency electrons having spherical distribution about the molecular 
center mass will undisturbed molecular rotation. Closed shells 
electrons having spherical distribution about nucleus will follow the 
nucleus during molecular rotation but their orientation will remain fixed 
space. However, ionic molecule, polarization will distort such spherical 
shells and the slip will reduced. 

have calculated the contribution the moment inertia the slip 
bound electrons use eq. 8-34 Townes and Schawlow. assuming 
that polarization exists and that all electrons are distributed spherical 
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shells about their respective nuclei are able find the maximum contri- 
bution this type slip. The equilibrium internuclear distance then becomes 


which differs from the previous value part 10*. Because the approxi- 
mations used the above treatment, great reliance can attached 
this value and the previous value obtained directly from will used 
the remaining discussion. 

higher than that calculated Honig the basis empirical relation 
which they had found accurate all other cases alkali halides. 
This unexpectedly large discrepancy led considerable waste time 
the search for transitions, did for Barrett and Mandel (1954) with 
the high-temperature microwave spectrometer. 

Pauling (1956) has published treatment interatomic distances alkali 
halide gas molecules which simpler than that Honig using the 
Born expression for the repulsive potential between two ions and neglecting 
any polarization effects, Pauling has derived simple formula for terms 
ionic radii and which had been previously obtained for crystals 


(Pauling 1948): 


where the Born exponent repulsive potential and 
the factor representing the ratio distances the vapor and crystalline 
phases (see Honig a/. 1954, 638). then assumes that which varies 
associated with the positive ion and with the negative. Apparently 
using out the experimental values which were known the time 
his writing, was able assign values the for all the alkalis and 
the for all the halides. With these was then able calculate values 
for the other seven molecules which turned out better agreement 
with experiment than Honig’s; this not surprising considering the number 
standards used. Pauling’s calculated values are shown Table along 
with the experimentally observed values. For comparison, values calculated 
Honig al. and Krasnov (1959) are also included. his standards, 
Pauling had presumably chosen the five bromides and CsF, KCI, and for 
were able reproduce his calculated values most closely using these 
standards. The and used are not identical with Pauling’s but 
believed that the latter had rounded his numbers off the first place after 
the decimal. The values used are 4.004, 4.690, 4.656, 4.520, 4.189 
for Lit, Nat, K+, and respectively and 2.073, 3.949, 4.500, 
5.299 for F-, and respectively. should pointed out that the 
and are not uniquely determined the experimental data; other sets 


are permissible. 
From Table can seen that Pauling’s calculated internuclear distances 
for the fluorides are not very good agreement with the experimental values. 
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TABLE 


Comparison between experimental and calculated equilibrium internuclear distances 
alkali halide gas molecules (in units cm) 


Present authors 


using 

LiF 1.51+0.08 1.527 1.518 1.527 1.584 
NaF 1.840 1.833 1.856 1.917 
2.129 2.130 2.154 2.166 
RbF 2.265 2.246 2.239 2.259 
CsF 2.345 2.346 2.370 
2.022 2.027 2.037 
NaCl 2.361 2.356 2.353 
2.787 2.793 2.790 
CsCl 2.906 2.911 2.916 
LiBr 2.170 2.166 2.170(s) 2.170(s) 
NaBr 2.502 2.500 2.503 
KBr 2.821 2.815 2.821(s) 2.821(s) 
RbBr 2.945 2.944 2.948 
CsBr 3.072 3.072(s) 3.075 
Lil 2.392 2.392 2.392 2.385 
Nal 2.712 2.727 2.719 
RbI 3.182 3.174 3.180 

3.315 3.306 3.315 


CsI 


Note: The quantities with (s) after them are those used as standards. 


| 


instead CsF use RbF the standard for the fluorides get the 
values given the fifth column Table The average agreement with 
experiment for the fluorides improved while for all the other alkali halides 
the agreement unchanged. More recently, Krasnov (1959) has adapted the 
empirical relation Honig and obtained the best agreement far 
between calculated and experimental values the interatomic distances 
alkali halide gas molecules. 
uses the equation 


where and are the polarizabilities the two ions, the Born exponent, 
and and are longer the crystal ionic radii but are the ionic radii 
free ions. set values and are obtained using certain experimental 
values standards much the same way Pauling obtained his sets 
and Once the and are known then the for the remaining mole- 
cules may calculated. These values are shown Table The improved 
agreement, particularly for the fluorides immediately noticeable, the values 
predicted for and agreeing within cm. therefore 
believed that Krasnov’s results would the most reliable for predicting 
transition frequencies the remaining alkali halides (LiF, NaF, and 


Dipole Moment 
was mentioned earlier, the line was chosen for the measurement 


the electric dipole moment the state. Static voltages volts 
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were applied the field which the gap was 0.69065 cm. The shift 
the frequency the line was measured with the static voltage first one 
direction and then the other. For +12 volts, 0.122+0.002 Mc/sec, 
for —12 volts, 0.103+0.002 Mc/sec. The difference between the 
two frequency shifts caused the contact potential the field. For 
this particular measurement the contact potential was 0.51+0.09 volt. The 


This value may compared with that e.s.u. obtained 
Grabner and Hughes and calculated from Rittner’s 
(1951) expression. 

The contact potential the field variable quantity dependent 
upon the immediate history the vacuum. Its highest measured value was 
0.7 volt. This will affect the energy levels nominally zero field. The 
level shows the strongest dependence the field and seen that for 
field corresponding contact potential 0.7 volt, the 
shift this level 0.00038 This considerably smaller than frequency 
measurement errors this experiment. 

Although the values the individual quantities and (or B,) given 
Grabner and Hughes have been shown considerable error, 
interest note that the value they get for the composite quantity 
excellent agreement with ours. They quote 
h?/2B, Hence, ignoring the stated uncertainties for 
reasons given below the Grabner and Hughes’ value 


while find (this work) 


7394 


This agreement not surprising. The frequency the transition they 
observed expresssible the form (Hughes 1949) 


where the strength the electric field which the transition takes 


Thus seen that the departure the observed frequency from quadratic 
greater accuracy than Since immediately yields value for while 
knowledge required for the extraction the individual quantities 
and B,, seen why the above agreements and discrepancies have occurred. 
For this reason too, the uncertainties quoted Grabner and Hughes for the 
individual quantities cannot carried over the quantity 


place: 
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substitute the accurate value found the present work into 
the quantity found Grabner and Hughes, find 


which agrees with our value very well. 

The dipole moments most the alkali fluorides and the alkali chlorides 
are now known. These are shown Table (following Kusch and Hughes 
(1959), but with some later entries). should noted that for the fluorides 
the dipole moment reaches peak RbF one goes down the column 
from Cs. The dipole moment CsF distinctly smaller than that 
RbF KF. This does not seem the case CsCl compared with 


TABLE 
Permanent electric dipole moments alkali fluorides and chlorides Debyes 
(for vibrational state 


Braunstein and Marple and 
Trischka (1955) Trischka (1956) 
8.5+0.2 Tate and 
Strandburg 
Present work Lee al. (1953) 
7.875+0.006 (1956) Trischka (1956) 


Relationships between the Kotational and Constants 
the Alkali Halides 


(1) 
assuming that the potential energy diatomic molecule may 
described Morse function, Pekeris (1934) obtained the relationship 


e 


We 


(2) 
this relation substitute values and w,x, tabulated Honig 
al. (1954) find for the values indicated crosses Fig. com- 
parison the experimental values are also plotted and are indicated 
seen that the experimental points fall very nearly straight line which 
satisfies the relation (with the variables megacycles per second) 


(3) 


the case the calculated points, however, least part the scatter must 

vibrational states once the first transition had been found, since 
simple theory rotational transitions 


fy. 1-0 = 2B. = [B.—a,(v 4)). 


a 
| 
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400 
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100 


uw 
oO 


plot vs. for the alkali halides. Except for RbF and KF, values are taken 
from Honig The values and used the calculated values are also from 


and 

seen from Fig. that the graph log versus log also roughly 
linear form but with considerable scatter. Expressing the variables 
megacycles per second, the plot may represented the relation 


the figure again can seen that the point deviates from the 
mean line much larger amount than any the other points. The possi- 
bility error suggested. examination the experimental data 
Honig shows that, with their errors from 0.2 1.0 sec the 
measured transition frequencies, the error the value can much 
larger than stated the authors. Whether the accurately linear dependence 
log log and the similar though less accurate one between 
and log has any fundamental significance not not known. 


4 
er 
uit vA 
35 x 
CsF 
x 
NoBr? 2.8! 
a x 
CsCl 
| 
x 
Cal | 


494 CANADIAN JOURNAL PHYSICS. VOL. 38, 1960 


1000 

800 

600 

400 
200 
80 NoBr?? F 
60} 


wd 


0.5 | e 2s 45 10 20 3040 
(Kmc/sec 


Fic. plot vs. for the alkali halides. Except for RbF and KF, values are taken 
from Honig al. 
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SPIN-LATTICE RELAXATION EFFECTS OBSERVED THE 
CONTINUOUS POWER SATURATION PARAMAGNETIC 


MARR AND PREM 


ABSTRACT 


The dependence the conventional saturation parameter the incident 
microwave power considered for Lorentzian-shaped paramagnetic lines and 
experimental observations may explained the basis 
transition probability which depends interactions. Values 
the effective spin-lattice relaxation times are compared with pulse technique 
determinations and estimates the corresponding phonon relaxation times are 
also given. 


With the advent the solid-state maser suggested Bloembergen (1956) 
and successfully demonstrated Scovil, Feher, and Seidel (1957) and others 
using magnetically dilute paramagnetic crystals, the need for clearer under- 
standing the processes involved the transfer energy from the para- 
magnetic ion the lattice required for the development other suitable 
maser materials. particular, the processes involved which determine the 
spin-lattice relaxation times are important. The experiments Bloembergen, 
Purcell, and Pound (1948) show that the relaxation nuclear magnetic states 
can adequately explained the Bloch equations (1946) with the aid 
spin-spin and relaxation times. The success these experiments, 
which the onset non-linearity saturation between the incident power 
and the absorbed power measured, led the application similar techniques 
electron resonance and Weidner 1953; Feher 
and Scovil 1957; Giordmaine, Alsop, Nash, and Townes 1958). 

alternative approach the problem has also been successtully emploved 
(Davis, Strandberg, and 1958), which the paramagnetic sample 
subjected high-power pulse radiation, saturate the sample, 
and the regrowth the absorption followed means small monitoring 
signal. The exponential rise the peak the absorption line allows value 
the relaxation time determined. Davis have pointed 
out that the value obtained their experiment decreases with increase 
monitoring power and that their low power values are much larger than 
continuous power saturation experiments under similar conditions. 
earlier paper, Strandberg (1958) indicated that may depend upon the 
degree saturation present phonons play significant part the relaxation 
process. 
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The purpose the present paper report the results measurements 
performed some paramagnetic resonance lines over range continuous 
power levels and show how the results may interpreted the basis 
transition probability which changes with the incident power the 


paramagnetic crystal. 


THEORY 
The magnetic properties paramagnetic substance can described 
terms the complex magnetic susceptibility, where the imaginary part, 
responsible for the absorption the microwave energy the medium. 


where the static susceptibility, g(v) the shape factor the line concerned, 
the saturation term, and the frequency peak absorption. Studies the 


paramagnetic spectra 


Swarup 1959) have shown that for small paramagnetic ion concentrations 
line shape Lorentzian. The shape factor may thus written 


the 


where the spin-spin relaxation time and related the width the 


insaturated line half-maximum absorption 
The saturation term defined 


and refer the ion populations the lower and upper 
states respectively and and are the equilibrium values the absence 


microwave radiation. 

the nuclear resonance experiments, the assumption made 
Bloembergen, and Pound that the may all times 
considered thermal equilibrium with the bath surrounding the 


these conditions, the peak value the differential 


ery low power and when decreased some traction its 


| 
| 
i 
4 
| 
‘43> 2 
LOW powe! Vi 


MARR AND SWARUP: SPIN-LATTICE RELAXATION EFFECTS 497 


While the above assumption applicable for nuclear magnetic resonance, 
not necessarily valid electron paramagnetic resonance. For this reason, the 
saturation term should derived terms which not invoke the concept 
spin and lattice temperatures initial assumption. 

The rate change the population difference under the influence micro- 
wave power given 


where the transition probability from the upper ionic state the lattice 
and the transition probability for the reverse process. the transition 
probability for the stimulated emission and absorption microwave radiation. 
The spontaneous transition probability can neglected that several 
orders magnitude less than this region the spectrum. Under steady- 
state conditions, the rate change population zero, thus enabling 
determined. standard radiation theory, the value for for the 


transition given 
(7) 
Substituting for Zg(v) equation (1) and differentiating w.r.t. obtain 


where and are the thermal equilibrium values the spin-lattice 


21 


Ie 


examination equation shows that least the early stages satura- 
tion, any significant change Must from within the square 
bracket and any suspected change the spin-lattice relaxation time with 
saturating power will indicated change itself since may 
written 

Van Vieck (1940) considered the spin the paramagnetic ion coupled 
the lattice through the mixing second- and higher-order perturbations 
the spin-orbit and coupling. Such mixing gives rise 
matrix elements which depend the parameters the spin and 
state the and addition, upon the phonon excitation expressed 


terms phonon quantum number Under these conditions, 
written the form 


(10) Wie 


Where the transition per unit phonon excitation, 


the mean phonon quantum number tor phonons having characteristic 
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The microwave frequency kMc/sec used the experiments well away 
from the upper limit for phonon frequencies imposed the interatomic 
distances the lattice and consequently p(v) can represented the Debye 
approximation 


(11) 


where and are the longitudinal and transverse phonon velocities respec- 
tively. the absence microwave radiation, energy exchanged between 
phonons due crystal imperfections, the coupling via the paramagnetic ion 
spins, etc., and the values tend fluctuate about their respective equilibrium 
values The probability phonon mode having energy 
governed distribution, and consequently the thermal 
equilibrium value (7) given 


No) I kT 
Under the influx the microwave radiation, energy transferred from the 
spins the lattice and the rate change can written (Strandberg 
1958) 


(13) 


where the excitation width for phonons direct contact with the spins 
and probably the same order less than the width the resonance line 
itself. this paper assumed that the two widths are equal; consequently 
the same symbol (Av) used for both widths. represents over-all relaxa- 
tion time for phonons way and interactions. 

For the steady state saturation conditions required here, zero 
and the departure from its thermal equilibrium value 


No) 
(14) 
and becomes 
FEU 


the initial stages saturation, value for may obtained assuming 

the spin populations not differ appreciably from their thermal equilibrium 
values, then substituting for may written 


saturation may therefore written 
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where 
1.2 
(18) K Ap(v){no) 27 ToT, 
and 


The difference between equations and that may not zero for 
paramagnetic resonance. 


EXPERIMENTAL 


conventional X-band paramagnetic resonance spectrometer utilizing 
impedance bridge was employed. Phase-sensitive detection and 
amplification the signal allowed the resonance line 
derivative Rectiriter chart recorder. The cavity operated the 
mode situated between the pole pieces 6-in. Varian magnet suitable 
low-temperature Dewars. The cavity and crystal were immersed liquid 
helium 4.2° order obtain good coupling between the crystal and 
the low-temperature bath possible. 1000 modulation the magnetic 
field was kept below peak-to-peak value 0.5 oersted. The steady magnetic 
field was swept the rate oersteds/minute and was controlled the 
driving motor the recorder. Pound stabilized X-13 klystron was used for 
the low-power measurements and Sperry for high-power saturation 
measurements. The power entering the cavity was monitored measuring 
the power Polytechnic 650-A universal power bridge. The 
power incident the cavity was variable over the range 
watts. 

conjunction with the observations other workers (Meyer 1955), 
initial experiment showed that the width half-maximum absorption 
powdered sample the free radical diphenyl did not change 
with and that the peak value the derivative, increased linearly with 
increase over the range power available. The under were 
all less volume and were placed the cavity along with small 
sample the powdered Values could then obtained for both 
the paramagnetic and under identical power conditions. 
Ratios for the paramagnetic crystal for made over the 
range available, enabled the saturation the paramagnetic crystal 
These ratios, normalized unity low power levels are repre- 
sentative equation Measurements the width the paramagnetic 
line low power enabled determined. The relation between and 
was obtained trom 


‘ 
(20) Hy (2 I) 


taken from (1955) for the cavity. the loaded 
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the cavity volume and the reflection coefficient. For the experi- 
ments reported here, 3000 and 0.66 giving 1.6 where 

The crystals studied were and and 
were prepared mixing known proportions weight the diamagnetic salt 
with its isomorphous paramagnetic salt, dissolving the mixture distilled 
water, and growing the combined crystals slow evaporation. this means, 
crystals were obtained where the nominal percentage the diamagnetic ions 
replaced paramagnetic ions were 0.1%, 0.5%, and for the ion and 
for the ion. determination the actual percentage made 
the case the crystals Atomic Energy Canada Ltd. X-ray 
fluorescence showed that the actual concentrations were somewhat less than 
the nominal concentrations which are used designate the crystals. 


IV. RESULTS 
The normalized peak absorption derivative ratios obtained for different 
values the microwave power incident the cavity are shown Fig. for 


100 


‘The saturation paramagnetic resonance lines 4.2° 
Experimental points for 1.0% 
Fitted theoretical curves 


the +3) transitions for the ion and the ion 
4.2° The full curves are plots equation with suitable values and 
constant value for was indicated from the experimental measurements 
with slight departure for the 0.5% ion the region 107! 
watt range. This departure can accounted for the change 
and The best fits the experimental data were obtained with 
and 5254-50 for the 0.5% ion and values 

Attempts were made extend the measurements the ion 
which contained and nominal paramagnetic ion concentrations, 
the former case, the very low powers needed for saturation and the decrease 
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signal-to-noise ratio reduced the accuracy obtainable. the case the 
crystal, extension the measurements higher values required 
determine accurately. 

Table shows the experimental quantities obtained 4.2° The values 
are the inverse the spin-lattice transition probabilities under thermal 
equilibrium conditions Estimates the relaxation 
times, are obtained from the values assuming that the populations the 
states concerned had not changed appreciably from thermal equilibrium values. 
Estimates these populations were made with the aid the actual concentra- 
tions available. 


TABLE 


Experimental results 4.2° 


AH, 


(.4) 20.8 0.6X10 140 


*Actual concentrations shown in parentheses. 


The change effective with saturation shown more clearly over 
limited power range Fig. For equal zero, and consequently constant 


in m wott 


value straight line would obtained. While this may justified 
ion curves. 

absolute values the relaxation times depend upon large number 
variables; however, they are probably correct within order magnitude. 
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0-5% Cr (7, in msec) 


\ 


effective T, 


0-00! 0-01 10 100 


in m watt 


Fic. effective with microwave power. 

tinuous power saturation experiments. 

Dotted curve: values obtained Davis al. (1958) for 


The change observed Davis al. for crystal containing nominal 
concentration 0.5% using their pulse 
technique and varying the monitoring power, shown Fig. with the power 
level arbitrarily placed the logarithmic scale (since the absolute power 
levels were not available). Also shown are the variations the effective 
values for the 0.5% and the obtained from the data Table 
appears likely that the same change being observed the two experiments. 
The absolute values depend rather critically upon the concentration 
the paramagnetic ion the crystal. Table indicates that decreases 
concentration increases. 

Since p(v) and are computable possible obtain estimate 
the spin-lattice transition probability per unit phonon excitation which 
should depend only the parameters the spin and electronic state the 
ion. These estimates are shown Table and appears that the concentration 
ions also affects the value The relaxation times 4.2° 
also appear change with concentration shown the table; whether they 
are fact independent the paramagnetic ion concerned, remains 
seen more crystals are examined. The results obtained are agreement 
with those suggested Strandberg (1958) less than seconds. 

Attempts were made observe saturation the liquid air 
temperatures but insufficient power was available obtain very reliable 
results. They did, however, indicate that decreases about The 
relaxation times were not determinable these temperatures. 
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NOTES 


Y-RADIATION AIR-GROUND INTERFACES WITH DISTRIBUTED 
SOURCES* 


For problems involving the y-ray dose due distributed sources impor- 
tant take account the scattering and absorption the 
interface. recent work the Defence Research Chemical Laboratories 
(Clifford al. 1959) measurements have been made that show the effects 
different materials the dose varying heights and radial dis- 
tances from the ground. was used the source since the 
y-ray energy (0.66 Mev) close the average energy y-rays from fission 

Some the more important results these measurements are given and 
also interest compare these results with calculations Berger (1957) 
for somewhat similar conditions but for higher photon energy 1.28 Mev, 
corresponding y-rays. 

Figure shows the results measurements made over smooth clay earth 
and compares the results Berger’s Monte Carlo calculations for 1.28 Mev 


6 
x K=90 

10 
x 


BASE DISTANCE METERS 


Fic. Boundary correction factor air above smooth clay earth. measured values 
with air-ground interface; ——from calculations for 


which ground and air were assumed have the same scattering properties 

per unit mass water. The parameter which appears the factor which 

the dose altered the presence the denser medium. Since both Co™ and 

y-rays are the energy region Compton scatter for low and medium 
*Issued D.R.C.L. No. 296A. 


address: Department Physics, McGill University, Montreal, Que. 
address: Ontario Cancer Institute, Toronto, 
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materials, might anticipated that values calculated for 1.28 Mev 
would little changed 0.66 Mev when distances are expressed mean 
free paths. energies below the region predominant Compton scatter, 
experimental values would expected vary from those calculated 
higher energies. Good experimental agreement with the calculations was 
reported Titus (1958) 1.28 Mev for source the plane boundary 
between media steel wool and steel. The agreement shown Fig. for the 
lower energy y-rays from indicates the accuracy with which 
values could used fall-out analysis. 

Similar values, obtained from measurements made region polyfoam 
(i.e. expanded polystyrene) density 0.072 g/cc with source located 
boundary, are shown Fig. Here the experimental 
values shown were obtained taking the ratio the dose measured the 
indicated base distance and height the dose observed equivalent 
distance when the source and detector were embedded polyfoam feet 
above the boundary. 

Because Monte Carlo calculations were available (Berger 1957) giving 
values for the limiting case source lying the surface completely 


HEIGHT 


100 200 300 
BASE DISTANCE 


HEIGHT— 


250 350 
BASE DISTANCE 


Fic. Boundary correction factor polyfoam above concrete. measured values 
with concrete-polyfoam interface; Berger’s calculations for Co®. 

Boundary correction factor above strongly absorbing medium. measured 
values with lead-foam interface; from calculations for 
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absorbing medium, the polyfoam was used obtain experimental values. 
This was done stacking the polyfoam above floor lead bricks and 
measuring the dose different distances and heights the foam when 
source was placed the surface between the lead and polyfoam. Since the 
ground now absorbing medium, the dose should something less than 
the dose infinite polyfoam medium. The results are shown Fig. where 
the lines represent the Monte Carlo calculations for Co® source in- 
finitely absorbing medium. The somewhat higher experimental values are 
probably due the use source lower energy y’s than and the 
extent which the lead floor not perfect absorber. The trend the 
values quite similar and further illustrates the reliability Monte Carlo 
calculations. 

The spectra the collimated radiation received various angles were 
measured and analyzed give the dose received through solid angle 0.1 
steradian function the direction angle collimated detector located 
meter above interface uniformly contaminated with The 
resulting polar distribution radiation shown Fig. where the doses 


4 
oN 
1 
/ 
4 
yak 


ANGLE FROM VERTICAL DEGREES 


Fic. distribution radiation. Dose rate per 0.1 steradian height meter 
above plane uniformly contaminated with 


calculated from the spectra are shown crosses. The solid line between angles 
100° and 180° measured from the vertical was calculated for primary radia- 
tion only since detector height meter will see only small amounts 
scattered radiation arriving these angular limits. When integrated over the 
steradians this polar distribution should lead the total dose actually 
observed height meter uniformly contaminated field. The dotted 
line was chosen satisfy this requirement. interesting note that the 
resulting curve indicates that 73% the total dose due radiation received 
within +20° the horizon. These measurements were made with contamina- 
tion smooth plane. the practical case fall-out appreciated that 
the polar distribution radiation would modified somewhat depending 
the extent ground roughness. 
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NOTE DIFFRACTION HALF PLANE 


WILLIAMS 


INTRODUCTION 


The problem the diffraction the field line source perfectly 
conducting half plane was first solved Carslaw (1899) using method similar 
that employed Sommerfeld (1896) considering diffraction plane 
wave. The problem was also considered Macdonald (1915) using com- 
pletely different approach and, particular, succeeded expressing the 
solution extremely simple manner. This latter form the solution was 
recently rederived Clemmow (1950) using yet another method. Though 
the form Macdonald’s solution extremely simple the methods used 
obtaining involve considerable amount analysis and the object the 
present note give simple method solving this problem involving 
very little analysis. For plane-wave excitation comparable type simple 
solution has been obtained Lamb (1906). 


FORMULATION THE PROBLEM 

Cartesian co-ordinate system Oxyz chosen that the half plane occu- 
source parallel the axis and its intersection with the plane will 
taken the point or, polar co-ordinates, Two other 
geometrical quantities which will required are the distances and 
(x,y) from and respectively. terms polar co-ordinates 


are defined 

The only component the electric field will parallel the axis and 
will denoted exp Maxwell’s equations show that 
(1) 
where the two-dimensional Laplacian operator, function, 
and denote the electric susceptibility and magnetic 
permeability (in m.k.s. units) free space. 
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finite except vanishes the half plane, and satisfies 
Sommerfeld’s radiation condition infinity. also have that satisfies 
the reciprocity condition, i.e. solving 
diffraction problem necessary make some assumption concerning the 
behavior the solution near sharp edges diffracting obstacles. the 
present case shall make the assumption, following Bouwkamp (1954), that 
behaves near All known solutions diffraction problems in- 
volving thin sheets do, fact, exhibit this behavior near the edge the 
may shown that the smallest permissible value shall make 
the assumption that, near has the form where function 

From can construct new function which also vanishes the half 
plane and satisfies the radiation condition. This function defined 
and is, fact, solution the homogeneous Helmholtz equation. 
Our assumption about the edge behavior shows that, near has 
the form sin Thus possible solution for sin 
exp This solution satisfies the homogeneous Helmholtz equation, the 
radiation condition, the boundary condition the half plane and possesses 
the appropriate edge singularity. The most general solution for consistent 
with our assumption about the edge behavior thus obtained adding the 
above particular solution second solution Helmholtz’ equation which 
non-singular the edge. This second solution also finite everywhere 
and satisfy Sommerfeld’s radiation condition. The uniqueness theorem 
Peters and Stoker (1954) shows that such function everywhere zero and 
hence sin exp —ikr. The form shows that satisfies the 
reciprocity condition and thus must have the form sin exp 
where numerical constant. 

The boundary value problem thus reduces integrating the partial differen- 
tial equation 


Ox 2+/(rro) 


exp 


not immediately obvious how this equation may integrated and the 
correct line attack obtained considering the particular case 
This corresponds determining solution Laplace’s equation with line 
charge the solution this problem may obtained conformal 
mapping (Jeans 1951). The solution given linear combination variables 
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The form and suggests writing the product sines equation 
difference cosines and writing where 


(3) “exp —ik(r+ro), 


The form the solution for suggests that choose and new 
variables and easily verified that 


The above relations indicate that function and and function 
equations and may integrated immediately give the following par- 
ticular integrals 


va b 


all the conditions imposed except that The behavior the source 
then the integral representation the Hankel function (Watson 1944) shows 

Thus the solution the boundary value problem has been obtained 
assuming the edge behavior the solution and employing the reciprocity 
condition. 

The above method may also employed determine time-dependent 
solutions the wave equation. easily verified that the two-dimensional 
time-dependent Green’s function obtained when the exponential equations 


where the wave velocity. The integration the equations elementary 
and the solution obtained agrees with that given Friedlander (1951). 
the exponential replaced 
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where the unit function, then the solution for cylindrical pulse ob- 
tained. The problem for the spherical pulse, solved Wait (1957), may also 
treated this manner. this case the exponential replaced 


The solution for time harmonic point source also easily obtained the 
present approach. 
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FRANCK-CONDON FACTORS AND r-CENTROIDS FOR SOME BANDS 
THE FOURTH POSITIVE BAND SYSTEM 


The factor arrays number band systems 
combustion interest discussed Gaydon (1957) were recently published 
(Nicholls, Fraser, and Jarmain 1959). The extent each the arrays was 
limited include only bands which were the readily accessible wavelength 
region (2000 10,000 standard photographic spectroscopy, and whose 
Values were amenable computation, desk calculator electronic 
computer, using the methods Fraser and Jarmain (Fraser and Jarmain 1953; 
Jarmain and Fraser 1953; Fraser 1954) for Morse molecules. The fact that 
the bands the fourth positive system (which whole extends 
wavelength between about 1000 and 2800 A), longer wavelengths than 
2000 arise from transitions between high vibrational levels, precluded the 
use the method for these bands although they are com- 
bustion interest. 

partial factor array was, however, com- 
puted some years ago for this system part routine program such 


*Now Yerkes Observatory, Williams Bay, Wisconsin. 
leave absence the National Bureau Standards, Washington, D.C. 
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calculations. gives data for bands whose wavelengths are less than 1988 
The recent growth interest intensity aspects vacuum ultraviolet 
spectra together with private requests for the array prompts present 
Table together with band head wavelengths and r-centroids. The wave- 


TABLE 


factors r-centroids and wavelengths for the fourth 
band system 


O.114 0.26: 0.285 0.196 0.096 0.03s 0.0lo (0.002)} 
0 1.181 1.211 1.241 1.271 1.302 1.334 1.366 1.398 
1545.3 1597.14 1653 .05 1712.19 1774.90 1841.47 1912.8 1988 .O8T 
0.216 0.155 0.003 0.059 0.195 0.185 (0.094) 
1 1.160 1.190 1.220 1.248 1.280 1.311 1.342 
1509.66 1560 .14* 1613 .99T 1669.68 1729.25 1792.38 1859.3 
0.230 0.01: 0.09% O. 116 0.005 (0.075) 
2 1.141 1.170 1.200 1.223 1.259 1.289 
1477.48 1525.75 1576.67 1630.40 1687 .44T 1747.20 
0.18) 0.020 O.116 0.000 (0.105) 
3 1.123 1.151 1.179 1.209 1.2388 
1447.27 1493.60 1542.34 1593 .97F 1647.90 
0.119 008s 0.034 (0.075) 
4 1.105 1.133 1.161 1.189 
1418.91 1463 .47* 1510.47¢ 1559.47 
0.069 0.12. (0.005) 
5 1.088 1.115 1.1438 
1391.97 1435.28 1480 
0.037 (0.115) 
6 1.072 1.099 
(O. O16) 
1.057 
1344.27 
LEGEND: 
(A) 


(A) 
*Bands observed overlapped, wavelengths thus calculated, not measured, by Read (1934). 
tBands not reported, wavelengths calculated. 
tFranck-Condon factors in parentheses not perturbation corrected. 


lengths were taken from the work Headrick and Fox and Read 
(1934). The r-centroids were evaluated graphically methods described 
Nicholls and Jarmain (1956), and the factors were evaluated 
the method Fraser and Jarmain (including perturbation correction) 
using the basic data for the transition tabulated Herzberg (1950). 

expected for band system which involves moderately large change 
internuclear separation (Ar, 0.107 A), the primary Condon 
the values fairly wide and does not quite include the (0,0) band. 
There also evidence secondary parabola whose vertex the region 
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the (2,2), (2,3), and (3,2) bands. Some recent studies one the 
analytical form the Condon cases where the potentials in- 
volved are not parabolic will published shortly. much more extensive 
r-centroid array for the fourth positive system than that presented 
Table available part comprehensive compilation such data for 
wide variety band systems, and the process publication elsewhere. 
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Under this heading brief reports important discoveries physics may published. These 
reports should not exceed 600 words and, for any issue, should submitted not later than six weeks 
previous the first day the month issue. proof will sent the authors. 


Frequency Measurement Standard Frequency 


Measurements are made Ottawa, Canada, using N.R.C. caesium-beam frequency resonator 
reference standard (with assumed frequency 192 631 770 
tions from nominal are quoted parts per negative sign indicates that the frequency 


below nominal. 


Date, 
December MSF, 
1959 60 ke/s average* 

1 —159 — 164 

2 N.M. —166 

3 N.M. —164 

4 N.M. — 163 

5 N.M. — 169 

6 N.M. N.M. 

7 NM. N.M. 

8 — 163 N.M. 

9 N.M. —164 

10 — 166 —160 

1l —165 —165 

12 —163 —164 

13 —173 — 164 

14 —170 — 160 

15 —172 —161 

16 —158 —160 

17 — 162 — 1638 

18 — 163 —161 

19 —165 — 163 

20 —165 —162 

21 N.M. —160 

22 N.M. — 162 

23 N.M. —160 

24 N.M. —164 

25 -171 —164 

26 N.M. — 163 

27 — 166 — 166 

28 —155 —161 

29 —163 —164 

30 —156 —161 

31 —162 —163 

Midmonthly 
mean —164 —163 
Midmonthly 


mean of —106 


Note: N.M. no measurement. 


*Time of observations: 00.00 to 07.00 and 22.30 to 24.00 U.T. 
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GBR, 16 ke/s 


24-hour 
average 
—I161 
— 166 
— 167 
— 166 
N.M. 
N.M. 
N.M. 
N.M. 
—162 
— 162 
— 164 
— 162 
— 164 
—161 
— 163 
—161 
—161 
—161 
— 162 
N.M. 
—161 
— 162 
—161 
— 1638 
— 164 
— 164 
— 163 
— 162 
— 161 
— 157 


—162 


KK2XEI, 


60 ke/s 
— 82 
— 80 
N.M. 
N.M. 
N.M. 
—78 
N.M. 
—79 
— 82 
N.M. 
N.M. 
—79 
—78 
—78 
N.M. 
N.M. 
—80 
—79 
—77 
—78 
N.M. 
N.M. 
N.M. 
—76 
—74 


—73 


Frequency devia- 
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